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ABSTRACT 


An  advanced  re-entry  systems  heat-transfer  handbook  for  hypersonic 
flight  has  been  developed  using  aerothermodynamic  prediction  methods  developed 
during  the  X-20A  (Dyna  Soar)  Program.  It  contains  l)  desigi  procedures  for 
computing  aerodynamic  heating  rates  to  re-entry  vehicle  configurational 
elements,  2)  discussion  of  differences  betveen  aerodynamic  heat  transfer 
and  pressure  distributions  observed  in  present  day  vind  tunnels  and  those 
vhich  would  occur  in  actual  free  fli^t,  3)  vind  tunnel  to  flight  extra¬ 
polation  factors,  k)  simplified  expressions  for  estimating  stagnation  point 
and  swept  cylinder  turbulent  stagnation  lire  heating  rates,  and  5)  graphs 
for  rapid  calculation  of  heating  rates  and  extrapolation  to  fli^t  factors. 

The  information  presented  is  applicable  to  complex  maneuverable  vehicles 
as  well  as  ballistic  bodies. 
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SYMBOLS  (continued) 
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SYMBOIB  (continued) 

transformed  heating  rate,  eq,.  (BUa) 
arbitrary  function  of  X,  eq.  (B13a) 
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transformed  velocity  component,  eq.  (BIO) 
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transformed  velocity,  eq.  (BIO) 
leakage  velocity 
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free  stream  velocity 
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Cartesian  coordinates;  curvilinear  coordinates  (see 
Appendix  B) 

heat- transfer  equivalent  distance 

surface  distance;  transformed  x  coordinate,  eq,.  (BIO) 

coordinate  normal  to  body  surface 

value  of  y  at  leading  edge  of  delta  wing 

transformed  y  coordinate,  eq^,  (BIO) 

coordinate  normal  to  x  and  y 

compressibility  factor,  ?/pWS}  transformed  z  coordinate, 
eq.  (BIO) 

angle  of  attack;  Prandtl  number  exponent,  eq.  (Bte) 
velocity  or  pressure- gradient  parameter,  eqs.  (3*6)  and  (b4i) 
ratic  of  specific  heats,  Cp/c^ 

vorticity- interaction  parameter,  eq.  (B31);  pressxu'e 
gradlem;  function,  eqs.  (B35)  and  (B36) 

flow  deflection  angle;  cone  semi -vertex  angle;  flap 
deflection  angle;  bouiidary- layer  thickness 

effective  wedge  ang^e 

boundary- layer  displacement  thickness 

total  streamline-divergence  parameter;  compressibility 
parameter,  eq.  (7.I);  finite  increment  operator 

shock  standoff  distance 

ficticious  length  (see  Appendix  C) 

emlsslvlty;  density  ratio  across  normal  shock;  surface  ray 
angle  measured  from  sharp  delta -wing  apex;  streamline  angle 
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SYMBOLS  (continued) 


extrapolation  factor 

similarity  parameter,  eq.*  (Blla) 

e 

boundary -layer  momentum  thickness,  eq.  (Bj);  angular 
coordinate 

e* 

leading  edge  streamline  angle 

shock  wave  angle 

0 

transformed  momentum  thickness,  eq.  (Blla) 

X 

canopy  angle 

A 

sweep  angle 

M 

absolute  viscosity 

V 

Prandtl-Meyer  function 

P 

density 

PrMr 

reference  density-viscosity  product  (see  Appendix 

b) 

a 

partial  Prandtl  number  for  translation,  rotation, 
vibration  only  (see  Appendix  B)j  Stefan-Boltzmaim 

and 

constant 

Or 

Prandtl  number  evaluated  at  enthalpy  and  pressure 
ponding  to 

corres- 

2 

boundary  layer  parameter,  eqs.  (B37)  and  (B38) 

T 

shear  stress;  model  skin  thickness 

7 

transformed  shear  stress,  eq.  (BU) 

^(nH-l)/ra^  eq.  (B13b) 

'I' 

’^w  angle 

CO 

coefficient  defined  In  eqe.  (2-12)  and  (2.13) 

Subscripts : 

ATM  ) 
atni,SL| 

evaluated  at  se^  level  condition 

av,AW 

adiabatic  wall 

BC 

blunt  cone 

BD 

blunt  delta  wing 

SYMBOLS  (continued) 


C 

c 

CYL 

DW 

e 

eff,EFP 

eq 

EXP 

F 

f 

FUT 

FP 

G 

HL 

1 

I-I 

L 

LE 

£1 

■ax^MAI 

aln^NHf 

N 

n 

NL 

o 


cone 

crossflow 
cylinder 
delta  wing 

at  boundary-layer  edge 

effective 

equivalent 

on  expansion  surface 
flap 

on  sphere  in  laminar  flow;  on  60*  swept  cylinder  in 
turbulent  flow 

at  flight  conditions 

flat  plate 

control  surface  gap 

at  hinge  line 

initial;  perfect  gas;  index  variable 

at  line  of  Intersection 

laminar^  leakage 

leading  edge 

at  X  ;  atean 
n 

maximum 

minimum 

normal 

in  normal  direction 
no  leakage 

tangential;  at  stagnation  conditions;  at  laminar  reference 
conditions 
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r 

BEF 

ref 

RAD 

BT 

s 

S 

S» 

SC 

SD 

SH 

SL 

sm^SM 

SP 

T 

vort 

w,W 

WT 

X 

e 

A 

0 

2 

to 

+ 


evaluated  at  enthalpy  and  pressure  corresponding  to 
conditions  used  to  evaloiate 

reference 

pertaining  to  reference  streamline 
radiated 

at  turbulent  reference  conditions 

streamwlse 

at  stagnation  point 

evaluated  at  stagnation  enthalpy  and  local  pressure 

sharp  cone 

sharp  delta  wing 

at  shoulder 

on  stagnation  line 

smooth 

stagnation  point 
turbulent;  total 
vorticlty  effects  included 
at  wall 

at  wl.nd  tunnel  conditions 

at  distance  x 

based  on  angulrr  location 

based  on  geometric  sweep  nnsrle 

at  Mach  number  zero 

behind  shook  wave 

free  stream  or  undisturbed 

Just  ahead  of  sudden  compression 
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STMBOIS  (concluded) 


Just  behind  sudden  cooipresslon 


Superscripts : 
* 


at  reference  ten^ierature  defined  by  eq[.  (3.10a) 


I.  IMTRODUCTION 

Successful  design  of  a  maneuverable  re-entry  vehicle  requires  detailed 
knowledge  of  aerodynamic  heating  rates  over  all  surfaces  of  the  vehicle. 
Performance  of  the  vehicle  is  limited  by  the  temperature  capability  of 
structural  materials  which  tend  to  establish  lower  limits  :)n  the  radii  of 
the  nose  and  leading  edges  and  upper  limits  on  the  angle  of  attack  at  which 
the  vehicle  can  operate.  Unfortunately,  design  features  that  are  favorable 
from  the  performance  standpoint  often  have  undesirable  heating  characteristics. 
Aerodynamic  heating  studies,  therefore,  are  required  to  determine  the  proper 
design  compromises. 

Aerodynamic  heating  rates  are  influenced  by  all  of  the  properties  of  the 
flow  field  in  addition  to  the  influence  of  body  surface  conditions.  In  the 
external  flow,  the  density-  velocity,  pressure  (both  level  and  gradient), 
streamline  pattern,  and  chemical  composition  must  all  be  determined  before  heat¬ 
ing  calculations  can  be  made.  At  the  body  surface  the  important  conditions  ere 
the  wall  temperature,  smoothness,  and  permeability.  The  boundary  layer  state, 
which  may  be  laminar,  transitional  or  turbulent ,  may  also  depend  upon  all  flew 
properties . 

The  ability  to  accurately  predict  Inflight  aerodynamic  heating  rates  varies 
with  geometry  and  flight  condition.  Analytic  methods  are  available  for  simple 
shapes  in  laminar  flowj  for  all  other  conditions  a  combination  of  theoretical 
and  empirical  methods  must  he  used.  A  ccanpletely  empirical  approach  is  precluded 
by  the  limitations  of  testing  facilities,  none  of  which  can  provide  completely 
satisfactory  simulation  of  the  re-entry  environment. 

This  report  presents  a  compilation  of  design  methods  for  computing  aero¬ 
dynamic  heating  rates  that  were  developed  during  the  X-20  (Dyna-Soar)  Program. 

The  results  of  more  than  five  years  of  aerodynamic  heat  transfer  testing  and 
analysis  are  reflected  in  these  methods.  Methods  developed  specifically  for  that 
program  are  presented  in  this  report  in  a  generalized  manner  Intended  to  supply 
information  useful  for  preliminary  design  of  advanced  re-entry  vehicles.  The 
methods  presented  were,  however,  developed  to  treat  the  critical  d':3ign  con¬ 
ditions  for  the  X-20,  and  are  not  intended  to  be  applied  at  velocities  greater 
than  26,000  .fps  nor  to  surfaces  cooled  by  mass  Injection.  The  flow  is  assumed 
to  behave  as  a  continuum,  ana  to  be  in  chemical  equilibrium.  With  the  exception 
of  vortlcity  Interaction  in  the  region  of  the  forward  stagnation  point,  low 
Reynoxds  number  effects  are  neglected.  No  mass  addition  effects  are  considered, 
but  some  information  is  given  regarding  leakage  which  may  occur  at  panel  Joints 
or  control  surface  hlngelinee.  The  effects  of  surface  roughness  were  found  to 
be  large  and  are  discussed.  Flow  field  interference  effects  on  X-20  fin  leading 
edge  and  canopy  heating  were  found  to  be  severe  and  an  spproxlmate  method  that 
defines  the  upper  limit  is  presented.  While  these  assumptions,  limits  and 
approximations  were  appropriate  and  sui’flcient  for  X-20  design  studies,  they 
should  be  reviewed  when  considering  other  lifting  re-entry  vehicles. 

The  material  that  is  presented  in  this  report  Is  divided  into  seven 
major  sections.  The  following  section  (Section  II)  describes  the  methods 
for  several  basic  two-dimensional  and  axisymmetrlc  shapes.  Combinations  of 
these  shapes  ere  discussed  in  Section  III,  and  surface  condition  effects  on 
the  level  of  basic  smooth  body  heating  rates  are  described  in  Section  IV. 
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Suggestions  pertaining  to  the  application  of  the  design  methods  to  com¬ 
plex  configurations  are  given  in  Section  V.  Tiiis  section  also  contains  two 
examples  from  the  X-20  progreim  that  illustrate  the  requirement  for  detailed 
analysis  tailored  to  the  re-entry  vehicle  being  studied  . 

Section  \T:  contains  charts  that  give  numerica],  results  for  several  basic 
geometric  shapes,  i.ncluding  the  hemisphere,  flat  plate,  swept  cylinder,  unyaw¬ 
ed  blunt  and  sharp  cones,  and  blunt  and  sharp  delta  wings  at  angle  of  attack. 
Heating  rates  for  all  configurations  are  given  by  distribution  cvirves  which 
are  to  be  multiplied  by  reference  heating  rates  to  give  absolute  values.  Ihia 

form  of  presentation  is  concise  and  simplifies  the  estimation  of  heating  rates 

for  geometries  not  specifically  presented.  Hemisphere  stagnation  point  heat¬ 
ing  rates  are  used  as  Isuninar  reference  values  for  all  shapes.  Turbulent 
heating  for  sharp  tip  configurations  are  referenced  to  the  turbulent  heating 
rate  for  a  60-degree  swept  infinite  cylinder.  For  blunt  bodies  it  was  found 
that  better  accuracy  was  obtained  by  referencing  turbulent  heating  rates  to 
laminai  hemisphere  stagnation  i>oint  values. 

Section  VII  discusses  interpretation  and  use  of  data  from  wind  tunnels 
and  other  ground-hased  test  facilities.  In  some  cases,  wind  tunnel  data  can 
be  applied  in  flight  by  merely  adjusting  the  level  hy  use  of  the  reference 

heating  rates  described  in  the  previous  paragraph.  However,  for  complex  geo¬ 

metries,  further  correction  is  often  necessary  becatise  of  real-gas  and  inter¬ 
acting  flow  field  effects. 

Appendices  are  included  that: 

1.  Give  simplified  expressions  that  approximate  some  of  the  more  ccxnplex 
heating  formulas. 

2.  Show  derivations  of  basic  methods. 

3.  Rrovide  a  method  for  calculating  laminar  boundary  layer  displacement 
thickness. 

4.  Give  examples  of  the  application  of  the  charts  presented  in  Sections 
VI  and  VII. 
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II.  METHODS  FOR  BASIC  SHAPES 

This  section  describes  the  methods  that  were  used  to  predict  heating 
rates  for  several  basic  shapes.  Including  flat  plate,  sharp  cone,  sphere, 
infinite  cylinder,  and  sharp  delta  wing.  These  geometries  were  selected 
because  relatively  detailed  and  accurate  calculations  can  be  made  for  each 
without  relying  heavily  on  experimental  results,  and  because  such  calculations 
provide  information  that  can  be  used  in  analyzing  more  practical  configurations. 

For  the  purposes  of  this  section,  heating  calculations  can  usually  be 
broken  Into  two  separate  steps:  l)  the  calculation  of  the  flow  properties 
neglecting  the  effect  of  viscosity,  and  2)  calculation  of  the  boundary  layer 
using  the  inviscid  flow  results.  The  inviscid  flow  calculations  must  pro¬ 
vide  a  cciuplete  description  of  conditions  at  the  boundary  layer  outer  edge, 
including  the  local  pressure,  velocity  (both  magnitude  and  direction),  tem¬ 
perature,  and  chemical  composition.  No  practical  completely  general  method 
has  been  developed  for  calculating  the  inviscid  flow  for  all  of  the  basic 
shapes.  It  was  therefore  necessary  to  use  several  different  methods  to 
ealcu3.ate  the  flow  about  the  basic  shax>es.  Including  normal  and  oblique  shock 
theory,  modified  Newtonian  theory,  Prandtl-iieyer  expansion  theory^  the  method 
of  characteristics,  and  the  blunt  body  inverse  method.  In  addition  to  these 
well  known  methods,  an  empirical,  pressure  distribution  was  used  for  infinite 
cylinders,  and  a  crossflow  method  was  used  for  predicting  delta  wing  stream¬ 
lines.  Except  where  otherwise  noted  real  gas  effects  are  Included  in  the 
flow  field  calcul-ations,  using  gas  properties  of  reference  1 .  Chemical  and 
themal  eqalllbrivun.  has  been  assxaned  in  all  calculations. 

In  contrast  to  the  variety  of  inviscid  flow  methods  used,  all  of  the 
boundary  layer  calculations  for  the  basic  shapes  were  made  with  a  single 
method,  the  "  Method".  This  method,  which  was  developed  by  R.  A. 

Hanks  in  the  course  of  the  X-20  program,  is  applicable  in  either  laminar 
or  turbulent  flow.  The  laminar  flow  version  of  the  p  j.Mr  is  based 

on  exact  similar  solutions  of  the  boundary  layer  equations.  The  exact 
solutions  were  used  to  construct  curves  of  boundary  layer  parameters,  in¬ 
cluding  form  factors,  (  8*/Q  ),  crossflow  momentum  thickness,  Reynolds 
analogy  factors,  and  effective  values  of  the  density -viscosity  product. 

The  boundary  layer  parameter  curves  are  then  applied  'in  nonsimilar  flows 
through  the  use  of  a  momentum  integral  equation. 

The  turbulent  flow  form  of  the  p  ^  fir  method  uses  the  same  momentiuu 
Integral  equation,  retaining  terns  that  are  often  neglected,  such  as  those 
related  to  the  effect  of  pressure  gradients  on  the  turbulent  boundajry  layer 
profiles.  Although  such  profile  changes  are  small,  and  cannot  be  exactly 
calculated,  their  effect  may  be  to  increase  heating  rates,  and  it  is  there¬ 
fore  unconservative  to  neglect  them  entirely.  The  evaluation  of  the 
turbulent  boundary  layer  parameters  was  guided  by  the  known  laminar  trends, 
end  in  some  cases  are  t<^ken  directly  from  the  aforementioned  laminar  curves. 

The  actual  values  used  were  selected  also  on  the  basis  of  agreement  with 
e:q>erlmental  results. 


3 


since  the  P  r  v  emhodlee  essentially  ell  laminar  similarity 

solutions  as  well  as  turbulent  flow  parameters^  it  is  necessarily  somewhat 
complex,  and  a  complete  description  is  given  in  Appendix  B.  For  the  follow¬ 
ing  discussion  of  methods,  it  is  sufficient  to  note  that  in  addition  to 
boundary  layer  edge  pressure,  density,  etc.,  which  are  required  by  all 
boundary  layer  methods,  the  P  j.  My  method  also  requires  knowledge  of  three 
flow  divergence  parameters,  r,  f,  and  A  .  The  parameter  r  accounts  for 
the  effects  of  streamline  divergence  on  heating  due  to  body  geometry,  which  is 
most  simply  evenpiified  by  the  Increase  in  heating  due  to  streamline  divergence 
on  an  tinyaved  cone.  The  parameter  f  accounts  for  the  effect  of  streamline 
divergence  due  to  transverse  pressure  gradients,  and  so  is  affected  by  body 
shape  only  Indirectly.  The  total  streamline  divergence  is  denoted  by  A  , 
which  Is  related  to  r  and  f  by 

A=  constant  *  r  •  f  (2.1) 


where  the  constant  is  arbitrary.  Referring  to  Appendix  B  it  will  be  seen 
that  r,  f,  and  A  can  always  be  normalized  so  that  absolute  values  are 
never  required.  Hence  it  is  sufficient  to  determine  the  ratios  r/r±, 
and  A  /  A  1  where  the  subscript  i  denotes  some  initial  or  reference  value. 
In  actual  practice  equation  (2.1)  is  usually  applied  In  tha  form: 


A  _  r  •  f 


(2.1a) 


The  distinction  between  r,  f,  and  A  is  necessary,  since  the  ultimate  effect  of 
streamline  divergence  on  heat  transfer  depends  on  the  presence  or  absence  of 
transverse  pressure  gradients. 

The  P  r  V.  method  has  been  used  for  all  heat  transfer  calculations 
presented  in  this  report  except  where  otherwise  noted.  In  particular, 
laminar  boundary  layer  calculations  involving  interactions  with  the 
external  flow  were  calculated  with  the  Boeing  Nonslmilar  Boundary  Layer 
Program,  which  is  described  in  Appendix  C.  This  program  numerically  inte¬ 
grates  the  boundary  layer  equations  in  partial  differential  form,  and  so  is 
free  of  similarity  restrictions.  This  program  was  required  for  only  three 
calciilatlons :  the  expansloi:  side  of  the  deflected  flap,  the  vortlclty  Inter¬ 
actions  at  the  stagnation  point,  and  surface  wave  calculations  described  In 
Section  IV.  The  particular  methods  used  for  each  of  the  basic  shapes  are  des¬ 
cribed  below. 

A .  Flat  Plate 

Exact  calculations  of  Inviscid  flow  properties  for  the  sharp-leading-edge 
flat  plate  were  made  using  oblique-shock  theory  and  real -gas  properties.  A 
description  of  real-gaa  oblique-shock  theory  Is  given  In  reference  2.  The 
actual  values  used  in  this  report  were  recalculated  using  gas  properties 


presented  in  reference  1.  The  three-dimensional  flov  i)arameter8  r,  f,  and  A 
are  constants  for  this  case,  so  that  r/r^,  f/f^  and  A/A^^  are  all  unity. 

B.  Hemisphere 

The  pressture  in  the  subsonic  portion  of  the  hemisphere  flov  field  vas 
calculated  by  Modified  Newtonian  theory.  In  the  Newtoniein  theory,  the  free 
stream  is  regarded  as  an  aggregation  of  particles  that  impinges  on  the  sur¬ 
face  of  the  body,  and  in  so  doing  changes  its  momentum  normal  to  the  surface 
into  pressure  acting  on  the  body.  This  idealization  of  the  flov  field  leads 
to  the  folloving  expression  for  the  pressure  coefficient 

Cp  *  2  sin^8  (2.2) 

vhere  S  is  the  true  flow  deflection  angle.  The  same  equation  is  obtained 
if  one  considers  the  impact  of  a  supersonic  fluid  stream  on  an  inclined 
surface,  provided  the  shock  wave  is  close  to  the  surface.  This  latter 
condition  is  approximately  satisfied  in  hypersonic  fli^t,  for  which  the 
density  ratio  P^/p^  is  large,  since  the  shock  standoff  distance,  A  , 
is  related  to  tne  density  ratio  by 

As  »  [o]  R 

R  »  Hemisphere  radius  (2.3) 

as  shown  in  many  references.  See,  for  example,  references  3  aiKi 

Equation  (2.2)  is  exact  only  if  A  *  0.  For  a  finite  density  ratio 
and  hence  finite  A  ,  improved  accuracy“can  be  obtained  with  Modified  New¬ 
tonian  theory 

Cp  =  K  sin^S  (2.h) 

where  K  is  a  constant  which  is  expected  to  be  near  2.0.  Its  actual  value 
for  the  hemisphere  can  be  calculated  by  normal  shock  theory  and  isentropic 
compression  to  zero  velocity.  The  result  is  very  well  approximated  by 

K  -  2  -  pjp^  (2.5) 

Equation  (2.5),  which  is  obtained  by  assuming  that  the  density  is  con¬ 
stant  on  the  stagnating  streamline  downstream  of  the  bow  shock,  is  derived 
in  Chapter  7  of  Truitt  (reference  5),  among  other  places. 

The  pressure  distribution  in  the  supersonic  flow  portion  of  the  hemisphere 
was  calculated  by  Prandtl-Meyer  expansion  theory  matched  to  values  given  by 
equation  (2.4)  in  the  vicinity  of  the  sonic  point.  The  actual  matching  point 
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was  selected  to  malce  both  the  pressure  distribution  and  its  derl  native 
continuous.  The  resulting  pressure  distribution  is  shown  in  figure  2''!. 

Velocity  distribution  can  be  calculated  fros  the  pressure  distribution 
either  by  Isentroplc  expansion  or  by  integrating  the  Bernoulli  equation 
equation 


^  =  1  dp 

d8  pg  Ug  de 


(2.6) 


where  u 

e 


and  6  are  defined  in  the  sketch  below 


p  la  the  boundary  layer  edge  density  and  P  ie  the  local  pressure.  Because 
e 

the  pressure  gradient  and  velocity  are  zero  at  the  stagnation  point  the 
Bernoulli  equation  must  he  applied  in  the  special  form 


in  s  P 

»x  R  "V’  Pe 


(2.7) 


where  x  is  measured  fr'om  the  stagnation  point  and  jB  is  the  velocity  gradient 
parameter.  Using  the  Hewtonlan  theory  pressure  distribution  and  the 
equation  of  state^  equation  (2.7/  leads  to 

^  V2  Z)o 

OD 

where  ^  is  the  gas  constant,  T  is  the  real-gas  temperature,  and  Z  is  the 
compresslhillty  factor. 

As  pointed  out  in  reference  3  the  use  of  modified  Newtonian  theory  for 
the  hemisphere  pressure  distribution  has  been  substantiated  empirically. 

The  Prandtl-Meyer  matching  method  has  also  been  verified  in  this  manner. 
However;,  it  is  not  necessary  to  only  rely  upon  experimental  data  for  their 
substantiation.  Comparisons  have  been  made  with  inverse -method  calculations 
(reference  6)  and  with  the  integral  method  calculations  of  Belotserkovskli 
(reference  7)  that  show  that  the  method  described  here  is  accunite  to 
within  a  few  percent. 

Streamline  divergence  on  a  hemisphere  is  due  entirely  to  body  geometry, 
and  is  given  by 


(ii.9) 
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where  A  is  the  local  streamline  divergence  parameter,  is  the  streamline 
divergence  i>arameter  at  8  a  0",  and  8  is  the  ti^,e  flow  deflection  angle. 


The  basic  hemisphere  heating  rates  for  both  laminar  and  turbulent  flow 

were  calculated  by  the  p  /*  method.  At  very  low  Reynolds  nxanbers  and  high 

r  r 

Mach  niimbers  a  correction  to  the  laminar  heating  rates  must  be  made  to  account 
for  the  effects  of  vorticity  in  the  external  flow,  which  may  become  important 
when  the  boundary  layer  is  not  thin  In  comparison  to  the  shock  layer  thickness. 
Analysis  and  experimental  data  have  been  published  (e.g.  references  8  -  10}  that 
show  Increases  in  heating  rate  resu.\tlng  from  vorticity  interaction. 

There  is  considerable  disagreement  between  the  various  published  theor¬ 
etical  predictions  of  the  vorticity  effect.  Many  of  the  predictions  contain 
rather  drastic  simplifications  and/or  patching  of  separate  approximate 
solutions  for  various  regions  within  the  shock  layer.  The  Nonsimilar 
Boundary  La;yer  Program  was  used  to  obtain  more  complete  numerical  solutions 
to  the  boundary  layer  equations  which  included  the  effects  of  external 
vorticity. 


A  comparison  between  results  obtained  using  the  Nonsimilar  Boundary 
Layer  Program  and  those  shown  in  reference  8  is  shown  in  figare  2-2  in 
terms  of  the  stagnation  point  vorticity  interact! on  parameter  F ,  defined  in 
reference  8  as; 


r  = 


(2.10) 


where  Re^ls  a  stagnation  region  Reynolds  number  defined  as 


Re 


0 


Op 

^o 


(2.11) 


and  0  ,  n  ,  and  i  are  I'espectively  the  density,  viscosity,  and  enthalpy 
0.0  0 

evaluated  at  stagnation  conditions,  llie  quantities  A  and  w  are  coefficients 
in  an  equation  for  the  inviscid  tangential  velocity; 


1  X 

=  a|  +  (0  (y/R)  +  . . .] 


(2.12) 


with 


A  = 


R  9u 


and 


(0 


where  x  and  y  are  curvilinear  coordinates  tangential  and  normal  to  the  body, 
R  is  the  body  radius,  and  u  is  the  velocity  normalized  with  respect  to  the 
boundary  layer  edge  velocity. 
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Physically  the  vorticity  parameter  F  represents  the  ratio  of  the  shear 
stress  in  the  shock  layer  to  the  shear  stress  at  the  body  surface,  i.e., 


r~ 


t 


Pe 

"r=  0 


where  is  the  potential-flow  velocity,  the  subscript  e  denotes  evaluation 
at  the  edge  of  the  boundary  layer,  and  is  the  vorticity-free  skin- 


friction  coefficient.  Numerical  values  of  F  can  be  obtained  if  u  and  «  are 
known.  The  potential-flow  velocity  u  can  be  obtained  from  the  Newoonlan  pres¬ 
sure  distribution  and  the  Bemo\illi  elation  as  defined  by  equation  (2.8).* 

The  vorticity  parameter  F  was  computed  using  the  following  equation  for  w  which' 
is  given 

r  -  .  (1-0  T/  .1/2  1 

(2.13) 


where  c  is  the  density  ratio  across  the  normal  shock  and  A  /R  is  the  dimen- 

3 

sionless  shock  standoff  distance. 


As  shown  in  figure  2-2,  Nonslmllar  Boundary  Layer  Program  calculations  are 
somewhat  lower  than  those  by  Ferrl,  et  al.  in  reference  8.  At  the  time  calcu¬ 
lations  were  made  by  The  Boeing  Company,  there  was  considerable  controversy 
over  the  magnitude  of  the  vorticity  interaction  effects,  with  the  majority  of 
the  published  methods  also  falling  well  below  the  reference  8  prediction. 

The  controversy  could  not  be  resolved  by  experimental  results,  due  to  the  ex¬ 
treme  difficulty  in  obtaining  reliable  data.  The  Nonsimilar  Boundary  Layer 
Program  calculations  Involved  fewer  simplifying  assumptions  than  those  in  re¬ 
ference  8,  which  could  be  the  source  of  the  difference  shown.  However,  an  ex¬ 
tensive  review  of  the  assumptions  involved  in  the  reference  8  analysis  was 
made,  using  the  detailed  Information  provided  by  the  Nonslmllar  Boundary  Layer 
Program,  without  discovering  any  serious  discrepancy. 

Van  Dyke  (reference  ll)  and  others  have  pointed  out  that  vorticity  inter¬ 
action  is  only  one  of  several  second-order  boundary  layer  effects,  and  suggest 
the  proper  consideration  of  these  additional  effects  would  reduce  Ferri's  (re¬ 
ference  8)  predictions  somewhat.  This  observation  does  not  explain  the  lower 
values  given  by  the  Nonsimilar  Boundary  Layer  Program,  however,  since  only  the 
classical  boundary  layer  equations  were  used.  The  effect  shown  is  entirely 
due  to  the  boundary  condition  applied  at  the  outer  edge  of  the  boundary  layer, 
namely, 


Constant 


*Note  that  p of  equation  (2.8)  is  very  nearly  equal  to  the  parameter  A  of 
equation  (2.10)  for  hypersonic  flight,  since  i^  ~ 

■iH^These  results  ere  not  conveniently  displayed  in  the  format  of  figure  2-2. 
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rather  than 


the  condition  applied  In  no-vorticity  calculations. 

The  present  authors  fesl  that  the  Nonsimilar  Boundary  Layer  Program  re¬ 
sult,  being  e.  more  complete  calculation,  is  less  likely  to  contain  any  hidden 
errors,  and  recommend  the  use  of  the  Nonsimilar  Boundary  Layer  Program  curve 
of  figure  d-2.  However,  once  P  is  calculated  the  reader  can  easily  obtain 
numerical  results  for  both  methods. 

C .  Infinite  Cylinder 

The  pressvire  in  the  immediate  vicinity  of  the  stagnation  line  of  an  un¬ 
swept  infinite  cylinder  was  obtained  by  modified  Newtonian  theory,  equation 
(2.^).  An  empirical  cxurve,  shown  in  figure  2-1,  was  used  for  the  pressure 
distribution  away  fiom  the  stagnation  line.  It  was  observed  that  there  was  very 
little  difference  be\';ween  the  data  from  many  different  tests;  hence  a  single 
curve  was  used  for  all  calculations. 

The  velocity  distribution  on  the  unawept  cylinder  can  be  calculated  using 
the  same  methods  as  were  used  for  the  hemisphere.  The  actual  velocity  dis¬ 
tribution  differs  slightly  from  that  on  a  hemisphere,  however,  since  the  two 
pressure  distributions  are  not  identical. 

The  invisoid  flow  over  a  swept  cylinder  can  be  calculated  by  resolving 
the  flow  into  components  normal  and  parallel  to  the  cylinder  axis.  The  flow 
normal  to  the  cylinder  can  then  be  calculated  by  the  same  method  used  for 
unswept  cylinders.  The  pressure,  temperature,  density,  and  circumferential 
velocity  depend  only  on  that  component  of  the  free  stream  flow  normal  to  the 
cylinder.  The  parallel  flow  component  is  constant.  The  streamline  pattern 
can  be  calculated  from  the  velocity  components  by  integrating  the  relation¬ 
ship 


y(x) 


(2.14) 


where  y  is  the  distance  of  the  streamline  from  the  stagnation  line 
(measured  along  the  cylinder  surface)  and  x  is  the  distance  measured  along 
the  cylinder  stagnation  line.  The  body  geometry  streamline  divergence 
parameter  r,  due  to  geometry  is  equal  to  the  cylinder  radius.  Since  r  is 
constant  for  all  x,  r/r^  =  1  and  there  is  no  direct  effect  of  r  on  the 

heating  rate.  There  is  of  course  a  pressure  gradient  effect  on  heating 
through  the  parameter  f.  The  total  streamline  divergence  is  then  equal  to 
that  due  to  pressure  gradients  alone,  so  that  equation  (2.1a)  reduces  to 


^  =  -L 

^  fi  (2.15) 
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where  d  and  f  can  both  be  obtained  from  the  integral  of  equation  (2.14). 
D.  Sharp  Unyawed  Cone 

Shock-wave  angles  and  flow-field  properties  for  a  sharp  unyawed  cone 
were  obtained  using  the  method-of -characteristics,  reference  12.  Since 
transverse-pressure  gradients  are  zero  for  an  unyawed  cone  the  value  of 
f/f^  =  1.  Hence,  the  expression  for  the  streamline  divergence  parameter 
reduces  to 

JL  =  A  _  iL 
ri  ~ 


E.  Sharp  Delta  Wing 


The  delta  wing  presents  a  more  diffic\tlt  analytic  problem  than  any  of 
the  other  basic  shapes  discussed  so  far,  and  rigorous  theoretical  methds 
are  available  for  only  a  few  restricted  cases.  For  sharp  leading  edge  delta 
wings  with  the  flow  superaonlc  everywhere  it  can  be  shown  that  all  flow 
properties  must  be  constant  along  any  straight  line  throu^  the  aipex.  At  low 
angles  of  attack  oblique  shock  theory  can  be  applied  to  the  region  upstream 
of  a  Mach  line  on  the  lower  surface  throu^  the  apex.  Oblique-shock  theory 
is  applied  in  the  plane  normal  to  the  leading  edge  at  an  effective  free- 
stream  Mach  number,  M  ,  and  wedge  angle  which  can  be  obtained  from  purely 
geometric  calculations  as  follows 


6'j' 


(2.17) 


and 

Mj^  =  V 1  -  sin^  A  cos^  a  (2.17a) 

where  a  is  the  angle  of  attack  end  A  is  the  geometric  sweep  angle. 

It  is  seen  that  for  very  highly  swept  delta  wings,  the  effective  wedge 
angle  will  be  much  larger  than  the  angle  of  attack,  so  that  shock  detachment 
will  occur  at  relatively  low  angles  of  attack.  For  very  high  angles  of 
attack  the  flew  at  the  leading  edge  is  outward,  (i.e.,  with  the  leading  edge 
behaving  as  a  trailing  edge),  and  the  flow  field  can  be  calculated  by  an 
extension  of  the  blunt  body  method  of  integral  relations  (reference  j)  as 
shown  in  reference  13 .  Such  calculations  have  thus  far  been  made  only  for 
ideal  gases  with  constant  specific  heat  ratios. 

For  most  of  the  conditions  of  practical  interest,  however,  neither  of 
the  above  methods  could  be  applied  and  approximate  methods  become  necessary. 
It  was  still  desirable  that  the  methods  be  as  valid  as  possible  in  the 
actual  flight  environment.  Accordingly,  emphasis  was  placed  on  developing 
approximate  analytic  methods  for  predicting  delta  wing  pressure,  velocity, 
and  streamline  patterns. 
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1.  Sharp  Delta  Wing  Pressure  Method 

The  approximate  analytic  method  devised  for  predicting  delta  wing  pressures 
was  a  modified  Newtonian  method.  In  which  the  constant  K  of  equation  (2.4)  was 
allowed  to  vary  so  as  to  obtain  the  best  agreement  with  known  exact  solutions. 
This  was  done  by  first  plotting  wedge,  cone,  and  blunt -body  stagnation  point 
theoretical  results  as  shown  in  figure  2-3* 

It  is  seen  in  figure  2-3  that  for  wedges  and  cones  the  coefficient  K  at 
first  decreases  with  normal  Mach  n\imber,  and  then  Increases  rapidly  to  a  limit 
fixed  by  shock  detachment.  The  effect  is  seen  to  be  less  pronounced,  and  to 
occur  at  larger  deflection  angles  for  cones  than  for  wedges.  For  both  wedges 
and  cones,  K  is  eilways  greater  than  2.0.  Bl\mt-body  stagnation  point  vedues 
for  K  are  seen  to  behave  in  an  opposite  manner,  being  always  less  than  2.0,  but 
increasing  slowly  with  normal  Mach  number.  (This  behavior  is  consistent  with 
the  previous  expression  for  K  given  for  blunt  bodies,  equation  (2.5))*  For 
delta  wings,  it  is  to  be  expected  that  the  wedge  curves  will  apply  in  the  lead¬ 
ing  edge  region  at  low  angles  of  attack,  as  previously  discussed.  At  very  large 
angles  of  attack,  approaching  90  degi*ees,  it  is  to  be  expected  that  the  blunt 
body  curve  would  be  more  applicable.  In  the  intermediate  range  K  would  be  ex¬ 
pected  to  vai^r  more  or  lees  smoothly  between  the  wedge  and  blunt  body  curves. 


Shock  detachment  is  caused  by  the  Inability  of  attached  flows  to 
satisfy  mass  conservation  requirements.  For  the  sharp-delta  wing,  the 
phenomenon  of  shock  detachment  would  not  be  expected  ..to  occur  in  the  seme 
sense  as  for  a  wedge  or  unyawed  cone  since  mass  in  conserved  by  flow  outwartl 
over  the  leading  edge.  Tlie  existence  of  outflow  at  the  delta  wing  leading 
edge  at  very  high  angles  of  attack  also  implies  that  at  some  Intermediate 
angle  the  flow  must  be  Just  parallel,  to  the  leading  edge,  and  so  would 
correspond  closely  to  cone  flow.  The  curve  of  K  for  the  delta  wing  would 
therefore  be  expected  to  begin  at  the  wedge  curve,  fair  throu^  the  cone 
curve,  and  terminate  on  the  blunt -body  curve. 


Such  a  family  of  cvtrves  was  drawn,  and  is  presented  in  figure  2-4  which 
was  used  for  all  delta  wing  calculations  in  this  report.  The  upper  branch  is 
based  on  a  formula  given  by  hypersonic  small- disturbance  theory, 


(2.18) 


which  was  used  primarily  because  it  gives  a  curve  of  the  desired  shape.  The 
curve  shown  was  calculated  for  V  »  1.1,  a  value  which  was  selected  because  it 
leads  to  the  desired  curve,  and  not  because  7  »  1.1  is  considered  realistic. 
(a  value  of  1.1  is  clearly  unrealistic  in  comparison  to  wind  tunnel  data.) 

The  lower  branch  is  faired  so  that  K  approaches  a  value  of  2.0  as  the  Mach 
number  approaches  infinity,  varying  as  follows 


K 


(2.19) 


for  Mach  numbers  greater  than  10 
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Althouc^  the  fairings  of  figure  2-k  are  sooewhat  arbitrary,  the  develoineDt 
just  given  is  preferable  to  a  purely  empirical  approach,  since  all  the 
required  exact  solutions  can  be  made  for  flight  conditions. 

The  major  objective  of  figure  2-k  is  to  predict  delta  wing  pressures; 
hovever,  it  is  also  expected  to  provide  reasonable  predictions  for  almost 
any  shape,  since  it  agrees  well  with  both  slender  and  bl\mt-body  results.* 

2.  Velocity  and  Streamline  Pattern 

The  Invlscld  velocity  component  parallel  to  the  delta  wing  centerline 
was  calculated  by  the  following  expression: 


—  =  1  -  ^2/5600 
u 


(2.20) 


where  a  is  the  wing  angle  of  attack  in  degrees.  Equation  (2.20)  was 
obtained  in  the  same  manner  as  the  pressure  curve  described  previously  by  an 
examination  of  cone  and  wedge  theory.  For  two-dJmenslonal  nonreacting  flow, 
the  Invlscid  velocity  depends  only  on  the  shock  wave  angle  and  density  ratio, 
and  is  given  by;** 


(2.21) 


where  -  is  the  shock  wave  ang].e  in  degrees.  The  wing  angle  of  attack 
is  found  by  the  following  equation 


08  -  O  = 


(2.22) 


*Figure  2-k  can  not  be  expected  to  apply  to  interfering  flow  fields,  however, 
since  one  of  the  effects  of  interference  can  be  to  Increase  the  local  dy¬ 
namic  pressure. 

♦•Equations  (2. 21)  and  (2.22)  are  obtained  by  resolving  the  flow  into  compon¬ 
ents  normal  and  parallel  to  the  shook.  In  crossing  the  shock  the  normal 
velocity  component  is  reduced  in  proportion  to  the  reciprocal  of  the  density 
ratio.  The  parallel  component  is  unaffected. 
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Some  typical  results  for  various  values  of  the  density  ratio  are  shovn  below: 


Equation  (2 .20)  Is  also  shown  for  comparison. 

As  shown  in  the  above  sketch  the  velocity  increases  with  density  ratio  at 
constant  O!  to  an  upper  limit  which  equation  (2. 21)  shows  tc  be  a  cosine  function. 
For  cones  no  rigorous,  simple  formula  is  available,  since  the  velocity  varies 
as  the  flow  compresses  within  the  shock  layer.  However,  tl  e  velocity  changes 
during  this  compression  are  slight,  and  equation  (2 .21)  is  nearly  correct  for 
cones  as  well  as  wedges.  Tlie  value  of!  ^a~a  for  cone  flow  is  only  about  one- 
half  the  wedge  value  given  by  equation  (2.22)*,  however,  so  that  the  cone  curves 
vo\ild  be  expected  to  lie  about  midway  between  the  wedge  and  cosine  curves. 

It  is  seen  that  the  velocity  ratios  are  fixed  within  fairly  nairow  limits 
for  density  ratios  greater  than  about  5.  In  hypersonic  flight  density  ratios 
less  than  5  can  occur  ouLy  at  small  wedge  deflection  angles  (less  than  about  11 
degrees  at  Mach  20)  where  the  effect  of  density  ratio  is  seen  to  be  small. 
Equation  (2 .20)  was  therefore  chosen  as  a  close  approximation  for  hypersonic 
fll^t  for  angles  of  attack  up  to  40  degrees  or  more. 


*  !Ihl8  result  is  obtained  by  applying  a  mass  balance  to  the  cone  shock 
layer,  assuming  that  8urf8u:e  velocities  are  equal  to  the  wedge  values. 
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The  spanwlse  velocity  component  for  sharp  delta  wings  is  obtained  by  a 
crossflow  theory.  As  shown  in  the  following  sketch,  in  the  crossflow  plane 
the  delta  wing  appears  as  a  sharp-cornered  blunt  body  for  which  the  velocity 
at  the  comer  must  be  equal  to  the  local  sonic  velocity,  a*  (the  detailed 

derivation  of  reference  13  leads  to  the  same  conclusion) . 


The  local  sonic  velocity  a^ 
appearing  in  the  crossflow  plane. 


is  calculated  from  the  Mach  number  component 
and  is  given  by 


* 

hi 


a 


00 


2  +  (y  -  1)  (Mo,  sino)^ 
.  2  +  (y^ 


(2.23) 


At  the  centerline  the  spanwlse  velocity  component  is  zero  but  Increases 
with  distance  away  from  the  centerline.  Superimposing  the  spanwlse  compon¬ 
ent  onto  the  axial  flow  component  (equation  (2.20))  leads  to  a  streamline 
qualitatively  similar  to  that  sketched  below. 


The  leading  edge  streamline  angle,  6^*,  with  respect  to  the  centerline  is 
calculated  from  the  foHowing  equation: 

tane*-^  (2 -S'*) 

n  ue 

using  equation  (2.20)  for  the  axial  velocity  component  Uq, 


Ik 


It  will  be  noted  that  this  method  predicts  outflow  even  at  zero  angle  of  attack. 
This  inconsistency  was  avoided  by  modifying  equation  (2.23)  to: 


>i< 

a„  =  a 


M  sin  a 


(2.25) 


This  modification  is  equivalent  to  neglecting  the  effect  of  the  free  stream 
static  temperat'ore  on  and  was  used  throughout  this  report. 

The  distribution  of  the  spanwise  velocity  component  cannot  be  calculated 
except  for  the  few  previously  noted  special  cases  for  which  complete  flow  field 
ceilculatlons  can  be  made.  However,  the  published  calculations  for  the  flat** 
faced  blunt  body  show  that  the  assumption  of  a  llneeu:  distribution  is  conser¬ 
vative,  i.e.,  overpredicts  the  spanwise  component.  Delta  wing  experimental  data 
were  similar,  but  more  nearly  linear  than  the  bl\iat  body  theory  would  indicate. 

A  linear  distribution  has  therefore  been  assumed  in  the  prex>aration  of  this 
report . 


-  \  ^y/yMAx^ 


(2.26) 


where  y  is  the  spanwise  distance  measured  normal  to  and  from  the  centerline 
and  yj^  is  the  leading-edge  value  of  y.  The  streamline  equation  was  obtained 

using  equations  (2.20)  and  (2.26)  and  is  given  by: 


y  =  C 

where  C  is  a  constant  for  any  particular  streamline,  and 

n  ^ 

■=  tan  6*  tan  A 


n  = 


Uetan(90'’  -  A) 


(2.27) 


(2.28) 


The  streamline  patterns  corresponding  to  equation  (2.27)  fall  into  two 
classes,  depending  on  the  veOLue  of  n.  For  n  less  than  1.0  the  streamlines 
originate  at  the  leading  edge;  for  n  greater  than  1.0  all  streamlines  origin¬ 
ate  at  the  apex.  The  two  types  of  patterns  are  sketched  below: 


n  <  1 


n  >  1 
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Since  the  streenllne  divergence  parameter  A  Is  proportional  to  y,  it  can 
be  expressed  using  equation  (2 .27) -as  _ _ 

&/A,  =  (2.29) 

'^1  +  (dy/dxp  i 

Since  the  delta  vlng  is  taken  to  be  flat,  it  is  assumed  that  there  is  no 
three-dimensional  effect  due  to  body  shape,  so  that  f/f^  “A/A^. 

3.  Heat  Transfer  Calcoilatlons 

The  sharp  delta  wing  heating  curves  of  Section  VI  were  calculated  using 
equation  (2.29)  together  with  pressure  given  by  figure  2~k,  the  velocity  com¬ 
ponents  gl«en  by  ecruatlon  (2 .20)  and  (2.26),  and  the  p  u  method  of  Appendix 

r  r 

B.  For  the  centerline  of  the  sharp  delta  wing  the  three-dimensional  flow  effect 
can  be  calculated  in  closed  form  for  laminar  flow  and  is  given  by 


h/hpp  =  V 1  +  2  n 


(2.30) 


where  n  is  defined  by  equation  (2.28),  and  hpp  is  the  heat  transfer 
coefficient  for  a  flat  plate  evaluated  using  delta  wing  surface  properties. 


A  similar  expression  can  be  derived  for  turbulent  flow  using  the  Blaslus 
skin-friction  lav: 


C 


f  " 


.0592 

Re*2 


(2.31) 


The  resulting  expression  for  the  three-dimensional  effect  is: 


Equation  (2.32)  is  an  approximation  in  that  equation  (2.3I)  is  considered 
to  be  less  accurate  than  the  Karraan-Schoenherr  or  Praadtl-Schlichting  equations. 
However,  as  shown  in  Appendix  B,  the  error  introduced  into  equation  (2.32)  is 
small. 

These  simple  forms  (equations  (2.30)  and  (2.32))  are  possible  only  because 
the  streamwise  distance  to  a  point  on  the  centerline  is  unaffected  by  streamline 
divergence.  For  locations  not  on  the  centerline  the  streamwise  distance  depends 
on  the  streamline  pattern,  as  may  be  seen  frcsa  the  previous  sketches.  The  com¬ 
bined  distance  and  flow  divergence  effects  must  be  evaluated  by  integrating  the 
boundary-layer  equations  along  the  streamline,  as  was  done  in  the  preparation  of 
the  design  curves,  figures  6-32  and  6=33 • 
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Figure  2 


Figure  2-2:  STAGNATION  POINT  VORTICITY  INTERACTION  EFFECT 


Figure  2-y  WEDGE,  CONE,  AND  STAGNATION  POINT  PROPORTIONALITY  CONSTANT 


Figure  2-4  BOEING  MODIFIED  NEWTONIAN  HYPERSONIC  PRESSURE  COEFFICIENTS 


III.  COMBINATIONS  OF  BASIC  SHAPES 

This  section  describes  methods  that  were  used  to  calculete  the  aero- 
thermodynamic  characteristics  of  three  combinations  of  basic  shapes,  the 
hemisphere-cone,  the  blunt  delta  wing,  and  a  flat  plate  with  a  trailing  edge 
flap.  The  blunt  delta  wing  is  made  up  of  a  flat  delta  wing  with  cylindrical 
leading  edges  and  a  spherical  cap.  The  spherical  cap  is  of  the  same  diameter, 
and  tangent  to,  the  cylindrical  leading  edges. 

The  examination  of  combined  simple  shapes  illustrates  the  adjustments 
that  must  be  made  in  applj^ing  basic  methods  to  complex  configurations.  It 
will  be  seen  that  in  some  cases,  as  for  example  the  slender  blunted  cone,  She 
combined  body  has  considerably  different  aerothermodynamic  characteristics 
than  would  be  Indicated  by  an  examination  of  each  element  separately. 

A.  Hemisphere-Cone 

1.  Local  Flow  Properties 

Nose  bluntness  affects  cone  heating  rates  by  altering  the  surface  pressure 
distribution  and  by  increasing  the  entropy  of  the  flow  at  the  boundary  layer 
edge.  Typical  Invlscid-flow  blunt-cone  pressure  distributions  are  presented 
in  figure  3“1*  The  values  shown  were  calculated  for  a  hemiphere-cone  using 
the  blunt  body  Inverse  method  (reference  ll)  for  the  subsonic  region  and  the 
method  of  characteristics  for  the  remainder  of  the  flow  field.  All  calculations 
are  for  real-gas  air  in  chemical  equilibritan.  Also  shown  ere  the  sharp-cone 
pressure  coefficients  calculated  using  figure  2-^. 


Aa  shown  in  figure  3“1^  bhe  presence  of  the  hemisphere  affects  the  pressure 
distribution  for  considerable  distances  downstream,  particularly  at  small  cone 
angles.  Theoretically,  one  expects  that  as  the  cone  angle  approaches  zero  the 
effect  of  nose  bluntness  will  persist  infinitely  far  downstream.  At  large  cone 
angles,  however,  the  effect  is  seen  to  be  confined  to  a  distance  of  a  few  nose 
diameters . 

Figure  3“1  is  apnliceble  only  for  cones  with  hemispherical  noses,  but 
qualitatively  siral-lar  behavior  is  to  be  expected  for  other  nose  shapes.  The 
blast  wave  analogy  proposed  by  Lin  in  reference  l4  and  applied  by  many  others 
(see  for  example,,  the  discussion  in  reference  5)  indicates  that  for  given 
free-stream  conditions  tne  distribution  on  a  blunt  cylinder  aligned  with  the 
flow  depends  onJ^  on  the  nose  drag  coefficient  and  not  on  the  details  of  the 
nose  shape.  This  conclusion  cannot  be  applied  in  the  vicinity  of  the  nose 
itself,  however,  since  the  analogy  breaks  down  as  the  shock  angle  approaches 
90  degrees.  Since  the  bluntness  effect  extends  such  a  short  distance  at 
large  cone  angles,  the  effect  of  nose  shape  on  the  details  of  the  pressure 
distribution  cannot  be  ignored  in  that  case.  However,  it  is  seen  that  the 
overall  effect  is  relatively  small  when  the  cone  angle  is  large. 
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In  contrast  to  the  effect  on  pressure,  the  effect  of  blxmtness  on  the 
entropy  of  the  inviscid  flow  theoretically  extends  infinitely  far  downstream 
for  all  cone  angles.  The  streamline  at  the  cone  surface  reaches  the  stagna¬ 
tion  point  after  having  passed  throu^  the  normal  shock  that  exists  ahead  of 
all  blunt  bodies  in  supersonic  flow.  In  crossing  the  shock  much  of  the 
kinetic  energy  of  the  stream  is  transformed  into  internal  energy,  with  the 
result  that  the  surface  streamlines  on  the  blunt  cone  have  a  lower  velocity 
and  density  than  at  corresponding  locations  on  sharp  cones,  even  though  the 
static  pressure  may  be  nearly  identical.  These  reductions  in  velocity  and 
density  lead  to  a  reduction  in  the  heating  rate  as  well.  The  changes  in 
local  flow  properties  and  the  reduction  in  heating  are  commonly  referred  to 
as  "bow  sliock  effects." 

In  calculating  the  bow  shock  effect  the  local  flow  properties  are 
obtained  by  isentropic  expemslon  fron  the  stagnation  point  to  the  local 
static  pressure.  Since  this  calculation  requires  knowledge  of  only  the 
total  enthalpy,  stagnation  point  pressure,  and  the  local  pressure.  It  is 
seen  that  the  magnitude  of  the  bow  shock  effect  does  not  depend  on  the  nose 
shape  or  size,  provided  only  that  the  nose  is  sufficiently  blunt  that  there 
exists  a  detached  shock  wave. 

The  existence  of  a  boundary  layer  on  the  cone  modifies  and  reduces  the 
bluntness  effect.*  A  streamline  actually  entering  the  boundary  layer  may 
originate  some  distance  away  from  the  normal  portion  of  the  bow  shock,  as 
shown  in  the  simplified  sketch  below: 


*  One  of  the  implications  of  this  statement  is  that  the  inviscid  flow  on  a 
bl\int  cone  does  not  approach  as  a  limit  the  conditions  of  a  local  flow  on  a 
sharp  cone  as  the  biuntness  Is  decreased  to  zero,  althou^  the  discrepancy 
is  confined  to  a  vanishingly  thin  layer  near  the  cone  surface. 
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since  the  mass  flov  within  the  boundary  layer  la  given  by 

Pe  “e 

It  is  seen  that  y^^  can  be  obtained  without  actually  tracing  the  streamline. 
Applying  a  mass  balance  leads  to: 

p  u  ffy^=7fr(R  +  6-  6*)^  -  R^l  p  u  =  2ir  R  (6  -  6*)  p  u  (3.1) 

00  00  00  L'  '  ^  06  '  ^  6*9 


or 


(6  -  6^) 


♦  Pe  ^e 


P  u 

CD  CD 


(3.2) 


If  the  variation  of  the  shock  wave  angle  with  3^  is  known,  the  entropy 
Just  downstream  of  the  shock  can  be  obtained  by  oblique  shock  theory.  Since 
the  static  pressure  at  the  edge  of  the  boundary  layer  is  equal  to  that  at  the 
cone  surface,  all  other  flow  properties  can  be  obtained  by  isentropic  expan¬ 
sion. 


The  actual  variation  in  heating  rate  along  the  surface  of  a  blunt  cone,  in¬ 
cluding  the  effect  of  the  boundary  layer  mass  absorption,  follows  the  blunt  cone 
theory  initially  and  then  Increases  toward  the  sharp  cone  theory.  Some  sli^t 
overshoot  probably  occurs  as  sketched  below: 


Distance 


Tlie  boundary  layer  thicknesses  8  and  8*  will  be  affected  by  the 
streanrwise  variation  in  edge  flow  properties,  restilting  in  a  difficult  } 

boundary  layer  calculation.  Very  little  has  appeared  in  the  literature  I 

regarding  the  effect  of  such  variation.  Since  the  boundary  layer  is  ustially  I 

assumed  to  be  "thin,"  it  is  to  be  expected  that  the  variation  in  local  flow  | 

properties  will  be  gradual,  and  adequately  accounted  for  by  using  the  (varying)  | 

local  properties  in  constant  property  boimdary  layer  formulas.  Such  a  procedure  | 

would  account  for  the  effect  of  changing  levels,  but  not  for  the  effects  of  | 

flow  property  gradients.  P'or  any  but  the  simplest  shapes  such  a  calculation 
is  very  difficult  and  the  results  far  from  certain. 

In  the  present  report,  blunt  cone  heating  rates  are  based  on  local  flow  * 

properties  calculated  by  Isentropic  expansion  frcan  behind  a  normal  shock, 
and  so  represent  the  lower  curve  of  the  above  sketch.  Further  information, 
including  some  experimental  results  for  flat  plate  turbulent  boundary  layers, 
may  be  obtained  from  reference  15 . 

2.  Streamline  Divergence  on  a  Blunt  Cone 

Bltintness  does  not  change  the  fact  that  transverse  pressure  gradients 
are  zero  if  the  cone  is  uuyawed.  As  a  result,  the  value  of  t/t^  =  1  and  by 
choosing  the  tangency  point,  or  shovilder,  of  the  hsmisphere-cone  combination 
as  the  location  of  the  reference  value  equation  (2.15)  can  be  written: 


sin  Q _ 

Bln  egg 

“  1  +  -^  tan  6 


X  <  0 


X  ^  0 


(3.3a) 

(3.3b) 


On  the  blimt  delta  wing,  as  on  the  blunt  cone,  the  hemispherical  nose 
affects  downstream  pressures  and  through  the  entropy  effect,  all  other  flow 
properties.  The  evaluation  of  these  effects  is  much  more  complex  than\for 
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the  cone,  varying  with  spanwise  position  as  well  as  chordwise  distance.  In 
addition,  hluntness  affects  the  delta  wing  streamline  patterns,  and  the  wing 
in  tvirn  has  ein  effect  on  the  nose  and  leading  edge  flow.  The  calculation  of 
tiie  flow  field  surrounding  a  blunt  delta  wing  is  more  difficult  than  for  the 
sharp  delta  wing,  since  the  previously  mentioned  condition  of  constant  flow  pro¬ 
perties  along  ray  lines  can  no  longer  be  applied.  In  principle,  the  flow  can 
be  calculated  by  the  method  of  characteristics  for  three  dimensions,  using 
blunt  body  theory  to  obtain  the  initial  values.  In  practice,  however,  it  was 
found  that  a  large  amount  of  computer  time  was  required,  and  the  computer 
results  obtained  showed  indications  of  numerical  Instabilities. 

This  lack  of  exact  theoretical  methods  for  the  blunt  delta  wing  has  re¬ 
quired  the  use  of  approximate  methods  and  empirical  results.  The  bases  for 
the  approximate  methods  are  the  basic  shape  methods  of  the  previous  section. 

The  empirical  information  is  based  on  the  results  of  an  extensive  series  of 
delta  wing  tests  that  were  conducted  during  the  X-20  program  covering  a 
Mach  number  range  from  6  to  22.  The  data  from  those  tests  has  been  reanalyzed 
under  NASA  contract  and  will  be  published  as  references  l6,  17  and  l8. 

The  effect  of  bluntness  on  centerline  pressure  distribution  is  shown  in 
figure  3-2.  These  curves  were  obtained  by  matching  empirical  results 
from  the  X-20  tests  to  the  hemisphere  method  previoiisly  given.  Also  shown 
in  figure  3-2  are  sharp  delta  wing  values  obtained  from  figure  2  -U.  The 
blunt  delta  wing  pressures  are  appreciably  higher  than  those  for  the  sharp 
delta  wing  for  at  least  three  nose  diameters,  and  for  angles  of  attack  greater 
than  approximately  30°)  do  not  fair  into  sharp  delta  wing  values.  The  solid 
curves  of  figure  3-2  have  been  used  for  the  heat-transfer  curves  of  Section 
VI. 


Near  the  centerline  the  effect  of  nose  bluntness  on  entropy  can  be 
calculated  as  described  in  the  discussion  of  hemisphere-cones.  However, 
calculations  show  that  the  effect  of  the  highly  swept  leading  edges  on  down¬ 
stream  entropy  levels  is  much  less  than  for  the  hemisphere,  and  has  a  much 
smaller  effect  on  the  heating  rates.  The  proper  evaluation  of  bluntness 
effects  on  a  delta  wing  therefoi'e  requires  a  detailed  know’ledge  of  the  be¬ 
ginning  of  each  streamline  that  wets  the  surface,  a  knowledge  that  at  present 
does  not  exist.  Since  the  bluntness  effect  reduces  downstream  heating  rates, 
it  is  conservative  to  neglect  the  bluntness  effect  except  at  the  centerline, 
where  the  surface  streamlines  are  known  to  have  come  through  the  hemisphere 
shock.  The  curves  presented  in  this  report  for  blunt  delta  wing  o>enterlirie 
heating  are  based  on  normal  shock  entrojir,  as  were  the  blunt  cone  curves. 

The  method  used  for  calcalatlng  spanwise  distributions  for  blunt  delta  wings 
will  be  described  following  a  discussion  of  streamline  methods. 


25 


Bluntness  effects  on  the  streamline  patterns  were  experimentally  Investi¬ 
gated  by  aoi  oil-flow  technique.  A  mixture  of  oil  and  lampblack  was  spread 
uniformly  over  the  surface  of  test  models  before  they  were  subjected  to  the 
hypersonic  tunnel  flow.  Patterns  were  seen  to  develop  that  Indicate  the 
direction  of  the  surface  shear  force.  These  patterns  were  then  either  photo¬ 
graphed  or  transferred  to  velltun  paper  by  direct  contact.  An  analys^.s  was 
made  of  the  relation  of  these  observed  shear  force  patterns  to  the  streamline 
directions  in  the  external  flow  which  led  to  the  conclusion  that  the  surface 
shear  patterns  always  indicate  more  outflow  than  actually  exists  in  the 
external  flow.  The  argument  leading  to  this  conclusion  may  be  STsnmarlzed 
as  follows: 

1.  A  fluid  particle  in  the  external  flow  describes  a  cvtrved  path  such 
that  the  centrifugal  force  acting  on  the  particle  is  in  equilibrium 
with  the  transverse  pressure  gradient. 

2.  The  same  transverse  press\ire  gradient  acts  on  particles  within  the 
boundaiy  layer.,  since  the  pressxire  is  constant  throu^  the  boundary 
layer . 

3.  Since  the  velocity  of  the  particles  within  the  boundary  layer  is  less 
than  those  in  the  external  flow,  the  centrifugal  force  will  not  match 
the  external  pressure  gradient  unless  the  path  of  the  part, ides  within 
the  boundary  layer  has  a  smaller  radius  of  curvature  than  that  in  the 
external  flow. 

Hence,  the  surface  oil  flow  patterns  (which  reflect  the  surface  shear 
forces)  will  indicate  more  curvatvire  than  exists  in  the  extenial  flow.  This 
conclusion  would  apply  in  either  laminar  or  turbulent  flow.  The  oil  flow 
results  showed  that  blrxntness  greatly  reduced  the  amount  of  outflow,  with  the 
streamlines  being  nearly  peurallel  to  the  centerline  in  proximity  to  the  nose 
even  at  high  angles  of  attack.  Since  outflow  increases  the  heating  rates  by 
thinning  the  boundary  layer,  it  is  conservative  to  use  sharp  delta  wing  stream¬ 
line  methods  on  the  blunt  delta  wing.  This  has  been  done  in  preparing  this 
report.  Since  the  heating  rate  calculation  Involves  strearawise  integration 
of  the  boundary-layer  momentum  equation  beginning  at  the  stagnation  point  (as 
described  in  Appendix  B)  delta  wing  streamline  assumptions  do  not  affect  the 
prediction  at  the  hemisphere  shoxuder.  It  is  also  to  be  expected  that  far 
downstream  where  the  wing  si)an  is  much  greater  than  the  leading  edge  radius 
the  s)iarp  delta  wing  streamline  method  should  be  applied.  Thus,  the  use 
of  s^arp  delta  wing  streamline  methods  in  the  region  just  downstream  of  the 
nose  results  in  a  prediction  that  begins  at  sharp  delta  wing  values.  Far 
downstream  the  prediction  approaclxes  the  same  asymptotic  values  as  would  be 
calculated  using  a  method  based  on  estimated  blunt  wing  streamlines,  but 
the  prediction  would  be  consei^ative  in  the  intermediate  region. 

Using  the  sharp  delta  wing  streamline  prediction  method,  the  streamline 
parameter  A  in  the  vicinity  of  the  centerline  is  given  by 

=  (1  +  I  cos  a)  (3.4) 
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where  R  is  the  nose  radius,  x  is  measured  along  the  centerline,  he  ing  at 
the  tangency  point  between  the  nose  cap  end  the  plane  surface  and  given 

by  equation  (2.28). 

2.  Blunt  Delta  Wing,  Leading  Edge  and  Nose  Heating 

The  blunt  delta  wing  discussion  has  thus  far  considered  nose  and  leading 
edge  effects  on  the  del-oa-wing  flow  field.  As  long  as  the  wing  flow  field  is 
entirely  supersonic,  both  normal  and  tangential  to  the  leading  edge  in  the 
plane  of  the  wing,  it  is  to  be  expected  that  the  flow  over  the  nose  and 
leading  edges  will  be  unaffected  by  the  lower  surface.  The  sonic  points  on 
both  the  leading  edge  and  nose  are  about  U5  degrees  from  the  stagnation  point, 
so  that  the  lower  surface  should  begin  to  affect  the  leading-edge  flow  field 
when  the  effective  angle  of  attack  given  by  equation  (2.16)  reaches  U5 
degrees.  Such  an  effect  on  the  location  of  the  leading-edge  stagnation  line 
was  observed  in  the  delta -wing  tests  as  illustrated  in  the  sketch  below: 


A=  73“ 


Angle  of  attack,  a,  degrees 


The  upper  curve  is  a  faired  curve  representing  oil  flow  data  measured 
on  an  un/awed  73-degree  swept  blunt  delxa  wing  at  Mach  numbers  of  6.0  and  8.08. 
The  lower  curve  can  be  computed  using  equation  (2.17)^  since  is  numeric¬ 

ally  equal  to  6_.  The  curves  show  that  at  angjies  of  attack  less  than  about 
15°  the  lower^suvface  has  no  influence  on  leading-edge  stagnation-line 
location.  However,  as  the  angle  of  attack  is  increased,  lower-surface  pres¬ 
sure  influences  the  location  of  the  stagnation  line  indicating  that  it  will 
reach  the  shoulder  at  a  lower  ang3.e  of  attack  than  predicted  by  the  infinite 
cylinder  theory. 
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During  the  X-20  program  an  approximate  method  was  developed  that  pre¬ 
dicts  the  angle  of  attack  for  which  the  stagnation  line  is  located  at  the 
shoulder.  The  method  is  based  on  the  observation  that  when  the  stagnation 
line  is  on  the  leading  edge  shoulder  the  velocity  in  the  plane  normal  to  the 
leading  edge  axis  is  zero  at  the  leading  edge  shoulder.  In  this  crossflow 
plane  the  velocity  at  the  shoulder  was  estimated  by  superimposing  the  previously 
described  delta  wing  outflow  velocity  onto  the  leading  edge  velocity  distri¬ 
bution.  Closest  agreement  with  data  obtained  during  the  X-20  experimental  program 
occurred  when  the  circumferential  velocity  distribution  on  the  cylinder  was 


calculated  by 


(3.5) 


where  6  is  measured  from  the  zero  angle  of  attack  stagnation  line, 
velocity  gradient  was  calculated  from 


(  ^\  -  al  \  fix 

R  /  de 


(3.6) 


a  -  R  .  SL  '  “  ' 

where  P  “  u  ^  (which  depends  only  on  the  normal  Mach  number)  was  calculated 
for  the  true  nofroal  Mach  number  as  given  by  equation  (2.17a)*  The  nomal 


velocity 


is  calculated  as  if  for  zero  angle  of  attack  by 


-  u  cos  A 
00 


(3.7) 


Setting  the  shoulder  value  of  v  equal  to  the  deltt-wing  outflow  velocity  as 
given  by  equation  (2.25)  leads  to: 


tan  a  "  I  A 


(3.8) 


Equation  (3.8)  gives  that  angle  of  attack  at  which  the  leading-edge 
stagnation  line  is  located  at  the  leading-edge  shoulder.  Stagnation  line 
locations  in  this  report  have  been  approximated  by  using  equation  (2.17)  for 
angles  of  attack  up  to  15*  and  then  fairing  a  line  to  the  value  obtained  using 
equation  (3.8). 


A  second  lower-surface  effect  on  the  leading  edge  that  was  first  noted 
early  in  the  X-20  program  was  a  reduction  in  stagnation-line  heating  rates. 

Such  a  reduction  is  to  be  expected  when  it  is  noted  that  at  very  high  angles 
of  attack  equation  (2.24)  predicts  that  the  flow  near  the  leading  edge  is  out¬ 
ward.  In  that  case  the  leading  edge  behaves  as  a  trailing  edge,  and  a  reduc¬ 
tion  in  heating  is  to  be  expected.  No  reliable  methods  were  developed  for  pre¬ 
dicting  the  magnitude  of  the  reduction,  however,  since  there  was  little  affect 
on  the  X-20  operating  envelope. 

At  about  45  degrees  angle  of  attack  the  delta-wing  lower  surface  influences 
the  subsonic  portion  of  the  nose  hemisphere  flow  field,  and  a  similar  reduction 
in  stagnation  point  heating  is  to  be  expected.  However,  the  effect  is  expected 
to  be  small,  and  no  evaluation  of  this  effect  was  made  during  the  X-20  program. 
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3.  Blunt  Delta  Wing  Spanwise  Heating  Rate  Distribution  -  Leading  Edge 

Stagnation  line  locations  shown  in  the  charts  of  Section  VI  were  deter¬ 
mined  using  the  method  Just  described.  Stagnation  line  heating  rates  were 
determined  by  infinite  cylinder  theory  for  flight  conditions  at  angles  of 
attack  less  than  20  degrees.  The  aforementioned  reduction  la  heating  rates 
at  the  stagnation  line  in  canparlson  to  infinite  cylinder  theory  is  based  on 
a  73 "degree  swept  delta  wing  data  and  is  purely  empirical. 

The  blunt  delta  wing  heating  curves  shown  in  Section  VI  from  the 
stagnation  line  around  the  leading  edge  and  onto  the  surface  were  defined  by 
t,he  data  of  reference  l6.  The  actual  magnitude  of  the  heating  rate  ratios 
presented,  however,  were  shifted  in  order  to  fair  smoothly  into  theoretical 
centerline  values. 

C .  Deflected  Flap 

1.  Compression 

Flaps  deflected  into  i  stream  compress  the  oncoming  flow  and  often  cause 
the  boundary  layer  ahead  of  the  flap  to  separate.  Effects  of  flow  separation 
and  reattachment  on  aerodynamic  heating  can  be  large  and  extremely  complex. 

In  most  cases  predictions  must  be  based  on  approximate  methods.  The  difficulty 
of  providing  useful  heat-transfer  information  for  actual  design  is  made  some¬ 
what  easier  by  the  fact  that  maximm  values  are  usually  of  most  interest. 

Aerodynamic  heat  transfer  through  a  laminar  separaced  region  was  analyzed 
in  reference  19  by  an  extension  of  boundary  layer  theory.  That  calculation 
indicated  that  the  average  heat  transfer  in  a  separated  region  is  reduced  by 
about  50  percent  as  compared  to  attached  flow  heat  transfer  for  the  same 
local  flow  properties  at  the  boundary  layer  outer  edge.  This  result  has 
been  verified  in  tests  of  reference  20  for  cavity-type  flow  and  tests  of 
reference  21  for  forward -facing  step  flow.  The  results  of  reference  20 
indicate  that  the  same  reduction  might  apply  to  turbulent  flow.  Test  data 
of  reference  21,  however,  indicate  an  increase  in  heating  for  separated 
transitional  and  turbulent  flow  over  forward-facing  steps. 

The  prediction  of  reattachment  heating  is  difficult  because  of  large 
pressure  gradients  that  exist  throughout  the  reattachment  region.  Calcula¬ 
tions  (reference  22)  have  been  made  for  laminar  boundary  layer  flow  at 
reattachment;  however,  a  prior  knowledge  of  the  reattachment  pressure 
gradient  is  required.  Since  regions  of  separated  flow  and  locations  of 
reattachment  are  practically  Impossible  to  predict  ior  complc’x  configurations, 
an  approximate  method  was  developed  during  the  X-20  program  to  predict  the 
upper  bound  of  reattachment  heating.  This  method  relates  the  maximum  inviscid 
pressure  rise  on  the  deflected  flap  to  a  maximum  heating  rise  at  the  point  of 
reattachment  as  follows: 

=  ^MAX  (3-9) 

^JWDISTURBED  ^UMDISTURBED 
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This  approximate  result  is  based  on  the  relation: 

X  i  P*  u 
,  constant  f  e 
h  = 


Pr 


r 


r  p*  M*  Ue 

jo 


l/(n+l) 


(3.10) 


which  is  a  slight  generalization  of  an  equation  shown  in  reference  23.  In 
equation  (3*10)  n  is  1  for  laminar  flow  and  4  for  turbulent  flow.  The  super¬ 
script  *  denotes  evaluation  at  the  reference  temperature,  defined  by  Eckert 
(reference  as; 


T*’  =  0.5  +  0.28  Tg  +  0.22 


To  evaluate  the  effect  of  a  sudden  compression  on  h  as  predicted  by  equation 
(3.10),  we  write 


h^  (p*M*l/«Ug)^ 


(p* 


r 


p*  p*  xXq  6x 


f  p*  fi*  Uq  dx 


1/n+l 


(3.11) 


where  the  subscripts  +  and  -  indicate  evaluation  Just  downstream  and  Just  up¬ 
stream,  respectively,  of  a  sudden  compression.  If  the  compression  occurs  over 
a  very  short  distance,  the  two  integrals  must  be  nearly  equal,  since 
x+  x_ 

J  p*  p*  vlq  dx  =  J  p*  p*  Uq  dx  +  (x+  -  X.)  p*  p*  Uq  . 

o  o 

and  for  small  values  of  (x^  -  x_,) 


so  that 


J  p*  p*  Ue  dx  =  J  p*p*nedx 

o  o 

!!t  .  (P*  ^  [Pe(M*^/»/T*)  uj^ 

h_  “  (p*  p*l/n^^)_  "  (p*1/YT*)  Ue]_ 


(3.12) 

(3.13) 


In  this  expression  the  change,  in  pressure  are  dominant.  The  changes 
in  the  reference  temperature  T  *  are  small  because  T^  and  T^^^  do  not  change 

appreciably.  In  addition,  the  changes  in  T*  and  p  *  tend  to  compensate. 
Refex'ring  to  equation (2 .20),  it  is  found  that  the  ratio  of  local  to  free- 
stream  velocity  is  approximately  vmity  for  angles  of  30  degrees  or  less. 
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Neglecting  the  "small"  differences  in  velocity  and  temperature  leads  to: 

^  ^  (3.,„ 

[p^(^*l/n/Tt)  uj_  (P). 

Immediately  downstream  of  the  compression  the  heating  rate  begins  to  decrease 
as  the  integral  in  equation  (3 <10)  increases.  Hence, 

^MAX  ^MAX  ^  ^ 

h_  ^UNDISTURBED  P_  ^UNDISTURBED 

Since  no  assumption  has  been  made  regarding  boundary  layer  state,  equation 
(3.9)  applies  in  either  laminar  or  turbulent  flow.  Equation  (3.9)  would  also 
be  applicable  in  the  presence  of  flow  separation  provided  that  no  appreciable 
increase  in  the  integral  of  equation  (3-11)  occurred  over  the  separated  region. 
Since  the  integral  represents  the  effect  of  wall  shear  on  the  boundary  layer 
growth,  it  seems  reasonable  that  the  small  shear  forces  in  the  separated 
region  are  also  negligible.  In  the  case  of  separation  the  subscripts  +  and 
would  refer  to  conditions  Just  ahead  of  separation  and  Just  downstream  of 
reattac)-^ent .  Although  the  development  of  equation  (3.9)  involves  several 
gross  approximations,  it  has  been  \i8ed  in  the  absence  of  any  better  method. 
Ccs!!parls.on6  of  equation  (3.9)  with  experiiriental  data  are  presented  in 
reference  1?. 

2 .  Expansion 

Aerodynamic  heating  rates  to  flaps  in  an  expanding  flow  field,  while 
less  than  those  obtained  on  compression  surfaces,  may  be  required  for  design. 
The  presence  of  the  boundary  layer  complicates  the  flow  field  by  creating  a 
displacement  effect  that  changes  the  effective  expansion  angle.  To  deter¬ 
mine  this  angle  in  a  two-dimensional  flow  field,  boundary-layer  growth  must  be 
taken  into  account. 


Consider  the  flow  illustrated  in  the  following  sketch: 
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In  the  absence  of  a  boundary  layer  the  effective  turning  angle  is 

and  the  flow  properties  at  the  wall  can  be  determined  using  the  Prandtl-Meyer 
equation 


(3.15) 


where  H  is  the  difference  between  total  and  local  enthalpy  and  a  is  the  speed 
of  sound  evaluated  using  local  enthalpy. 

In  the  presence  of  a  boundary  layer,  equation  (3 •15)  is  valid  only  for 
that  portion  of  the  turning  angle  external  to  the  boundary  layer.  The 
effective  expansion  angle  is  given  approximately  by 


c  d6* 

"  *off  "  dx 


(3.16) 


Equation  (3«l6)  was  used  In  preparing  this  report.  However,  this 
equation  is  exact  only  if  the  boundary-layer  edge  flow  properties  are 
constant.  No  exact  simple  formulation  exists  if  the  external  flow  properties 
are  not  constant  due  to  the  existence  of  nonisentropic  flow  within  the 
boundary  layer. 

Consistent  values  of  Ai>  and  8*  were  obtained  by  an  interative  process 
in  which  equation  (3.15)  was  first  used  to  calculate  pressure  neglecting  the 
boundary  layer  effect  on  turning  angle.  The  displacement  thickness  was  then 
determined  using  the  Nonsimilar  Boundary  Layer  Program.  Equation  (3 *10) 
was  then  used  to  determine  a  new  turning  angle.  This  process  was  repeated 
until  further  changes  were  negligible. 

Since  the  iteration  process  is  lengthy  the  following  approximate  closed 
form  I'elationship  was  fitted  to  the  numerical  results; 


(3.17) 


where  h^^^  is  the  heat  transfer  coefficient  on  the  expanded  surface,  h^j^^ 

is  the  heat  transfer  coefficient  at  the  hinge  line,  1®  inviscid 

pressure  ratio  across  the  expansion  wave,  is  the  ratio  of  the  distance 

to  the  hinge  1^  .e  to  the  distance  from  the  leading  edge  of  the  flat  plate. 

For  laminar  flow  n  is  0.5  and  for  turbulent  flow  n  is  0.2.  The  exponent  m 
has  a  value  of  approximately  0.8  for  laminar  flow  and  approximately  O.9  for 
turbulent  flow,  However,  due  to  the  fact  that  equation  (3 -I?)  involves  several 
approxirations,  a  was  conser/atively  assumed  equal  to  0.8  for  both  laminar  and 
turbulent  flow. 


Use  of  equation  (3.1?)  is  limited  to  two-dimensional  shan)-flat-plate-flap 
combinations  at  angles  of  attack  greater  than  zero  degrees. 
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DISTANCE  FROM  STAGNATION  POINT,  S/R 


Figure  3-2:  BLUNT  DELTA  WING  LOWER  SURFACE  CENTERLINE  PRESSURE 


IV.  SURFACE  CONDITION  EFFEJCTS 

In  general  it  is  not  to  toe  expected  that  the  surface  panels  of  a  re-entry 
vehicle  will  have  the  Ideally  smooth  surfaces  assumed  in  the  toasic  heat  trans¬ 
fer  analysis.  Roughnesses  resulting  from  manufacturing  tolerances  which  may 
^llow  protruding  fasteners,  panel  mismatch,  and  panel  surface  curvatxxre  can 
cause  flow  field  disturbances.  The  severe  thermal  environment  may  lead  to 
buckling  of  surface  panels]  further,  shallow  corrugations  or  surface  waves  may 
toe  required  to  control  the  orientation  of  skin  buckling  that  results  from 
thermal  expansion.  It  is  also  unlikely  that  perfect  sealing  can  toe  provided  at 
panel  Joints  and  particularly  at  the  hingelines  of  movable  control  surfaces. 
Leakage  at  such  points  can  bleed  off  the  relatively  cool  air  in  the  lower  part 
of  the  boundary  layer  and  so  cause  increased  heating  rates. 

In  the  X-20  program  the  Nonslmllar  Boundary  Layer  Program,  Appendix  C, 
was  used  to  calculate  the  effects  of  leakage  and  shallow  surface  waves  on 
iPTniTiar  heating  rates.  The  results  of  those  calculations  are  described  in 
this  section.  For  all  other  roughness  effects  it  was  necessary  to  rely  on 
experimental  results.  The  surface  rouj^ness  tests  were  not  completed  at  the 
termination  of  the  X-20  program,  nor  was  the  analysis  completed  of  those  tests 
that  had  been  made.  The  results  available  at  the  time  led  to  the  following 
conclusions : 

1.  Forward  facing  steps  are  not  acceptable  for  desiga. 

2.  An  allowance  of  20  percent  for  the  effect  of  aft  facing  steps  was 
required . 

3*  W'nere  notches  were  unavoidable,  they  should  be  combined  with  aft 
facing  steps. 

A  detailed  test  of  aft  facing  step  and  notch  configurations  was  in  progress 
at  the  time  of  the  termination  of  the  X-20  contract.  The  results  of  the 
s'orface  roughness  tests  are  being  reaneJ-yzed  by  The  Boeing  Company  under 
contract  to  NASA,  and  will  be  published  in  reference  l8. 

A.  Shallow  Surface  Waves 

Calculations  of  viscous  interaction  effects  on  the  pressure  and  boundary 
layer  over  a  shallow  wave  were  made  by  ccmibining  the  results  of  linearized 
theory  and  the  Nonsimilar  Boundary  Layer  Program.  An  Inverse  method  was  em¬ 
ployed  using  linearized  theory  for  the  pressure  distribution,  and  the  Non¬ 
similar  Boundary  Layer  Program  to  csdculate  laminar  boundary  layer  profiles. 

A  method  essentially  the  same  as  the  oiie  used  in  this  report  was  developed  by 
Baxter  and  Flugge-Lotz,  reference  25* 
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A  distribution  of  presenre  was  assuned  and  Ay  found  fron  linear  theory  with  the 
aid  of  equation  (4.l): 


Cp  = 


i/Me^  -  1 


(4.1) 


Usings  the  Nonsimilar  Boundary  Layer  Program  the  effect  of  the  assumed  pressure 
distribution  on  the  boundary  layer  can  be  accurately  calctilated.  Then  the 
height  of  the  surface  wave  corresponding  to  the  assumed  pressure  distribution 
is  given  by: 


«  =  yref  +^y  -  6*1 


^e^FP 
3  ^  e 


(yref  ■  ^*FP) 


(4,2) 


When  the  resxxlts  of  many  such  calculations  were  examined  it  was  found  that  the 
followj.ng  equations  could  be  used  to  predict  the  maximum  laminar  heating  rates: 


^MAX 


j  |-  Ki  +  [ki2  f  4K2]^^^| 


(4.!,) 


where 


Ki  =  [i  +  ABC  +  (U/6pp)  (BCD  -  1)]  j -t  R/6^jjBC)] 


K2  =  (R/fipp)  B^Cac  +  R/6pp(CD)] 


(4.4) 


k,  B;  C,  and  D  are  defined  as 
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Th«  folloving  numerical  evaluations  of  the  parametera  A,  B,  and  C  were /made 
through  correlation  of  Nonaimilar  Boundary  Layer  Program  results. 


A  « 


_ 0.3 

0.78  +  0.84 


B  =  2.5  (0.78  +  0.84  I^/Iq) 


(0.78  +  0.84  i^/io) 


Cujrves  are  given  in  Section  VI  that  provide  numerical  results  in  laminar 
flow  based  on  equations  (**^.3)/  (^•^)  and  (4,5). 


No  analytical  method  had  been  developed  to  predict  the  effect  of 
surface  waves  in  tvirbulent  flow. 


B.  Leakage 

The  effect  of  leakage  on  laminar  heating  rates  was  also  calculated  with 
the  Nonsimilar  Boundary  Layer  Program.  Calculations  were  made  for  constant, 
parabolic,  and  sinusoidal  leakage  velocity  distributions. 
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I'he  leokai^e  velocities  assuiuea  are  illustrated  by  the  following  sketch: 


where  ra  is  the  rate  at  which  mass  leaves  the  boundary  layer,  is  the 

leakage  velocity,  andu^is  the  boundary  layer  edge  velocity.  Downstream  laminar 

heating  rates  were  computed  using  the  Nonslmllar  Boiindary  Layer  Program, 

The  results  were  well  correlated  by  the  following  equations  independently  of 

the  distribution  of  : 

m 


_ ^  ^  0.494  L 

^NL 


h  -  rii  (Pw  Mw  «e  «l) 


-1/3 


*NL 


where  H,  =  - — = - — 

^  ^AW  - 


(4.7) 

(4.8) 


The  definitions  of  the  above  terms  are  shown  in  figure  4-1.  Equation 
(4.6)  was  used  to  prepare  the  leakage  effect  curve  of  Section  VI. 


Fo"  turbulent  flow  no  analytical  method  existed  for  nredicting  the 
effect  of  leakage  on  heating  rates. 


I 
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Figure  4-1:  EFECT  OF  LEAKAGE  ON  LAMINAR  HEAT  TRANSFER 


V.  .APPLICATION  TO  COMPLEX  CONFIGURATIONS 

The  aerotheraodynamlc  environment  of  a  re-entry  vehicle  is  highly 
configuration  and  trajectory  dependent,  and  no  set  of  generalized  rules 
can  be  expected  to  insure  an  accurate  and  thorou£^  emalysis.  Rather, 
each  vehicle  must  be  examined  indlvi.dually  and  the  adequacy  of  the  analysis 
will  depend  in  good  part  on  the  experience  and  care  of  the  engineer  making 
the  «tudy.  However,  when  properly  applied  the  charts  of  the  following 
section  can  provide  good  Ini^'lal  estimates  of  heating  rates  for  typical 
hypersonic  vehicle  configurations  csing  considered  today.  As  a  minimuin 
the  following  steps  should  be  taken: 

1.  Make  preliminary  estimates  by  dividing  the  vehicle  Into  geometric 

elements  corresponding  as  closely  as  possible  to  the  shapes  discussed 
in  Sections  II  and  III,  euid  apply  the  coiresponding  curves  of  Section  VI. 

a.  Fjcamine  the  relation  of  each  element  to  the  others.  If  an  element  is 
within  the  flow  field  of  another  a  detailed  calculation  of  the  local 
flciw  should  be  made  to  determine  the  local  pressure  increase.  Then 
rauj.xlply  the  previous  estimate  of  the  local  heating  rate  by  the  presture 
ratio  30  obtained.  (This  calculation  is  particularly  applicable  to 
locsil  re-entrsint  comers,  such  as  deflected  trailing  edge  flaps.) 

3.  Estimate  the  location  of  shock  waves  to  deteimine  if  any  will  impinge 
on  any  other  part  of  the  configuration.  There  is  no  simple  method  for 
calculating  the  effects  of  such  iiapingeraents  in  all  cases,  but  in  the 
absence  of  more  definitive  information  it  should  be  assumed  that  a 
significant  heating  rate  increase  will  result. 

Examine  each  location  at  which  the  heat  protection  system  changes  in  any 
way,  and  look  for  any  design  feature  that  could  interfere  with  the 
action  of  the  heat  protection  system.  For  example,  on  a  radiantly  cooled 
vehicle,  examine  each  location  at  which  the  material  temperature  limit 
changes,  and  look  for  radiation  blockage  by  nearby  configuration  com¬ 
ponents  . 

5.  Estimate  the  effect  of  the  surface  condition  during  re-entry  on  smooth 
body  heating.  No  simple  theory  exists  for  any  roughness  element  dif¬ 
ferent  than  a  wave.  The  nile  followed  for  the  X-20  was  to  try  and 
eliminate  grooves  tnd  steps  where  possible  and  attempt  to  control  the 
orientation  of  surface  buckling. 

It  is  emphasized  that  the  above  are  only  the  preliminary  steps,  with  the 
course  of  further  analysis  dependent  on  the  results  of  these  initial  calcu¬ 
lations.  For  example,  step  2  may  lead  to  unnecessarily  conservative  values, 
but  win  highlight  potential  problem  areas.  Thus,  on  the  highly  swept  delta 
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ving  glider  shovm  in  figure  5“1>  step  2  in  the  above  procedure  would  direct 
attention  to  the  elevens  and  canopy,  while  step  3  would  indicate  that  the  fin 
may  be  critical  at  low  angles  of  attack.  On  the  hypersonic  airbreather  shown 
in  figure  5-2  the  critical  regions  indicated  would  be  the  wing  leading  edge 
at  its  Junct\ire  with  the  body  and  at  the  impingement  i>oint  of  the  body  shock 
wave,  the  entire  cowl,  the  tip  fins,  and  the  interior  of  the  engine  where  re¬ 
radiation  of  heat  is  blocked  by  the  enclosed  geometry.  Two  examples  from  the 
X-20  program  that  illustrate  the  need  of  detailed  analysis  will  now  be 
described . 


A.  Control  Surface  Gap 

If  aerodynamic  flaps  are  to  be  used  for  control  they  must  be  allowed  to 
move  freely.  It  was  found  in  the  X-20  studies  that  free  movement  could  be  pro'^ 
vided  only  by  preventing  contact  between  the  outboard  edge  of  the  eleven  and 
the  fin.  When  due  allowance  for  thermal  and  loeid  deformations  had  been  made 
it  was  found  that  a  sizeable  gap  existed  between  the  two  surfaces,  and  that 
flow  through  the  gap  would  occur  at  many  fli^t  conditions.  Although  the 
resulting  aerodynamic  heating  was  not  severe,  excessive  temperatures  occurred 
because  the  heated  surfaces,  being  in  close  proximity,  were  unable  to  re- 
radiate  effectively.  Thus  it  became  necessary  to  obtain  accurate  heating 
rate  information  for  the  inside  of  the  gap. 

Since  it  was  unlikely  that  ao  analyt.ical  investigation  alone  would  lead 
to  an  acceptable  fin-elevon  design,  a  series  of  tests  were  conducted  specifi¬ 
cally  to  solve  this  problem.  A  large  number  of  fin  lower  edge  shapes,  eleven 
outboard  edge  shapes,  and  gap  widths  were  tested  to  determl’"'  the  combined 
effects  of  aerodynamic  heating  and  improved  radiation  view  factor  on  gap  tem¬ 
peratures.  The  analysis  of  experimental  data  and  radiation  view  factors  in¬ 
dicated  that  the  point  of  maximum  surface  temperature  occurs  very  near  the 
point  of  minimum  fin-elevon  separation  for  fins  and  elevens  having  cylindrical 
surfaces,  and  near  the  Junction  of  the  entrance  emd  plane  sections  for 
parallel  fin-elevon  surfaces.  Gap  width  was  found  not  to  affect  pressure 
distributions  as  long  as  the  flow  was  not  choked  by  boundary  layer  growth. 

The  above  information  forms  a  basis  for  the  analysis  in  this  report 
which  considers  the  flow  fields  associated  with  the  three  gap  configurations 
presented  in  figure  5“3- 

Maximum  temperature  will  occur  on  these  gap  configurations  approximately 
90*  from  the  shoulders  Indicated  on  figure  5"3*  The  pressure  distribution 
used  in  calculating  a  heating  rate  distribution  by  the  p  method  is  shown 

in  figure  2-1.  Stagnation  conditions  were  based  on  the  normal  velocity  com¬ 
ponents, 

^n  =  ^oo«^”“  (5-1) 

where  is  the  stream  velocity  and  a  is  the  vehicle  angle  of  attack.  Tiie 
heating  rate  90*  from  the  shoulder  is  used  for  the  design  chart  shown  in 
Section  VI. 
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B. 


Fin  Leading  Edge  and  Canopy  Interference  Effects 


The  vehicles  illustrated  in  figures  5“1  and  5-2  have  vertical  fins  that  are 
to  be  used  for  aerodynamic  control  and  canopies  forming  the  cockpit  area. 
Experimental  data  obtained  during  the  design  of  the  X-20  indicated  that  at 
lov  angles  of  attack,  fin  leading-edge  and  canopy  heating  rates  would  b" 
considerably  higher  than  predicted  by  swept -infinite-cylinder  theory  or  local 
oblique-shock  theory  based  on  free-streajc  conditions.  ^Hie  higher  heating  rates 
are  attributed  to  the  shock  wave  generated  by  the  wing  surface  interfering 
with  the  leading  edge  and  canopy  surface  flow  field. 

The  exact  nature  of  the  Interference  effect  has  not  betn  established,  even 
at  the  present  time.  However,  based  on  the  results  of  the  calculations  de¬ 
scribed  below,  it  appears  that  the  primary  effect  is  the  change  in  local 
flow  properties  caused  by  the  wing  shock  envelope.  The  flow  field  is  iUus- 
ti>ated  below: 


The  effect  on  the  canopy  heating  rate  was  estimated  by  using  oblique  shock 
theory  to  calculate  local  flow  properties  Just  downstream  of  the  shock.  The 
canopy  heating  rate  was  then  calculated  with  infinite  cylinder  theory  using 
the  local  flow  properties.  It  was  found  that  for  each  angle  of  attack  the 
effect  of  the  oblique  shock  varied  such  that  there  was  a  most  unfavorable 
angi.e,  as  Illustrated  below: 

LUm  of  MiTllI 

,  bottlaf  r«t«a 
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Since  the  oow  shock  is  cui*ved,  it  was  assumed  that  at  some  point  the  most 
unfavorable  angle  would  occur,  and  that  therefore,  the  maximum  canopy  heating 
rate  would  be  given  by  the  maxima  of  the  above  family  of  curves.  Good  agree- 
■nent  with  test  data  was  obtained  by  this  method. 

It  will  be  noted  that  for  given  free- stream  conditions  the  final  result 
depends  only  on  the  true  sweep  angle,  X,  of  the  canopy,  (since  the  angle  of 
the  shock  is  assumed  to  take  on  all  possible  values)  and  thus  can  be  plotted 
as  a  single  curve.  The  curve  so  obtained  was  also  applied  to  fin  leading 
edges,  and  again  good  agreement  with  test  data  was  obtained.  The  agreement 
with  the  fin  data  may  be  fortuitous,  however,  since  the  effect  of  sweeping 
the  initial  shock  in  the  plt  nform  view  was  not  considered. 
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Ftoure  5-3:  CONTROL  SURFACE  GAP  GEOAffiTRY 


VI.  DESIGN  CHARTS 


This  section  contains  charts  that  provide  numerical  values  of  heating 
rates  to  the  geometric  elements  and  methods  discussed  in  Sections  II,  HI,  IV 
and  V.  As  stated  in  Section  II,  calculations  vere  made  assuming  the  flow  to 
behave  as  a  continuum  and  to  be  in  chemical  equlllbritim.  All  charts  showing 
altitude  were  developed  using  the  U.  S.  Standard  Atmosphere,  I962.  Other  gas 
properties  have  been  determined  using  references  1,  2,  and  26. 

A.  Method  of  Presentation 

Heating  rates  are  piesented  in  terms  of  equivalent  heat-transfer  coeffi* 
dents  and  reeoveiy  temperattires  defined  by 


h  = 


(6.1) 


where  4  is  the  local  heating  rate  Is  the  local  wall  temperatvire  and  T^^  Is 
the  recovery  temperature  defined  as 


T  =  IaW  _  W  (6.2) 

AW  cp  “  0.24 

where  1^^  is  the  recovery  enthalpy  and  Cp  is  the  specific  heat  at  constant 
pressure  for  a  perfect  gas.  Defining  h  and  in  this  manner  is  equivalent 
to  assuming  that  the  specific  heat  of  air  at  the  wall  temperature  is  always 


a,  =0.24 


Btu/lbm 


(6.3) 


While  in  reality  this  assumption  is  usually  not  true,  the  error  in  making 
this  approximation  is  small  because  of  the  large  difference  between  recovery 
enthalpy  and  wall  enthalpy  for  hypersonic  fli^t  conditions.  In  addition, 
the  assumption  is  always  conservative,  in  that  its  use  always  results  in 
slightly  higher  values  of  heating  rate.  Finally,  use  of  equation  (6.2) 
allowed  structural  designers  to  easily  obtain  equilibrium  wall  temperatures 
for  radiation  cooled  structures  using  the  simplified  thermal  balance  rela¬ 
tionship 


^RAD  ^  (6.4a) 

-  h  <Taw  -  Tw,,)  (6M) 

where  c  Is  the  surface  emissivity,  o  is  the  Stefan- Boltzmann  constant  and  4j^ 

is  the  heat  flux  lost  throu^  radiation  when  the  wall  temperature  has  reached 
eq^illlbrlum. 
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Heating  rates,  4>  were  calculated  assiaaing  a  wall  temperature  of 
2000  ®R.  This  temperature  was- chosen  because  Lt  represents  the  order  of 
temperature  one  can  expect  on  nonrefrectojTr  materials  during  the  most 
critical  part  of  lifting  re-entry.  Correcxion  charts  are  presented  that 
allow  heat-transfer  coefficients  to  be  obtained  for  other  wall  temperatures. 

B.  Reference  Conditions 

Heat  transfer  coefficients,  as  represented  by  equation  (6.l),  were 
calculated  for  all  of  the  elements  considered  in  Section  II,  III,  IV  and  V  at 
the  two  flight  conditions*  shown  in  Table  I. 


REFERENCE  FLIGHT  CONDITIONS 

Lacdnar 

Turbulent 

Altitude 

180.000  ft. 

20.VQ0  fps 

IliHRSNii 

T/3LE  I 

This  table  represents  the  design  conditions  most  critical  for  the  X-20. 

Distribution  functions  were  formed  for  each  geometry  by  normalizing  local 
values  of  heat  transfer  coefficient  by  a  reference  value  calculated  at  the 
appropriate  flight  condition  for  laminar  or  turbulent  flow.  Distribution 
functions  formed  in  this  way  are  relatively  weak  functions  of  altitude  below 
about  280,000  feet  and  velocity  between  8,000  to  26,000  fps.  For  laminar  flow, 
the  reference  value,  h  ,  chosen  was  the  stagnation  point  of  a  1- foot-radius 
hemisphere.  Laminar  heating  rates  for  all  of  the  geometric  elements  (i.e., 
sharp  plates,  sharp  and  blvint  cones  etc.)  considered  in  Sections  II,  III,  IV 
and  V  have  been  normalized  by  this  value  and  are  presented  in  terms  of 
attitude.  Other  altitude  and  velocity  combinations  are  then  obtained  by 
multiplying  distribution  functions  by  the  appropriate  reference  value. 


*For  laminar  flow  the  "altitude-velocity"  chosen  is  in  the  region  of  cheiilcal 
nonequilibrium  flow.  The  effect  of  chemical  nonequilibrium  on  heating  is 
not  taken  into  account  in  this  document.  However,  it  is  expected  to  be  small 
and  favorable. 


No  unique  vay  has  been  foxind  to  normalize  all  turbulent  beating  rates  by 
the  same  reference  value  and  keep  the  resulting  ratio  relatively  Insensitive 
to  altitude  and  velocity.  As  a  result^  In  turbulent  flov  tvo  stethcds  vere  used 
to  find  reference  values  for  normalizing  loced  heat  transfer  coefficients. 

The  first  Is  used  In  conjunction  with  sharp  flat  plates,  sharp  cones, 
sharp  delta  wings,  yawed  Infinite  cylinders  and  control  surface  gaps.  Refer¬ 
ence  values  correspond  to  turbulent  stagnation  line  heating  rates  for  a  60* 
swept  infinite  cylinder  with  a  1-foot  radius.  This  reference  value  was  used 
exclusively  during  the  design  of  the  X-<?0  and  was  designated  at  that  time  as 
the  "Turbulent  Reference  Heat  Transfer  Coefficient,"  h^.  This  nomenclature 
A.as  been  retained  in  this  report. 

Development  of  the  second  method  was  completed  during  vhe  compilation  of 
this  report.  It  ia  used  In  determining  reference  values  to  normalize  turbulent 
heating  rates  on  the  hemisphere,  unyawed  Infinite  cylinder,  blunt  cone  and 
blunt  delta  wing.  This  method  utilizes  the  relationship  between  laminar  and 
turbiilent  flow  inherent  In  the  p  (i  method.  This  expression,  developed  in 
Appendix  B,  Is 


0.185  R 

_ r _ 

0.332  [log  10  (Rr+  3000)] 2- 584 


(6.5) 


where  R  Is  the  reference  Reynolds  number,  defined  in  Appendix  B  as 
r 


R 


r 


Pr  ^*r  %  ^eq.L.o 


’‘eq.L 


(6.6) 


Equation  (6.5)  Is  shown  In  Appendix  B  to  be  valid  in  the  Reynolds  number 
remge  between  O.OO3  million  and  100  million. 

Absolute  values  of  laminar  or  turbulent  heat-transfer  coefficients  are 
then  obtained  at  altitudes  and  velocities  different  than  shown  in  Table  I 
through  the  use  of  the  appropriate  reference  value  using  t le  relationship 


where  the  subscript  ref  signifies  "o"  for  laminar  flow  and  "RT"  for  turb^ilent  flow. 


•♦Called  reference  because  of  Its  special  significance  to  the  P  j.  method. 
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C.  Reference  Condition  Charts 

1.  Recovery  Temperature 

The  equivsilent  recovery  temperature,  equation  (6.2)  has  been  eveduated  over 
a  range  of  velocities  between  8,000  and  26,000  fps.  It  Is  presented  In  figure 
6-1  as  a  function  of  free- stream  velocity  and  the  ratio  of  the  recovery  temper¬ 
ature  to  the  stream  total  temperature 


Tt  =  T„ 


2  c, 


P-.I 


(6.8) 


The  ratio,  T^/T^,  is  shown  over  the  range  O.85  (corresponding  to  flat 

plate  incompressible  laminar  flow  vedues  at  d  »  0;  reference  27*)  to  1.0,  the 
latter  value  being  at  the  stagnation  point. 


2.  Laminar  Reference  Heat -Transfer  Coefficient 

The  reference  heat-transfer  coefficient  for  laminar  heating  is  the  laminar 
heat  transfer  coefficient,  h^,  at  the  stagnation  point  of  a  1  foot  radius  hem¬ 
isphere.  It  is  shown  in  figure  6-2.  Wall  temperature  corrections  are  shown 
in  figure  6-3  as  a  function  of  free-streaa  Mach  nuaber,  M,p  .  Radius  corrections 
may  be  made  with  the  equation: 


ho 


.  R/ 


(6.9) 


where  h  is  the  stagnation  point  heat  transfer  coefficient  for  a  sphere  of 

radius  R  and  h  is  the  reference  laminar  heat  transfer  coefficient  for  a  sphere 
0 

of  radius  R  ■  1  foot, 
o 

3.  Turbulent  Reference  Heat  Transfer  Coefficient 

Absolute  values  of  the  reference  heat-transfer  coefficient  for  turbulent 
heating,  h^,  are  presented  in  figure  6-U.  These  curves  represent  the  turbulent 

stagnation-line  heat-transfer  coefficient  for  e  60°  swept,  1-foot-radius  in¬ 
finite  cylinder.  Wall  temperatui-e  corrections  are  shown  in  Figure  6-5. 

Figure  6-4  can  be  used  to  determine  turbulent  stagnation  line  heating  rates 
for  60®  swept  cylinders  of  radii  different  than  1  foot  by  using  the  correc¬ 
tion  factor  shown  in  figure  6-6.  This  correction  is,  however,  dependent  on 
the  reference  Reynolds  number  at  the  stagnation  line.  To  avoid  confusion, 
this  particular  evaluation  of  equation  (6.6)  was  designated,  R_,  and  is  pre¬ 
sented  in  figures  6-7  and  6-8  as  a  function  of  altitude  and  velocity. 
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D.  Geometric  Elements 
1.  Hemisphere 

The  laminar  design  curve  for  the  hemisphere  is  shown  in  figure  6-9.  Wal3 
temperature  correctirus  may  be  obtained  from  figure  6-3  and  a  radius  correc¬ 
tion  by  equation  (6.9) • 

Turbulent  design  curves  for  the  hemisphere  are  normalized  by  the  laminar 
value  as  described  above. 

Local  reference  Reynolds  number,  equation  (6.6),  at  the  location  of  maxi¬ 
mum  turbulent  heating  Is  presented  in  figure  6-10  in  terms  of  altitude  and 
velocity  for  a  hemisphere  of  R  *  1  foot.  The  reference  Reynolds  number  for 
a  hemisphere  of  radius  R  can  be  found  by  the  equation 


=  (R)  (Rr 
MAX  MAX  R»ift 


(6.10) 


where  R  is  obtained  from  fig<are  6-10.  The  location  of 

^MAX.R  =  1  foot 

maximum  turbulent  heating  is  shown  in  tenns  of  local  reference  Reynolds 
number  in  figure  6-11.  Local  heat  transfer  coefficients  normalized  to  the 
maximum  turbulent  values  are  presented  in  terms  of  the  local  angular  location 
in  figure  6-12.  A  slight  Mach  niunber  dependency  is  noted.  Finally,  the  level 
of  turbtilent  heating  is  determined  by  figure  6-13  which  relates  the  maximum 
txxrbulent  heat-transfer  coefficient,  the  laminar  stagnation-point 

value,  h_.  Figure  6-I3  is  based  on  a  correlation  that  is  a  function  only  of 

ex' 

the  local  reference  Reynolds  number,  R^  ,  provided  that  the  pressure  distri¬ 
bution  in  the  vicinity  of  maximum  turbu¥ent  heating  is  based  on  modified 
Newtonian  theory,  equation  (2.U).  For  local  reference  Reynolds  numbers  in 
the  range  0.01  million  to  1  million,  the  correlation  for 
pressed  by  the  equation. 


^MAX 

hsp 


0.0844Rx^AX-283 


(6.11) 


Vorticity  Interaction  Effects  on  Stagnation  Point  Heating 

Design  curve*  for  modifying  stagnation  point  heating  rates  for  vorticity 
interaction  are  shown  on  figures  6-l4  and  6-15. 

The  vorticity  interaction  parameter,  F  ,  is  obtained  using  the  relatlon- 

ship 

r  =  Fj  F2/VR80  (6.12) 
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where  and  are  defined  by  figure  6-l4  and 


Re 


o 


Mo 


(6.13) 


The  parameters  in  equation  (6.13)  were  defined  In  Section  II  as  follows: 
p  is  density, M  _  is  viscosity  evaluated  at  stagnation  conditions,  1  is 
stream  total  enthalpy,  and  R  is  the  nose  cap  radius. 

Figure  6-15  presents  the  ratio  of  heating  with  vortlclty  to  heating 
without  in  terms  of  the  vortlclty  interaction  parameter  f.  Both  the  Non- 
slmllar-Boundary- Layer  results  and  Ferrl's  results  (reference  8)  are  shown 
on  this  figure  allowing  the  reader  to  apply  either  method. 

2 .  Infinite  Cylinders 

a)  Swept  Infinite  Cylinders 

Design  curves  for  laminar  stagnation  line  heat-transfer  coefficients  are 
shown  in  figure  6-16.  Circumferential  heat-transfer  distributions  for 
laminar  flow  are  presented  in  figure  6‘-17  for  sweep  angles  between  zero  and 
80  deg?«es.  Turbulent  stagnation  line  heat-transfer  coefficients  presented 
in  figure  6-18  are  normalized  using  the  turbulent  reference  heat  transfer 
coefficient  h^^,.  Circumferential  heat -transfer  distributions  are  shown  in 

figure  6-19  over  the  same  range  of  sweep  angles  as  in  the  laminar  case  ex¬ 
cept  that  the  zero  degree  sweep  (\inswept  cylinder)  case  has  been  deleted 
from  this  curve. 

b)  Unswept  Infinite  Cylinders  in  Turbulent  Flow 

The  method  described  under  turbulent  heat-transfer  distribution  on  a  hem¬ 
isphere  also  applies  to  the  unswept  infinite  cylinder  althoutpi  the  charts  are 
different.  The  reference  Reynolds  number  at  the  location  of  maximum  turbu¬ 
lent  heating  on  the  unswept  cylinder  is  shown  in  figure  6-20  for  an  unswept 
Infinite  cylinder  1  foot  in  radius.  The  differences  between  the  infinite 
cylinder  and  the  hemisphere  are  primarily  due  to  changes  in  x  and  are  on 

the  order  of  2  to  2.5  times  larger  for  the  unswept  cylinder.  The  local 
.  »^erence  Reynolds  number  for  an  vinswept  cylinder  of  radius  R  may  be  found 
iag  equation  (6.10)  after  R  is  obtained  from  figure  6-20. 

*'mAX  R  =  1  foot 

Figure  6-21  shows  the  variation  in  the  location  of  maximum  turbulent  heat¬ 
ing  with  R  .  From  figure  6-22  the  turbulent  heating  distribution  can  be 
MAX 

obtained  in  terms  of  the  meximum  turbulent  heat-transfer  coefficient 


The  lower  curve  in  figure  6-I3  relates  the  maximum  turbulent  heat-trans¬ 
fer  coefficient  on  an  unswept  infinite  cylinder  to  the  stagnation  point  heat- 
ti^insfer  coefficient  of  a  hemisphere  having  the  same  radius  as  that  of  the 
cylinder.  For  reference  Reynolds  numbers  in  the  range  o.Ol  million  to  1 
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million  the  correlation  of  expressed  by  the  equation 


0.0613  R,  <'•283  (6.14) 

hsp  MAX 

Ccmparing  equations  (6.II)  and  (6.IU),  and  remembering  that  is  2  to 

2.5  times  larger  for  the  cylinder  than  for  the  hemisphere,  we  see  that  turbu¬ 
lent  heating  rates  on  an  unswept  cylinder  are  epproximetely  105^  less  than  on 
a  sphere  of  the  same  diameter. 

3.  Sharp  Flat  Plates 

Design  curves  for  laminar  and  turbulent  heating  on  the  sharp  flat  plate 
are  shown  in  figure  6-23.  The  laminar  heat -transfer  coefficient  hpp  has 

been  computed  for  x  =  1  foot.  For  turbulent  flow  the  heat -transfer  coeffi¬ 
cient  hjp  ^  is  presented  for  x  -  10  feet. 

The  wall  temperature  correction  is  shown  in  fi-?ire  6-24  and  applies  also 
to  sharp  unyawed  cones  and  sharp  delta  wings. 

The  distance  correction,  shown  in  figure  6-25  except  for  low  angles  of 
attack,  is  given  by: 


“m 

where  h  is  the  heat-transfer  coefficient  at  a  distance  x  euid  h  is  the  heat- 
X  m 

transfer  coefficient  at  distance  x  .  For  ].aminar  flow  n  =  O.5,  x  »  1.0 

n  m 

foot  and  h  »  fa_  ^  •  For  turbulent  flow  n  =  0.2,  x  *10  feet  and 
ni  m 

\  "  ^,T* 

The  exponent  on  the  ratio  (xj^/x)  is  a  variable  that  is  dependent  upon  the 
reference  Reynolds  number,  equation  (6.6). 


Figure  6-25  also  applies  to  unyewed  sharp  cone  and  sharp  delta  wing 
distance  corrections. 

k.  Unyawed  Cones 
a)  Sharp  cones 

Des.tgn  curves  for  laminar  and  turbulent  heating  on  imyawed  sharp  cone 
surfaces  are  shown  in  figure  6-26.  The  larolneu:  heat-transfer  coefficient 
hgp  on  the  sharp  cone  is  for  x*  1  foot  and  for  turbulent  flow  the  heat- 
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transfer  coefficient  h  is  for  x  =  10  feet.  Bie  wedl  tonperature  and 

oC^T 

distance  corrections  are  shown  in  figures  6-24  and  6-25  respectively. 

b )  Blunt  cones 

Design  curves  for  laminar  heating  on  an  unyawed  blunt  cone  having  a 
hemispherical  nose  cap  are  shown  in  figure  6-27*  A  noticeable  "bunrp"  is 
noted  in  this  curve  just  aft  of  the  nose  cap- cone  Junction.  Ihis  "bump" 
corresponds  to  the  region  for  which  perturbations  in  the  pressure  distribu¬ 
tion  have  occurred  6is  was  shown  in  Section  III,  figure  3-1*  lliese  effects 
damp  out  more  rapidly  at  large  semi-vertex  angles  of  the  cone. 

Wall  ten^jerature  effects  are  given  by  figure  6-3  for  the  stagnation 
point  region  of  the  hemispherical  nose  cap  and  figure  6-24  for  the  coni¬ 
cal  surface. 

Turbulent  heating  rates  are  presented  using  the  second  method  of  Sec¬ 
tion  VI-B.  The  local  reference  Reynolds  number  at  the  point  of  maximum 
turbulent  heating  on  a  1-foot-radius  nose  is  dei-ermined  using  figure  6-10.  For 
nose  cap  radii  different  than  1  foot,  modify  R  using  equation  (6.10). 

^MAX 

Figure  6-28  presents  local  reference  Reynolds  niunbers  at  a  distance. C 
from  the  stagnation  point.  The  ratio  defined  by  equation  (6.5)  for  blunt 
cones  is  shwm  in  figure  6-29*  Absolute  turbulent  heat-transfer  coeffi¬ 
cients  at  the  location  S  are  then  determined  from  the  relationship. 


-te)  (^) 


where  the  ratio  (h^^  defined  by  figure  6-29,  the  ratio 

(hRc  ij/hgp)  is  defined  by  figure  6-27,  defined  by  equation  (6.9) 

for  the  radius  used  to  determine  ,  and  h^  is  the  laminar  reference 
heat-treuisfer  coefficient  shown  in  ^e  6.2: 


5 .  Delta  Wings 
a)  Sharp  Delta  Wings 

Design  curves  for  laminar  and  turbulent  heating  on  a  70°  swept 
sharp  delta  wing  are  shown  in  figure  6-30.  The  laminar  heat-treinsfer  coeffi¬ 
cient  h^,  ^  is  for  x  =  1  foot,  and  the  turbulent  heat-transfer  coefficient, 
Dw,L 

hn(,  „  is  for  X  =  10  feet.  The  wall  temperature  and  distance  correction 

im  f  JL 

factors  are  shown  in  figures  6-24  and  6-25,  respectively.  Corrections  for 
the  effect  of  sweep  on  centerline  l-eating  can  be  made  using  the  curves  of 
figure  6-31*  Sp&nwise  heat  transfer  distribution  curves  Me  given  in 
figures  6-32  and  6-33  for  laminar  and  turbulent  flow,  respectively,  in  terms 
of  a  ratio  of  local  heat-trsuisfer  coefficient,  h,,  to  the  coefficient  at 
the  centerline,  h^_  The  effective  gamma  used  for  streamline  calculations 
was  equeil  to  1.2.' ~ 
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b)  Blunt  Delta  Wlnga 

Laminar  heating  design  curves  are  presented  in  figure  6-3^  for  obtaining 
values  on  the  centerline  of  a  73"*  swept  blunt  delta  wing.  As  noted  in  Sec¬ 
tion  nose  blimtness  effects  were  neglected  and  are  not  Included  on  this 
figure.  Correction  factors  for  sweep  effects  on  centerline  heating  are  pre¬ 
sented  in  figure  6-35* 

Spanwise  laminar  heetlng  distributions  are  shown  in  figures  6-36  and  6-37* 
As  noted  in  Section  III,  these  are  semi-empirical  distributions  that  have 
been  determined  using  data  (reference  I6)  measured  in  proximity  to  equiva¬ 
lent  centerline  locations  of  S/R  =  I7  and  2Tj  in  conjimetion  with  theoretical 

centerline  heating  rates  calculated  using  the  fi  U  method,  Api)endix  B,  and 

r  r 

semi-empirical  leading-edge  stagnation-line  locations  determined  using  the 
faired  curve  shown  in  Section  III-B-2. 

Leading-edge  sweep  angle  influences  spanwise  heating  diiStributlons. 

For  the  range  of  sweep  angles  of  most  interest  in  hypersonic  flight  (usually 
between  65®  and  80® )  a  first  order  approximation  can  be  made  using  informa- 
v.ion  available  on  design  charts  appearing  in  this  document.  The  following 
procedure  is  suggested: 

(1)  At  the  Leading  Edge 

Establish  a  new  stagnation  line  location  using  figure  6-38  and  effec¬ 
tive  sweep  angle  using  figure  6-39*  Determine  the  percentage  difference 
between  laminar  stagnation  line  heating  rates  for  swept  infinite  cylinders, 
figure  6-17,  having  effective  sweep  angles  corresponding  to  geometric  sweep 
angles  of  73"  and  the  one  being  investigated.  Shift  the  stagnation  line 
location  shown  on  figure  6-36  or  6-37  to  the  new  value  and  modify  the  stag¬ 
nation  line  heating  rate  ratio  by  the  percentage  difference  in  heating  rates 
previously  determined  for  swept  infinite  cylinders.  Constiruct  a  new  dis¬ 
tribution  curve  in  proximity  to  the  stagnation  line  using  the  slopes  shown 
on  figure  6-36  or  6-37* 

(2)  At  the  Centerline 

Modify  the  values  at  the  centerline  using  the  sweep  correction  curves 
of  figure  6~35*  Construct  a  new  distribution  curve  similar  t>o  the  one  shown 
on  figure  6-36  or  6-37  depending  on  S/R. 

For  leading-edge  sweep  angles  less  than  73*  centerline  heating  rates 
decrease  whereas  the  stagnation-line  heating  rates  increase.  For  leading- 
edge  sweep  angles  greater  than  73"  centerline  heating  rates  Increase  but 
stagnation  line  heating  rates  decrease.  For  either  situation  modify  the 
leading  edge  distribution  cxurve  of  figure  6-36  or  6-37  to  intersect  the  new 
lower-surface  heating  distribution  curve. 

To  establish  laminar  spanwise  heating  distributions  for  centerline  values 


^<Thi6  sweep  angle  was  chosen  because  the  bulk  of  the  delta  wing  models  used 
during  the  X-20  test  program  had  this  sweep,  and  the  data  obtained  from  that 
pi’Ogram  was  used  to  substantiate  the  method  and  provide  spanwise  distributions. 
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of  S/r  different  than  those  shown  on  figures  6-36  and  6-37  the  following 
procedure  is  suggested.  The  value  at  the  centerline  for  the  s/R  of  interest 
is  established  using  figure  6-3*^.  Construct  a  new  lower  surface  heating 
distribution  curve  parallel  to  the  one  shown  in  figai'e  6-36  or  6-37  (using 
the  figure  having  an  S/R  closest  to  the  one  desired)  until  it  intersects 
the  leading  edge  distribution  curve  shown. 


Turbulent  heating  rates  to  the  centerline  of  the  delta  wing  can  be 
obtained  using  the  second  reference  method  which  is  slmlleur  to  the  procedure 
used  with  a  blunt  cone.  The  reference  Reynolds  number  at  the  location  of 
maximum  turbulent  heating  on  the  hemispuerical  nose  is  obtained  from 
figure  6-10  and  is  used  in  conjunction  with  figure  6-40  to  determine  tnr' 
local  reference  Reynolds  ntmber  on  the  lower  surface  centerline  at  a  dis¬ 
tance  S  from  the  stagnation  point.  The  relation  between  the  turbulent  and 
laminar  heat  transfer  coefficient  is  determined  using  figure  6-29*  The 
laminar  surface  heat  transfer  coefficient  is  related  to  the  laminar  stagna¬ 
tion  point  heat  transfer  coefficient  by  figure  6-34.  Finally  the  turbulent 
heat  transfer  coefficient  h^  is  calculated  by  the  relationship 


h,j.  =  O'sp^o)  **0  (6.17) 


where  h  is  obtained  from  figure  6-2.  The  value  for  (h_v^h  )  is  determined 
o  SF  o 

using  equation  (6.9) ♦ 

6.  Deflected  Flaps 

Design  curves  for  maximum  laminar  heating  on  deflected  flap  surfaces 
attached  to  sharp  flat  plates  are  shown  in  figure  6-^4l.  That  portion 
of  the  curves  that  is  dashed  is  outside  the  range  for  which  the  "upper 
limit"  theory  discussed  in  Section  HI  is  applicable.  Use  of  these  curves 
within  the  dashed  range  should  be  made  with  caution.  For  a  positive  deflec¬ 
tion  the  maximum  heat-transfer  coefficient  h^,  may  occur  anywhere  on  the  flap; 

generally,  if  flow  separation  has  occurred,  hp  is  near  the  point  of  flow  re¬ 
attachment.  The  heat-transfer  coefficient  h^  is  the  laminar  heat  transfer 

coefficient  calculated  just  ahead  of  the  hinge  line  of  the  flap  with  the 
assumption  that  no  flow  separation  has  been  caused  by  the  deflected  flap. 

The  value  of  the  hinge  line  heat-trsuisfer  coefficient  hg^^  may  be  evaluated 

by  sharp  flet-plate  heat-transfer  methods,  figures  6-23  through  6-25* 

For  a  negative  deflection,  hp  is  an  approximate  heat  transfer  coeffi¬ 
cient  at  the  hinge  line  assuming  that  the  expansion  has  already  occurred. 

Aft  of  the  hinge  line  the  heating  decay  follows  the  relation  on  the  follow¬ 
ing  page 
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h  ~  (l/x)” 


(6.18) 


where  n  *  0.5  for  laminar  flow  and  0.2  for  turbulent  flow. 

Maximum  turbulent  heating  rates  vere  normalized  using  the  laminar  values 
for  the  same  local  boundary-layer  conditions.  This  methodt  la  consistent 
with  the  second  turbxxlent  reference  method  described  in  Section  VI-B.  Ved- 
ues  are  preserted  in  figure  6-42. 

7.  Surface  Conditions 

a)  Waves 

Heat  transfer  design  curves  for  a  particular  type  of  roughness^  shallow 
surface  waves,  are  shown  in  figures  6-43  *uad  6-44,  in  teims  of  a  wave  height 
parameter,  E/5*.  The  displacement  thiclmess  5*  may  be  eveduated  by  the 
use  of  the  displacement-thickness  heating  parameter  h6*  shown  in  figure 
6-45  and  a  vaj.ue  for  the  heat -transfer  coefficient  h  for  a  smooth  surface. 

The  latter  is  found  using  figure  6-23.  The  derivation  of  the  product  h8* 
can  be  found  in  Appendix  D. 

b)  Leakage  (Mass  Removal  Through  an  Orifice  at  the  Surface) 

Leakage  effects  on  sharp  flat  plate  laminar  heating  at  all  angles  of 
attack  in  terms  of  the  distance  downstream  of  a  leak  are  shown  on  figure  6-46. 
The  heat-transfer  parameter  is  shown  in  terms  of  a  nondimensional  heat-trans¬ 
fer  ratio  and  a  nondimensional  mass  flow  ratio,  L,  as  defined  in  Section  IV. 

8.  Control  Surface  Gaps 

Design  curves  for  laminar  and  turbulent  heating  at  the  most  critical 
location  on  a  fin  for  representative  tip-fin-elevon  control -surface  gaps 
(figure  5”3»  Section  V)  are  shown  in  figure  6-47.  This  location  is  90*  from 
the  tip-fin  shoulder.  For  the  gapj  nhown  in  figure  5”3  this  is  the  point 
of  minimum  fln-elevon  separation.  Meximum  heating  rates  do  not  occur  at 
this  location.  However,  maximum  temperatures  do  occur  if  the  structure  is 
radiation  cooled,  because  the  radiation  view  factor  is  at  its  minimum.  For 
all  of  the  gap  shapes  shown  in  figure  5-3  one  design  curve  is  applicable  at 
each  of  the  two  flight  reference  conditions.  The  leading  edge  radius  of  the 
fin  shown  in  this  figure  is  equal  to  1  foot.  For  radii  different  from 
1  foot,  corrections  can  be  made  to  the  heating  rates  using  figure  6-6  for 
tvirbulent  flow  and  equation  (6.9)  for  laminar  flow  where  h  is  the  reference 

heat- transfer  coefficient  for  a  radius  of  R  =  1  foot,  and  R  is  the  fin 
radius  in  feet. 


fUse  of  this  method  for  this  combination  of  two-dimensional  bodies  consid¬ 
erably  simplified  the  extrapolatlon-to-fllght  methods  for  this  geometry  that 
are  presented  in  Section  VTI. 
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9*  Canopy  and  Fin  Leading  Edges 

Laminar  and  turbulent  flow  design  curves  for  predicting  canopy  stagna¬ 
tion  line  heating  rates  are  shown  In  figure  6-48.  The  curves  shown  were  de- 
veloi>ed  using  the  method  described  in  Section  V-B.  Heat -transfer-coefficient 
ratios  are  shown  in  terms  of  the  time  sweep  angle  (a  +  \) . 

Values  for  h  are  obtained  from  fig*are  6-2  and  hj^^  from  figure  6-4. 

The  radius  of  the  leading  edge  shown  is  l  foot.  For  radii  different  from 
thls^  a  correction  factor  can  be  obtained  using  figure  6-6  for  turbulent  flow  or 
equation  (6.9)  for  laminar  flow. 

As  noted  in  Section  V,  heating  rates  obtained  using  the  above  method 
were  also  applied  to  tip-fin  leading  edges  during  the  X-20  program.  Figure 
6-48  can  be  used  directly  for  this  geometry  by  allowing  the  indepehdent  var¬ 
iable,  (a  +  X),to  represent  the  effective  sweep  angle,  A  of  "the  fin 
leading  edge  measured  with  respect  to  the  vertical  axis  Sormal  to  the  free 
stream  velocity  vector.  This  angle  is  defined  by. 


cos  Ag££  =  cos  a  ^cos  cos  (0  +  /5)  cos  A  +  sin  Sgj^  sin  (||)  +  j8)^ 
+  sin  a  cos  ^ 


(6.19) 


where  A  is  the  fin  leading-edge  geometric  sweep  angle  and  is  the  location 
of  the  stagnation  line  as  defined  by: 


tan  6, 


sin  (<f)  +  j3) 


SL  cos  (</)  +  P)  cos  A  +  tan  a  sin  A 


(6.20) 


The  remaining  angles  are  defined  in  the  following  sketch. 
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Figure  6-3:  WALL  TEMPERATURE  CORRECTION  FACTOR  FOR 

LAMINAR  REFERENCE  HEAT  TRANSFER  COEFFICIENT 
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Figure  6-9:  HEMISPHERE  LAMINAR  HEATING  DISTRIBUTION 
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Figure  6-13:  HEMISPHERE  AND  UNSWEPT  INFINITE  CYLINDER 
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Figure  6-15:  STAGNATION  POINT  VORTICITY  INTERACTION  EFFECT 
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Figure  6-17:  SWEPT  INFINITE  CYLINDER  LAMINAR 

HEATING  DISTRIBUTION 
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AT  LOCATION  OF  MAXIMUM  TURBULENT  HEATING 
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Figure  6”23:  SHARP  FLAT  PLATE  HEATING 
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Figure  6-26:  UNYAWED  SHARP  CONE  HEATING 
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6-27:  UNYAWED  BLUNT  CONE  LAMINAR  HEATING 
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Figure  6-30;  SHARP  OaiA  WING  LOWER  SURFACE  CENTERLINE  HEATING 
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HEAT  TRANSFER  COEFFICIENT  RATIO 


Figure  6-31:  EFFECT  OF  SWEEP  ON  SHARP  DELTA  WING 
LOWER  SURFACE  CENTERLINE  HEATING 
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Figure  6-32:  SHARP  DELTA  WING  SPANWISE  LAMINAR  HEATING  PARAMETER 
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SPANWISE TURBULENT  HEATING  PARAMETER 
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Figure  6-35;  EFFECT  OF  SWEEP  ON  BLUNT  DELTA  WING 
LOWER  SURFACE  CENTERLINE  HEATING 
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Figure  6-36:  BLUNT  DaiA  WING  LOWER  SURFACE 
SPANWISE  LAMINAR  HEATING 
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Figure  6-37:  BLUNT  DELTA  WING  LOMR  SURFACE 
SPANWISE  LAMINAR  HEATING 
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Figure  6-43:  SURFACE  WAVE  LAMINAR  HEATING 
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Figure  6-44:  SURFACE  WAVE  LAMINAR  HEATING 
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Figure  6-45:  SHARP  FLAT  PUTE  LAMINAR  DISPLACEMENT 
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Figure  6-46:  EFEa  OF  LEAKAGE  ON  LAMINAR  HEAT  TRANSFER 
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igure  6-47:  CONTROL  SURFACE  GAP  HEATING  (TIP  FIN  -  ELEVON) 
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VII.  USE  OF  GROUND-FACILITY  DATA  IN  DESIGN 


Groimd  facilities  presently  being  used  in  the  design  of  hypersonic 
flight  vehicles  cannot  sijnultaneously  duplicate  velocities  and  Mach  numbers 
associated  with  re-entry.  Ground-test  data,  however,  are  required  to  obtain 
detailed  aerodynamic -heating  distributions  oh  complex  configurations  for 
which  completely  analytic  solutions  are  not  possible.  To  apply  ground -test 
data  in  actual  design,  it  is  necessary  to  correct  for  the  effects  of  any 
differences  that  exist  between  ground  facility  environment  and  flight  con¬ 
ditions.  Several  examples  of  effects  that  ml^t  be  expected  are  discussed  in 
reference  2d;  a  few  specific  examples  are  discussed  in  this  section. 

Also  Included  in  this  section  is  a  description  of  methods  for  correcting 
wind  tunnel  data  to  flight  conditions.  Charts  are  presented  that  provide 
the  required  correction  factors  for  those  geometric  elements  which  could  be 
investigated  analytically.  The  purpose  of  these  charts  is  to  serve  as  a 
guide  in  extrapolating  ground-test  data  to  flight  conditions  for  complex 
geometries.  Since  these  charts  provide  a  means  of  extrapolating  wind-tunnel 
data  to  fll^t  conditions,  they  will  be  referred  to  as  "extrapolation-factor 
charts . " 

All  charts  presented  in  this  section  are  based  on  the  assumptions  dis¬ 
cussed  in  Section  II  and  do  not  account  for  chemical  nonequlllbrlum  effects. 
In  addition,  chemical  equilibrium  is  assumed  to  exist  in  ground -facility  test 
sections  both  in  the  free  stream  and  within  the  body  shock  wave. 


A.  Real-Gas  Effects  on  Shock-Layer  Properties 

1.  Pressures 

Some  effects  of  the  differences  between  ground-test  and  flight-test 
conditions  on  inviscld-flow  shock-layer  properties  can  be  determined  from 
oblique-shock  calculations  for  a  two-dimensional  flat  plate.  Real- gas  effects 
on  sharp-flat -plate  pressure  coefficients  are  shown  in  figure  7"1»  Pressure 
coefficients  in  ground  facilities  having  free-stream  trmperatures  on  the 
order  of  100  ®F,  are  compared  to  pressure  coefficients  that  would  be  obtained 
at  an  altitude  of  240,000  feet  at  the  same  Mach  nianbers,  Viscous-inviscld 
interaction  and  separated-flow  effects  are  not  considered  in  Figure  7“1* 

At  high  angles  of  attack,  pressure  coefficients  for  flight  are  seen  to 
be  lower  than  those  in  ground  facilities  at  the  same  Mach  number.  This 
reduction  in  pressure  is  primarily  caused  by  the  higher  enthalpy  achieved  in 
fll^t.  At  the  higher  enthalpies  of  fll^t  real-g..v3  effects  increase  the 
density  Jump  across  shock  waves  which  thereby  tend  to  lie  closer  to  the  body. 
Since  the  pressure  experienced  by  the  body  is  the  reaction  to  flow  deflections 
caused  by  the  shock,  the  smaller  shock  envelope  occurring  in  fli^t  leads  to 
lower  pressures  on  the  body. 

Of  course,  for  highly  swept  shock  waves,  the  downstream  temperature  may 
be  relatively  low,  even  at  hypersonic  speeds,  and  so  fail  to  induce  any 
appreciable  real-gas  behavior.  This  explains  the  absence  of  any  apparent 
differences  in  pressure  coefficients  for  angles  of  attack  less  than  about  10 
degrees. 
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The  conclusions  drawn  from  figure  7-1  are  not  applicable  to  ail  geom- 
tries,  however.  An  opposite  trend  can  occur  on  a  trailing-edge  flap  deflected 
into  the  stream,  as  illustrated  in  figure  7“2.  At  a  10-degree  angle  of  attack, 
flap  pressure  coefficients  for  wind-tunnel  conditions  are  shown  to  agree  with 
those  for  fli^t  conditions.  When  the  flat  plate  is  inclined  at  larger  angles 
of  attack,  however,  real-gas  effects  become  noticeable  and  are  often  opposite 
of  those  for  a  flat  plate.  The  departure  from  the  flat-plate  trend  reflects 
the  effect  of  shock  layer  flow  properties  on  the  pressure  jump  at  the  flap. 

The  effective  free  stream  for  the  flap  is  of  course  the  flow  wj.thin  the  plate 
shock  J.ayer.  Although  the  static  pressure  of  the  local  flow  is  reduced  by 
real-gas  effects,  the  aforementioned  Increase  in  the  local  density  leads  to 
a  corresponding  increase  in  the  local  dynamic  pressure.  As  a  result,  the 
pressure  jump  at  the  flap  is  increased,  causing  the  trends  shown  in  figure 
7-2. 


The  combined  effect  of  Mach  nvimber  and  enthalpy  differences  on  deflected 
flaps  is  illustrated  in  figure  7"3.  The  curves  shown  are  for  two  specific 
wind  tunnels  and  particular  flight  conditions,  as  noted.  Negative  flap 
deflections  are  seen  to  cause  lower  pressure  coefficients  in  flight  than  in 
the  wind  tvinnels,  while  positive  deflections  lead  to  higher  values.  Thus 
it  is  seen  that  the  effect  of  flap  deflection  is  always  larger  in  flight  than 
in  the  wind  tunnel.  Although  figure  7-3  is  drawn  specifically  for  flat- 
platf  flap  combinations,  similar  trends  are  to  be  expected  for  flaps  attached 
to  sharp  cones  or  delta  wings. 

Real-gas  effects  on  stagnation  point  pressure  are  relatively  small,  as 
shown  in  figure  7-^'* 

2 .  Streamlines 

The  lack  of  enthalpy  simulation  in  groxmd  facilities  has  an  effect  on 
streamline  patiems  on  swept  infinite  cylinders.  The  differences  in  the 
local  streamline  angle  at  each  circumferential  location  on  the  surface  of  a 
50®  swept  infinite  cylinder  are  illustrated  in  figure  7-5  as  a  function  of 
total  enthalpy  ratio.  The  local  streamline  angles  were  calculated  using 
methods  described  in  Section  II. 

3.  Boundary  Layer  Properties 

In  flight,  dissociation  may  occur  within  the  boundary  layer,  altering  the 
temperature,  density,  and  velocity  profiles.  A  comparison  of  flight  and 
ground-facility  gas-temperature  profiles  similar  to  those  presented  in 
reference  28  are  shown  in  figure  7-6  for  a  flat  plate  at  a  20-degree  angle 
of  attack.  The  profiles  are  shown  as  a  function  of  tj,  the  nondlmensional 
distance  normal  to  the  wall.  Considerable  differences  are  seen  not  only  in 
the  levels  of  temperature  but  also  in  the  ratio  of  peak  temperature  to  boun¬ 
dary-layer-edge  temperature.  Part  of  the  difference  is  due  to  the  relatively 
higher  wall  temperature  in  the  groimd  facilities  and  part  to  dissociation 
in  the  boundary  layer  in  the  flight  case. 
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Displacement  thickness  is  significantly  affected  by  the  changes  in  the 
density  and  velocity  profiles  in  the  boundary  layer.  The  effect  on  displace¬ 
ment  thickness  is  shown  in  figure  7"7*  For  the  conditions  indicated,  the 
displacement  thickness  in  the  wind  tunnel  is  approximately  50^  higher  than 
in  flight. 

4 .  Heating  Rates 

The  resultant  effect  of  the  changes  in  boxmdary-iayer  characteristics  on 
heat  transfer  shown  in  figures  7“6  and  7"7  can  be  illustrated  by  comparing 
a  normalized  surface  heat-transfer  coefficient  at  ground-facility  conditions 
to  a  similar  ratio  at  flight  conditions.  The  resulting  I'atio,  denoted  by  f 
and  specifically  defined  in  figure  has  been  calculated  for  several 

geometries.  Values  of  ^  for  a  flat  plate  are  shown  in  figure  7"8  as  a 
fimction  of  angle  of  attack  for  two  different  facilities.  Free-stream 
temperatiires  in  these  two  facilities  are  nearly  eq^ual  but  the  total  enthalpy 
and  the  Mach  Number  are  higher  in  the  shock  tunnel. 

Similar  information  for  a  deflected  flap  is  shown  in  figure  7“9»  This 
figure  illustrates  differences  in  the  heat-transfer  parameter  for  a  flap 
that  is  attached  to  a  sharp  flat  plate  at  a  20-degree  angle  of  attack. 

Partial  simulation  effects  on  stagnation-line  heat  transfer  for  infinite- 
swept  cylinders  are  Illustrated  in  figure  7~10  fis  a  function  of  sweep  angle 
for  two  different  ground  facilities.  Curves  shown  Indicate  that  at  large 
sweep  angles  opposite  trends  in  the  heat-transfer  parameter  can  occur  in 
tunnels  having  different  energy  levels  and  Mach  mmibers. 

B.  Extrapolation  of  Ground-Facility  Heat- Transfer  Test  Data  to  Flight 
Conditions 

Calculations  of  the  type  just  described  have  been  made  in  order  to  provide 
extrapolation  factors  that  allow  the  correction  of  experimental  ground -test 
heat-transfer  data  to  flight  conditions.  The  factors  presented  are  calculated 
for  the  basic  shapes  discussed  in  Section  II,  but  should  serve  as  guides  to 
the  collection  of  data  for  more  complex  shapes. 

The  previously  defined  heat-transfer  parameter  ^  can  be  used  to  correct 
ground  facility  data  to  flight  conditions.  For  example: 

•■l.FLT”  (^)VFLT(t^)^^ 

-T.FLT  -  (^)  -RT.FLT 

For  several  of  the  basic  shapes  it  was  found  that  the  major  effect  of 
Mach  number  on  f  is  due  to  its  effect  on  pressure.  Accordingly,  a  compress- 
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ibllity  parameter  A  has  been  defined  as  follows : 


A  = 


(7.2) 


where  P  is  the 


local  surface  pressure,  P_  is  the  stagnation-point  pressure 

^2 


on  a  reference  body,  the  subscript  WT  denotes  ground -facility  stream 
conditions  and  the  subscript  FLT  refers  to  the  fli^t  reference  condition. 
The  subscript  f  indicates  a  sphere  in  laminar  flow  or  a  60®  swept  cylinder 
in  tvirbulent  flow.  The  exponent  n  is  (by  definition)  O.5  in  laminar  flow 
and  0.8  in  turbulent  flow.  Using  A  ,  the  correction  of  wind  tunnel  data  to 
flight  conditions  is  made  for  example  as  follows: 


^L.FLT  "  ^o 


A^ 

Ol 


(4)©., 


(7.3) 


Curves  of  A  and  ?/A  for  several  basic  shapes  are  presented  in  this 
section.  For  control  surfaces  extrapolation  to  fli^t  conditions  is  made 
with  a  single  chart  without  the  use  of  a  compressibility  parameter.  The 
basis  for  these  charts  is  discussed  in  detail  in  the  following  sections. 


1.  Compressibility  Parameter  (A) 

Compressibllity-par&meter  charts  applicable  to  hypersonic  ground  facil¬ 
ities  are  presented  in  figures  7“11  through  7"l8*  Free  stream  temperatures 
from  120® R  down  to  the  llqulfaction  limit  (which  ranges  from  60  ®R  to  90  “R 
depending  on  the  pressure)  were  considered.  The  effects  of  enthalpy  differ¬ 
ences  on  A  for  the  range  of  ground  facilities  indicated  above  is  small  and 
therefore  not  shown. 

2 .  Extrapolation  Factor  ( C  / A ) 

Extrapolation  factors  corresponding  to  figures  7-II  through  7-18  are 
presented  in  figures  7-19  throu^  7”29.  Also  included  in  the  grov.p  is  the 
extrapolation  factor  for  the  hemisphere,  figure  7“25»  For  the  hemisphere, 

A  is  unity  and  C  /  A  =  C  •  This  is  consistent  with  the  assumption  made  in 
Section  II  that  the  pressure  distribution  is  constant  for  all  conditions. 

The  extrapolation  factors  were  computed  for  the  same  range  of  conditions 
as  for  A  .  The  effect  of  Mach  number  differences  on  C/A  for  the  range  of 
conditions  studied  is  small  and  therefore  not  shown. 
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3.  Correction  Factors  for  Deflected  Flaps 


The  heat-transfer  method  employed  in  the  analysis  of  deflected  f.^aps  is 
approximate  and  the  extrapolation  procedvire  described  in  the  above  paragraphs 
does  not  apply.  Instead  the  extrapolation  is  made  with  the  equation 


The  extrapolation  factor  ?  for  specific  wind  tunnel  Mach  numbers  is 
presented  in  figures  7-3O  and  7-3I. 
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FLAT  PLATE  PRESSURE  COEFFICIENT 
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STAGNATION  POINT  PRESSURE  RATIO 


STREAMLINE  ANGLE,  t —DEGREES 


ANGULAR  LOCATION,  ^—DEGREES 


Figure  7-5:  REAL  GAS  EFECT  ON  STREAMLINE  ANGLE  FOR 

INFINITE  SWEPT  CYLINDER 
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DISPLACEMENT  THICKNESS 


HEAT  TRANSFER  PARAMETER,  C 
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FLAP  DEFLECTION  ANGLE,, ^-DEGREES 

Figure  7-%  FLAP  SURFACE  HEATING  COMPARISON 


CYLINDER  HEATING  COMPARISON 


LAMINAR  COMPRESSIBILITY  PARAMETER 


TURBULENT  COMPRESSIBILITY  PARAMET 


Figure  7-12:  SHARP  FLAT  PLATE  TURBULENT  COMPRESSIBILITY 

PARAMETER 
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LAMINAR  COMPRESSIBILITY  PARAMETER 


Figure  7-13;  UNYAWED  SHARP  CONE  LAMINAR  COMPRESSIBILITY 

PARAMETER 
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TURBULENT  COMPRESSIBILITY  PARAMETER,  ^SC,T 


Figure  7-i4:  UNYAWED  SHARP  CONE  TURBULENT  COMPRESSIBILITY 

PARAAAETER 
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TURBULENT  COMPRESSIBILITY  PARAMETER,  j  LAMINAR  COMPRESSIBILITY  PARAMETER 


LAMINAR  COMPRESSIBILITY  PARAMETER 
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Figure  7-16:  SWEPT  INFINITE  CYLINDER  LAMINAR 
COMPRESSIBILITY  PARAMETER 


TURBULENT  COMPRESSIBILITY  PARAMETER 


TURBULENT  COMPRESSIBILITY  PARAMETER,  LAMINAR  COMPRESSIBILITY  PARAMETER, 
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Figure  7-18j  CONTROL  SURFACE  GAP  COMPRESSIBILITY  PARAMETERS 
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Figure  7-19:  SttARP  FtAT  PLATE  LAMINAR  HEATING  EXTRAPOLATION  FACTOR 
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Figure  7-20:  SHARP  FLAT  PLATE  TURBULENT  HEATING  EXTRAPOLATION  FACTOR 
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Figure  7-22:  UNYAWED  SHARP  CONE  TURBULENT  HEATING  EXTRAPOLATION  FACTOR 
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Figure  7-23:  SHARP  DELTA  WING  LAMINAR  HEATING  EXTRAPOLATION  FACTOR 


EXTRAPOLATION  FACTOR 


#  A  -  COMPRESSIBILITY  PARAMETER 

.  T  -  turbulent  center 

LINE  HEAT  TRANSFER 
COEFFICIENT  ON  70'-^ 
SWEPT  SHARP  DELTA 
WING  AT  X  =  10  FT 
«  h  -  TURBULENT  REFERENCE 
HEAT  TRANSFER 
COEFFICIENT 

•  TURBULENT  REFERENCE 
CONDITION 

VEL  =  18,700  FPS 
ALT  =  180.000  FT 


TOTAL  ENTHALPY  RATIO,  ^ 


Figure  7-24:  SHARP  DELTA  WING  TURBULENT  HEATING 

EXTRAPOLATION  FACTOR 
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Figure  7-26:  SWEPT  INFINITE  CYLINDER  LAMINAR  KEATING  EXTRAPOLATION  FACTOR 
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Figure  7-27:  SWEPT  INFINITE  CYLINDER  TURBULENT  HEATING  EXTRAPOLATION  FACTOR 
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Figure  7-29:  COfsTTROL  SURFACE  GAP  TURBULENT  HEATING  EXTRAPOLATION  FACTOR 
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Figure  7-30:  FLAP  SURFACE  HEATING  EXTRAPOLATION  FACTOR 
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Figure  7-31;  FLAP  SURFACE  HEATING  EXTRAPOLATION  FACTOR 
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APPENDK  A 


SIMPLIFIED  EQUATIONS 


Simplified  equations  have  been  developed  for  modified  Newtonian  pressure 
coefficients,  figure  2-h,  reference  heat  transfer  coefficients,  figures  6-2 
and  6-U,  reference  Reynolds  number,  figures  6-7  and  6-8,  and  the  local 
reference  Reynolds  numbers  at  the  location  of  maximum  turbulent  heating,^ 
figure  6-10  and  6-20.  Heat-transfer  coefficients  have  the  units,  Btu/fu  -sec-*R. 

The  pressure  coefficient  equatio'n  was  cieveloped  for  Y  =  1.1  and  does  not 
apply  to  wind  tunnel  conditions.  All  other  equations  for  both  wind-tunnel 
and  flight  conditions  were  developed  from  numerical  results  using  the  Pj, 
method,  Appendix  B.  The  I962  U.  S.  Standard  Atmosphere  defines  the 
flight  conditions. 

The  accuracy  of  the  equations  in  general  is  +  or  better.  Simplified 
equations  applicable  to  wind  tunnels  have  been  developed  in  terms  of  free- 
stream  conditions  Instead  of  the  more  desirable  stagnation  chamber  conditions 
because  difficulty  was  experienced  in  obtaining  a  simple  form  of  the  equation. 
Free- stream  conditions,  however,  can  be  accurately  and  easily  converted  to 
stagnation  chamber  conditions  using  figure  A1  and  equation  (AI5). 


A,  Modified  Newtonian  Pressure  Coefficient 

For  a  ratio  of  specific  heats  equal  to  1.1  the  obarte  on  figure  2*4  can  be 
approxliaatad  by  the  following  equation 


P 

sin^S 


mm 

li 

1/2" 

1.05  + 

1.1025  +  -"5 - oT 

,mt 

sin*^  0 

,  -  1.278 

Mm 


(Al) 


B.  Reference  Heat  Transfer  Coefficient 
1.  Laminar 

The  reference  laminar  heat-transfer  coefficient  h^  can  be  approximated 
for  wind-tunnel  and  flight  conditions  by  the  following  equations. 


(a)  Wind  Tunnel 
V.  0.004 
‘‘o  •  rO.5 


{■ 


\0.5 


l.io  „0.12 

M  T 


(A2) 
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(b)  night 


2 .  Turbulent 

The  reference  turbulent  heat-transfer  coefficient,  h^^^,  can  be  approxi¬ 
mated  for  wind-tunnel  and  flight  conditions  by  the  following  equations. 

(a)  Wind  Tunn^ 


0.0415 

/  p  \o.a  i.*»6  -0.114 

** 

(b)  Flight 

hj^T  - 

( 

0.2  1 
R  \ 

^’ato  / 

m 

2:hIL 

\  ®ATM- 

>0.7®/  V  \l.54 

)  (iS^) 

(A4) 


(A5) 


Equations  {Ak)  and  (A5)  are  applicable  only  for  flow  conditions  where  tur¬ 
bulence  is  expected,  and  where  the  cylinder  radius  is  near  1  foot.  If  the 

radius  is  much  different  from  1  foot  the  (i/R  “^)  term  creates  lar^e  errors. 

An  alternate  method  has  been  developed  that  uses  the  above  equations  ns  a  base 
and  increases  the  accuracy  to  within  +  O.S^.  This  method  follows: 

1.  Calculate  h^^  with  R  =  1  foot  using  equation  Ah  or  A5. 

2.  Calculate  a  reference  Reynolds  number,  using  the  method  des¬ 
cribed  in  G.  below. 

3.  Correct  h_n  to  the  desired  radius  by  using  figure  6-6  of  Section  VI. 

nT 

C.  Reference  Reynolds  Number  (Rj^) 

The  reference  Reynolds  number,  R^^,  is  approximated  for  wind  tunnel  and 
flight  conditions  by  the  following  equations: 


(a)  Wind  Tunnel 


P  \  0.937  -0.679 

R„  =  2.65  X  lO^^f-.^  )  M  T  R 


(A6) 


(6)  Flight 


V»  <  13,000  ITS 


R  =  1.55  X  10* 
R 


(AT) 


v«,  a  13,000  FP3 


R^  =  1.266  X  10^ 


1.02 


-0.046 


-1.0 


/  ‘’o.  \ 

1  \ 

(  100  ) 

(A8) 
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0,  Reference  Reynolds  Number  at  Location  of  Maximum  Turbulent  Heating  (R  ) 

’^MAX 

The  significance  of  the  local  reference  Reynolds  number  at  the  location  cor¬ 
responding  to  maximum  turbulent  heating  on  a  hemisphere  and  unswept  cylinder 
is  that  it  is  unique  and  can  be  defined  in  terms  of  free  stream  velocity  and 
altitude.  The  approximate  equations  for  R  for  both  fli^it  and  wind  tunnel 

’^MAX 

are  accurate  to  +  15')^  •  At  first  this  appears  as  a  poor  approximation,  but 
R  is  only  required  to  determine  the  heat  transfer  ratio 

^‘max  h — 

"SP 

location  of  maximum  turbulent  heating,  The  ratio  is  not  sensitive 

^SP 

to  R  .  For  example,  in  the  range  where  turbulent  flow  can  be  expected,  a 

^■max 

10^  change  in  R  will  result  in  a  ^  change 
MAX 

location  of  maximum  turbulent  heating,  is 

independent  of  R  . 

MAX 


in  the  ratio 


for  practical  design  purposes 


1.  Hemisphere 

Approximate  equations  for  the  local  reference  Reynolds  number,  R„  for 

MAX 

a  hemisphere  have  an  accuracy  of  +  15^^.  These  equations  are: 


(«)  Wind  Tunnel 


MAX 


(b)  Flight 


■ 


(A9) 


\L  <  13,000  FPB 


MAX 


’  -  V 

>  100  / 

(AlO) 


‘^MAX 


v«  »  13,000  Fpa 

■ -  1-02  1.232  (v  /p  X-O.OW  -1.0 

,  /  P  \  /V  \  ^  'V W  ,  T  , 

(i?)  (®)" 


(All) 
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2.  Unswept  Infinite  Cylinder 

The  fpproximate  equations  for  the  local  refertace  Reynolds  number, 

R  ,  for  an  unswept  cylinder  have  an  acc^-racy  of  +  15^.  These  equations  are « 
^MAX 

(l)  Wind  T<innol 


R 

*'MAX 


=  a. 3  X  10® 


“  0-659 


(A12) 


<2)  Flight 

v.<  13,000  ?ra 


i 


n 

^MAX 


.  1.70  X  loi'  (pyp^^) 


I.i96(ryp^) 


0.046 


-0.97 

R 


0.97 

(A13) 


\i.i  13,000  FPa 


N.  1  .;0  V  Iv  /P  V  • 


(^)  (^) 


0.771 


0.97 


(Al4) 


i 

► 

3.  Conversion  of  Free- Stream  Conditions  to  Stagnation  Chamber  Conditions  | 

The  free  stream  conditions  P«,  and  are  generally  difficult  to  measure 
in  a  wind  tunnel  whereu.=  stagnation  chamber  conditions  can  be  obtained  with 
relative  ease.  The  free  stream  temperature,  Tg» ,  can  be  related  to  stagnation 
point  enthalpy  through  the  energy  equation  and  becomes  ' 


T 

m 


C 


p 


{A15) 
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where 


c„  =  0.2it  BTU/lb  “R 

p  m 

Equation  (AI5)  applies  if  the  wind  tiinnel  free  stream  flow  has  reached 
chemical  equilibrium.  The  stagnation  point  pressure  can  be  related  to  the 
free- stream  pressure  by  figure  Al, 

Environmental  limitations  of  approximate  expressions  for  wind  tiumel  and 
flight  conditions  are  shovm  in  table  Al. 


TABLE  A1 


PRESSURE 


APFEHOIX  B* 


LAMINAR  AND  TURBULENT  Prt*r  TRANSFER  METHOD 


The  method  used  for  theoretics!  predictions  throughout  this 

report  was  developed  by  Richard  A.  Hanks  in  the  course  of  the  X-20  prograa. 
The  method  is  based  on  the  Integral  form  of  the  boundary  layer  momentum 
equstion.  This  equation  is  transformed  into  an  equivalent  incompressible 
form  that  is  then  solved  to  yield  a  generalleed  equation  that  Includes 
the  usual  boundary  layer  thickness  and  form  factors  as  undetermined  funo- 
tions.  For  laminar  flow  these  functions  were  evaluated  by  equating  the 
generalized  equation  to  exact  numerical  solutions  of  the  differential 
equations  for  self-similar  boundary  layers.  Two  apparently  universal 
boundary  layer  functions  were  evaluated  that  allow  a  general  heat  transfer 
equation  to  be  written  that  agrees  with  essentially  all  of  the  exact 
similarity  solutions  to  within  about  3  percent,  including  the  effects  of 
fluid  property  variations,  finite  streamwise  and  crossflow  pressure 
gradients,  and  streamline  divergence.  The  two  functions  are! 

1.  a  reference  value  of  the  density-viscosity  product  Pj.  Pj. 

that  depends  only  on  the  density-viscosity  products  evaluated 
at  the  wall,  edge,  and  stagnation  enthalpies,  and 

3,  a  boundary  layer  profile  parameter,  F  ,  that  depends  only  on 
the  density  evaluated  at  a  linear  combination  cf  the  wall,  edge, 
n  .id  stagnation  enthalpies. 

No  analytic  derivations  for  these  functions  have  been  found,  and  all  results 
In  this  report  were  calculated  with  the  aid  of  curves  given  in  this  appon* 
dix.  Recently,  however,  simple  expressions  have  been  found  that  agree 
closely  with  the  plotted  curves. 

The  extension  to  turbulent  flows  was  guided  by  the  laminar  results, 
physical  considerations,  and  comparisons  with  experimental  results.  The 
functions  PrMr  T  ■*'®  retained  In  the  turbulent  flow  method,  and 

are  equal  to  the  laminar  va.iues.  However,  the  expressions  used  to  calcu¬ 
late  the  heating  are  of  course  somewhat  different,  and  so  the  effects  of 
Pf  heating  rate  are  also  somewhat  different.  The 

general  form  of  the  basic  momentum  Integral  equation  allows  turbulent 
flow  heat  transfer  data  from  different  sources  to  be  compared  on  a  cona la¬ 
tent  ahd  e.'stematic  basis.  During  the  X-20  program  extensive  comparisons 
to  experimental  data  were  made  in  which  no  data  were  consciously  ignored. 

The  method  described  here  reflects  those  comparisons.  Although  the 
derivation  given  is  restricted  to  the  vicinity  of  a  plane  of  symmetry, 
the  resulting  expressions  have  been  applied  to  infinite  cylinder  dlstrl- 
butlona  with  good  results. 

The  following  discussion  describes  the  method  as  It  la  now  programmed 
for  digital  computers  and  as  It  was  used  In  the  preparation  of  the  present 
report.  Some  modifications  of  the  method  have  been  made  unde?  NASA  contract 


♦This  appendix  la  based  on  Appendix  B  of  reference  16. 
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NA88-11321  which  are  not  Included  la  the  present  c^loulntlonn.  Theap 
■odlfloatlonn  are  primarily  for  flow  conditions  other  than  those  of  the 
data  discussed  in  this  report,  and  do  not  lead  to  appreoiable  numerloal 
differences  here.  Further  information  regardlnc  the  later  modifloatione 
nay  be  obtained  from  reference  29.  The  following  desorlptlon  of  the 
derivation  and  application  of  the  nethod  is  In  six  parts i 

1.  Derivation  of  a  general  form  of  the  boundary  layer  nonentun 
Integral  equation. 

5.  Transformation  of  the  Integral  equation  to  an  equivalent 
Incompressible  form. 

3.  Correlation  of  exact  lanlnar  solutions. 

4,  Evaluation  of  turbulent  boundary  layer  parameters. 

B.  Combined  laminar  and  turbulent  method. 

6.  Summary  of  nethod. 


Derivation  of  the  Itomentum  Integral  Equation 

A  derivation  of  the  boundary  layer  momentum  Integral  equation  la  s 
general  curvilinear  coordinate  system  will  now  be  given.  The  derivation 
is  restrictod  to  the  viol.'«ity  of  a  plane  of  symmetry  as  well  as  by  the 
usual  boundary  layer  assumptions.  A  control  volume  is  defined  so  shoni 
la  the  sketch  below t 


r.,w 


1?2 


The  length  elements  In  the  x  and  s  directions  are  unity.  However,  the 
length  element  for  j  is  determined  by  the  function  g  =  g  (x),  which 
is  considered  arbitrary  (subject  to  the  restriction  that  dg/dx  remains 
finite).  Later  it  will  be  seen  that  in  some  cases  the  most  convenient 
choice  of  g  is  determined  by  the  shape  of  the  body  under  consideration. 

The  height  of  the  control  volume,  h,  is  constant  end  must  be  larger 
than  the  boundary  layer  thickness  but  is  otherwise  arbitrary.  The  surface 
y  B  o  coincides  with  the  straight  streamline |  hence  v  s  o  when  y  >  0. 
However,  v  is  not  necessarily  zero  nor  even  constant  on  any  other  surface 
of  constant  y. 

Maas  conservation.-  The  mass  entering  the  control  volume  through  the 
surface  x  «  x^  la  given  by: 


L  y 


*1 


The  mass  leaving  at  x 
Ixpanding  in  a  Taylor 
yields  t 


B  x.-f  Ax  is  given  by  a  similar  expression, 
series,  and  retaining  only  the  first  order  tern 


Ay 


g  f  p  u  dzj  +  Ay  Ax  ^  g  P  u 


dzl 


(Bl) 


so  that  the  mass  remaining  within  the  control  volume  is: 

JL 


Ay  Ax 


fK 


-h 

g  J  P  u  dz 
0 

mm  mm 


(B2) 


Applying  this  technique  over  all  six  surfaces  of  the  control  volume,  and 
requiring  steady  flow,  yields:  (B3) 

AxAy^j^g  J  pu  dzj  +  AxAy^|j  p  v  dzj  +  gAxAy[p  w]'  =  o 

Since  w(0)  «  0  there  results  in  the  limit  as  Ax  and  Ay  approach  zero: 
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x-moaentua. -  In  a  similar  manner  the  following  expression  for  x-aomen- 
turn  in  obtained; 


XI 

»x 


J  P  <lz  <  J  p  u  V  dz^  +  g  [p  U  wj 


z=h 

z=0 


L  J  z=0 


h  ^  [gp]  +  h  p 


<B5) 


Combining  equation  (B5)  with  the  previous  result  for  mass  conservation, 
and  noting  that  T(k)  s  0: 

^  [k  (  P  ^  [jf  P  ’  I"  -  "e'  <*“1 


=  -  g  -  h  g 


IE 

lx 


(B0) 


Introducing  th«  usual  boundary  layer  thickness  parameters; 

h 


■omentum  thickness 


9 


dieplacement  thickness 


crossflow  momentum  thickness  ratio 

.-h 


'I  ft [i- fell" 

•■■■!•(' -it)' 

«•«*  f  A  '  f 


(B?) 


I  =  E  »  1  f"  -fiJL. 

®  ®  •'O  ^8^0  ' 


dz 


leads  to  the  following  expression  for  Ax  and  Ay  approaching  seroi 


w 


P. 


*  -L  ,  i/lK  ,  I  A  !lo) 

ax  j^Ug  lx  \  e  /  Pe  >X  ^ix  e  Uo  ly  / 


(B8) 
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Definition  of  g(x).-  Since  the  definition  of  g  In  ntlll  arbitrary,  It 
seems  desirable  to  make  a  definition  that  will  simplify  equation  (B8)  If 
possible.  It  might  appear  that  If  g  were  selected  such  that  ■  0 

(l.e. ,  streamline  coordinates)  the  last  term  would  be  made  eero  also.  Such 
Is  not  necessarily  the  case  however,  since  Ve  Is  a  divisor  in  the 
definition  of  E  .  It  will  be  seen  that  In  the  limit  for  small  y,  the 
product 


becomes 


E 


dz 


An  additional  condition,  fv  /  fv  •  0  does  cause  this  term  to  go  to  aero. 
An  examination  of  the  complete  boundary  layer  differential  equation  shows 
that  Iv  /  ly  ■  0  occurs  only  If 

»p/ty  =*  0 


and 

9^  p/iy^  “  0 


These  conditions  are  met  only  If  the  body  is:  1)  two-dimensional,  or  2) 
axlsymmetric  and  at  aero  angle  of  attack.  In  both  cases  the  surfaces  y  • 
constant  follow  streamlines  If  g  «  r  ,  the  local  body  radius.  In  the 
case  of  an  axlsymmetric  body  at  angle  of  attack  the  streamlines  will  not 
follow  surfaces  of  constant  y  however.  The  additional  dlv‘ rgence  is 
denoted  by  f,  defined  by 


1  3  J.  !!e 

f  #x  Uog  »y 


(B9) 


where  y  is  measured  with  respect  to  the  y,  x  coordinate  system.  In  terms 
of  r  an8  t,  equation  (B8)  becomes 


1  it  .  E  /l  lf  \ 
r  lx  6  \f  |x/ 


(BIO) 
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Physically,  the  term  r  may  aI..o  be  thought  of  bh  ntreamllne  divergence  due 
to  body  ahape  while  the  term  f  repreaentn  atreamllne  divergence  due  to 
transverse  pressure  gradients.  The  quantities  r  and  f  are  related  by 

rf  «  A  (B9a) 

where  A  1*  the  total  distance  between  any  two  streamlines.  In  the  absence 
of  transverse  pressure  gradients  the  final  term  in* equation  (BIO)  will  be  aero. 


Transformation  of  the  Momentum  Integral  Squat ion 

In  order  to  ol)taln  n  more  useful  form  of  the  momentum  equation  (BIO)  a 
modified  Stewartson  transformation  suggested  l)y  Mngcr  (ref.  30)  is  adopted  in 

-  r*  PrMr 
*  =  I  F  dx 

•^0 

^ "  y  z  (Bii) 

Z  =  F  f  ^  dz 
•'o 

U  =»  u/k 


V  =«  v/F 


where  X,  Y,  Z,  U  and  V  are  the  transformed  coordinates  and  velocities.  The 
stagnation  values  of  density  and  viscosity^  Po  and  are  required  to  be 
constant,  and  F  is  an  unspecified  function  of  x  only.  With  these  definitions, 
the  momentum  thickness,  skin  friction  at  the  wall,  and  heat  transfer  rate  in 
the  transformed  coordinate  system  are  respectively: 


-  _  ^w/^o 
F  VPi-Pr/ 

1  /Po  Po\ 

^  *  „2  1 0  “  ' 

^  F  \f'r  '-r/ 

Thw  traasformed  iioaontum  Intogral  equation  becomes 

_  p  A  /l  ^  a  f-A_  JL  i^.  ,rv  ^  A 

PeUe*  '  ^  [u^  F  »X  <  e  '  p„  5X 

^9  ml 

^  1  dr  -=  1  (A  1 
r  dh  ®  f  ilX  J 


(Blla) 


J  dr  -= 
+  -  — -  ♦  E 
r  dh 


(B12) 
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in  which  A  -  2  +  —  (A  and  B  are  unchanged  by  the  transformation). 

Solution  of  the  transformed  Integral  equation.-  It  la  assumed  that 
the  local  friction  coefficient  In  the  transformed  plane  Is  given  by 


(B13) 


Equation  (B13)  is  substituted  Into  equation  (B12)  yielding  an  equation 
of  the  form: 


X 

^  +  P(X)  0  =  e  Q(X)  (B13a) 

which  becomes  a  linear  first  order  equation  with  the  change  of  variables 


0 


m+1 


(B13b) 


The  solution  obtained  las 


™  in 
PoPo^e 


la  the  untranafomed  physical  plane,  (B14)  la; 


Jffll 

Ue 


1-m 

m 

^  m  P  o  Pr  P  r  **8 


A  4m-) 


(rfS)' 


l--m 

m 


PrPr  «e 


m 


m4l 

m 


Im4l 


dx 


(B16) 


Neither  p  nor  F  appear  In  this  aquation,  and  their  definitions  are 
therefore  Immaterial. 
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ReaeabarliiK  that  and  #1q  are  aaauned  to  bo  independent  of  x  , 

and  defining 


(s)' 


(B16) 


0  =  Prl^r%^ 


(BIB)  can  then  be  reduced  to: 


1-m  m 

^  1+iM  _m+l 

CyXo  yr*‘r“el 


0 


/  - 


(B17) 


Noting  that  the  quantity  within  the  bracketa  in  the  denMlnator  on  the 
right  hand  aide  haa  the  unite  of  length,  we  can  de^:^np 

~  /’'*  o(f® 


oUE  A-1 
G\f 


where  the  aubacript  1  Indicatea  evaluation  at  the  current  point  of 
Intereat  xi  •  All  effects  due  to  flow  three-dimenalonality ,  atreamwlse  and 
tranaverae  preaaure  gradlenta,  and  upatrean  blatory  are  now  Included  in 
Beq  ,  if  independent  of  these  phenonena  (the 

validity  ox  this  assunptlon  will  be  demonstrated  subsequently).  Thus 
can  be  conaldered  to  be  "the  equivalent  flat  plate  distance''  for  skin 
friction  at  the  point  of  Interest  z^. 


With  (B18),  (B17)  can  be  rewritten: 


C  u  p  U  u  S  m+1 
X  ’^r  e  eg  i 


158 


which  is  IdcntJcfll  In 
flat  plate  flow, 

T 

w 


One  approach  to  the  solution  of  a  boundary  layer  Integral  equation,  such 
ns  developed  In  the  preceding  analysis.  Involves  the  assumption  and  Inte¬ 
gration  of  boundary  layer  profiles  to  obtain  the  required  boundary  layer 
thickness  parameters  (see,  for  Instance,  Reckwlth  and  Gallagher  (ret.  31^. 
In  that  approach  it  Is  necessary  to  derive  the  energy  equation  corres¬ 
ponding  to  (1U5).  An  alternative  method  Is  used  here,  wherein  a  general 
form  of  Heynolds  analogy  Is  assumed 


form  to  the  corresponding  expression  for  low  speed 

m 


Re,e 


in+1 


(n20) 


H  = 


_ SL _ 

*ttW  "  W 


(021) 


It  la  of  course  weJ 1  known  that  the  Reynolds  analogy  factor  J  has  the  value 
unity  for  constant  property,  unity  Prandtl  number,  flat  plate  flow.  It 
will  be  subsequently  demonstrated  that,  In  the  presence  of  more  realistic 
gas  properties,^  is  for  laminar  flat  plate  flows  still  a  function  only  of 
the  Prandtl  number  and  (In  dissociated  flow)  the  Lewis  number.  For  conve¬ 
nience  denoting  these  flat  plate  flow  functional  relationships  by  Fp^ 
and  ^  respectively,  (B21)  Is  rewritten 


H  = 


/ 

Fpp  8  u^ 


(1122^ 


where  the  factor  S  Incorporates  all  effects  of  flow  three-dimensionality, 
atreamwlse  and  transverse  pressure  gradients,  and  upstream  history  on 
Reynolds  analogy. 


Combining  (D17),  (D18),  and  (B22), 


H 


(D23) 


Since  Snq  Includes  all  effects  of  flow  history  and  pressure  gradient  on 
skin  friction,  and  3  performs  a  similar  function  on  Reynolds  analogy, 
the  definition  of  an  equivalent  distance  for  heat  transfer  suggests  Itself. 
Accordingly,  we  define 

m+1 
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^eq  “  S 


1  ^ 
A  r 


o(fE 


im-1 

m 


m+1 


= —  dx 


o(®  »/-*>)  "■ 

L  '  '  J 


<B24) 


leading  to  a  general  expreaslon  for  heat  transfer  of  the  fora 

|m/(m+l) 


H  = 


C  jf  U 

X  ^ 

^Pr  ^cq 


Pr  Pr  "e  ^ 


eg. 


Mo' 


(B25) 


Means  of  evaluating  the  various  paraaeters  appearing  In  (B26)are  presented 
In  the  following  sections. 


1 

j 

! 
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Evaluation  ol  Laminar  Boundary  Layer  Parameters 


Exact  solutlone  of  fhe  .nlmllarJty  form  of  the  laminar  boundary  layer 
equatlona  were  uned  to  evaluate  the  parameters  appearing  In  equation  (B25). 
This  was  done  In  sn.  orderly  manner,  beginning  with  two-dlmcHHlonal  constant- 
property  constant-pressure  flow,  and  progressing  to  the  most  complex  condi¬ 
tions  for  which  exact  solutions  are  available.  The  evaluations  determined 
from  the  simpler  cases  were  retained  or  amplified  in  analyzing  the  more 
complex  cases.  Thus,  the  constant  C  for  laminar  flow  Is  always  taken  to 
be  0.33206,  the  value  given  by  Howarth  In  reference  32  for  Incompressible 
flat  plate  flow.  The  effects  of  pressure  gradients,  wall  cooling,  etc. , 
are  accounted  for  In  other  terms  of  equation  (B25). 

In  some  cases  alternative  definitions  were  possible.  For  example, 
the  authors  of  references  35  and  38  Incorporated  (In  effect)  pressure 
gradients  Into  the  term  pf /ir  appearing  In  equation  (DIO),  while  In 
the  present  formulation  such  effects  appear  In  the  equivalent  distance, 

Xeq.  The  latter  definition  Is  to  be  preferred  as  the  former  cannot  be 
made  consistent  with  the  results  of  reference  33,  which  presents  solutions 
for  various  pressure  gradients,  but  with  pp  held  constant.  The  defini¬ 
tions  used  here  were  adopted  only  after  an  examination  of  several  possible 
alternatives.  The  criteria  for  selection  were  consistency  l)etween  the 
results  of  the  various  special  cases,  consistency  with  physical  considera¬ 
tions,  accuracy,  simplicity,  and  freedom  from  Interdependencies. 

Ooneral  considerations.-  As  a  matter  of  physlcrl  consistency,  It  Is 
required  that  If  the  fluid  properties  p  and  p  are  constant  through  the 
boundary  layer,  the  reference  values  of  the  fluid  properties  l)c  equnl  to 
those  constant  values.  This  principle  la  extended  to  constant  products 
as  well,  l.e.  ,  It  Is  required  that  when  In  a  git'en  numerical  calculation, 
e.g. ,  references  33  and  34,  the  product  of  density  and  viscosity  Is  held 
constant  at  sor.'>  base  value  (usually  the  wall)  the  reference  density 
viscosity  product  P|.  must  also  be  equal  to  that  base  value.  The 
functions  Fpr  and  ^  are  equnl  to  1.0  when  a  and  Lo  are  equnl  to  i.O,  and 
^  e  1,0  for  ideal  gases.  Also,  In  flat  plate  flow  the  equivalent  distance 
Is  equal  to  the  physical  distance  from  the  lending  edge. 


Two-Dimensional  Flat  Plate  Flow 

The  special  case  of  two-dimensional  flat  plntc  flow  Is  examined  first 
since  the  effect  of  fluid  property  variations  within  the  boundary  layer 
can  bo  examined  without  the  additional  complexity  of  strcnmwlse  variations. 
For  the  cn.se  of  constant  fluid  properties  the  snlntlons  of  Howarth  show 
that  m  =  1  and  -  .332,  so  that  equation  <D25)  becomes 


l6l 


(D26) 


(B27) 


where  (B27)  follows  from  the  prlticlplee  stated  under  "Oeneral  Considerations**. 
Tor  this  special  case  the  only  undetermined  quantity  Is  th*  Reynolds 
analogy  factors  ^  and  Fp,  .  Kote  that  the  reference  stagnation  i^lsooslty, 

,  no  longer  a:;>peara. 

Reynolds  analogy  factors,-  The  Prandtl  number  affect  on  Reynolds 
analogy  In  flat  plate  flow,  usually  given  an  Fpr  >=  a  for  constant  O 

la  slightly  better  represented  by  a  as  may  be  seen  In  figure  (Bl). 

Following  the  practice  of  reference  35,  for  example,  the  Prandtl  number 
effect  is  correlated  In  teres  of  (T  ,  the  partial  Prandtl  number  for 
translation,  rotation,  and  vibration. 


For  variable  Prandtl  number  there  Is  an  uncertainty  as  to  which  value 
should  be  used  In  correlating  Its  effect.  All  aolutlous  In  the  literature 
for  which  the  Prandtl  number  is  variable  also  involve  variable  pil, 
so  that  p  JU  1*  not  necossarlly  equal  to  •  For  such  caaea  It 

waa  found  that  the  Prandtl  number  should  be  evaluated  at  the  enthalpy  an<^ 
pressure  corresponding  to  PfUr  •  ‘^1*  value  of  the  Prandtl  number  is 

hereafter  denoted  as  Oj,  ‘me  adequacy  of  this  avalustlon  is  demonstrated 

by  the  agreemeut  of  the  three  aeta  of  calculationa  presented  in  figure  (Bl), 
which  also  serves  to  establish  the  lack  of  dspendanoy  of  Fpr  on  anything 
other  than  Of. 

With  the  Prandtl  number  effect  oorrelnted  In  terms  of  the  partial 
Prandtl  number  the  effect  of  energy  transport  by  diffusion  must  be  treated 
separately.  This  effect  waa  first  calculated  by  the  outhora  of  raferenoa  35, 
wherein  the  expression 


^  Le 

/ - 

^  Le  s  1 


1  +  (  Le -  1) 


(B28) 


was  found  to  agree  well  with  exact  aolutlons  for  Le  •1.4,  in  staicnation 
point  flow.  In  high  liaoh  nuaber  flows,  however,  equation  (B28)  may  predict 
a  significant  diffusion  effect  under  conditions  for  which  no  dlssoolstloa 
actually  exlsta,  since  the  temperatures  within  the  boundary  layer  are  always 
well  below  the  stagnation  value.  To  avoid  this  ir.«ounMistanoy,  equation  (B28) 
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was  modified  to  operate  on  the  local  static  enthalpy,  rather  than  the 
stagnation  value.  The  modified  expression. 


/-=  1  +  (Le  -52  _  (1329) 

*e 

of  course  reduces  to  (B28)  for  stagnation  point  flow.  Equation  (B29) 
was  used  for  all  calculations  in  the  present  report,  although  later 
publications,  reference  37,  for  example,  indicate  that  equation  (B29) 
overestimates  the  heating  rate  hy  5  to  10“  in  some  cases. 

Reference  density-viscosity  product.-  The  reference  density-viscosity 
product  was  first  evaluated  for  zero  Mach  number  with  various  degrees  of 
wall  cooling  using  the  solutions  of  references  36,  38,  and  39,  and  some 
unpublished  solutions  by  Malvorson  and  Cassmeyer  of  The  Boeing  Company, 
as  shown  ;tn  f  Igure  B2. 

For  edge  Mach  numbers  greater  than  zero  it  was  found  that  the  reference 
density-viscosity  product  Mr  represented  as  a  function  only 

of  Pq  Pp  ,  p  and  pg,  |igi  ,  where  the  latter  is  the  density-viscos¬ 

ity  product  evaluated  at  stagnation  enthalpy  but  the  local  pressure.  Using 
the  solutions  of  references  36  and  38  an  effective  edge  value  of  p  p  was 
determined  that  allots  the  use  of  figure  B2  for  Mach  numbers  other  than 
zero.  The  effective  p  p  product  (pp  Mc'gff  found  to  be  a  function 

of  Pat  Mqi  >  Pe  Mo  >  only.  The  curve  that  defines  this  relation  is 

given  In  figure  B3  (a).  All  of  the  solutions  discussed  so  far  are  well 
represented  by  the  faired  curve  of  figure  B2  when  plotted  against  Pg  Mg  eff 
PS  may  be  seen  in  figure  B3  (b). 

Subsequent  Investigations  described  below  have  shown  that  Pj- pj- 
is  independent  of  pressure  gradients.  The  values  of  Pr  Mr  obtained  from 
figures  B2  and  B3  were  used  for  all  calculations  appearing  in  this  report.* 


Pressure  Gradient  Effects  -  Similar  Flows 

Referring  to  equation  (B25),  and  recalling  the  earlier  comment  that 
Cjj,  m  ,'Fpr  and  ^  are  by  definition  taken  as  the  flat  plate  values  it  is 
seen  that  all  pressure  gradient  effects  are  reflected  in  Pj.  p^  and  Xgq 
These  effects  can  be  evaluated  for  similar  flryws  from  the  solution  publlsned 
(for  example)  in  references  33  and  38,  which  consider  streamwlse  pressure 
gradients,  and  references  34  and  37  which  consider  cross-flow  pressure 
gradients. 

♦Recently  some  simple  expressions  have  been  found  which  approximate 
curves  of  figures  B2  and  B4  closely.  The  expressions  are  given  In  the 
final  section  of  this  appendix. 
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Beglnnlns  with  the  elnpleet  poeelble  caee,  two-dlwenaion«l  flow  of 
an  ideal  gas  with  Pfandtl  number  of  one,  and  the  viscosity  proportional  to 
temperature,  the  equivalent  distance  effects  can  be  Isolated.  Since 
la  alwaya  equal  to  PePe*  PrMr  Is  also  equal  to  Pe  Me  *  (Note 
that  p  |l  la  not  necessarily  constant  through  the  flow  field,  but  varies 
with  the  local  boiuidary  layer  edge  pressure. )  With  these  values  Incor¬ 
porated,  the  equivalent  distance  expression  (B24),  Is  reduced  to 


(B30) 


In  equation  (B30  the  term  (P  u,,)  reflects  thg. effects  of  upstream  variations 
In  Pe  Me  ^  •  '"'Nile  the  terns  8*  and  Ug  ’  account  for  local  pressure 

gradient  effects  on  the  boundary  layer  profiles. 

Equation  (D30)  can  be  evaluated  If  S  and  A  are  known,  and  although 
Ijtborlous,  they  could  be  determined  from  the  numerical  solutions.  For¬ 
tunately,  specific  evaluation  of  these  parameters  has  proven  to  be  unnec¬ 
essary,  since  a  convenient  simple  correlation  has  been  found  for  the 
oonblned  effect  of  8  and  which  may  be  written  as 


a  f  ‘  (P  u,l  ■!«  1  P 

“  0  [<PUe)Ue2(A-»]  <P 


(B31) 


where  fi  is  the  dimensionless  pressure  gradient  parameter  similar  to  that 
defined  by  the  authors  of  reference  33.*  The  profile  parameter  F  Is  a 
single  valued  function  of  •  neau  boundary  layer  density,  p^^  ,  calculated 

by 


P  2T 

^  =  lie (332) 

Pe  (Tw  T«,) 


•ubaequent  Investigations  of  exact  solutions  for  nonunity  Prandtl  number  md 
nonlinear  viscosity  laws  have  shown  that  expressions  of  the  form  of  (B31) 
are  valid  for  these  more  complex  conditions  as  well,  either  for  two-dlmon- 
slonnl  flows  with  streamwlse  pressure  gradients,  or  for  yawed  cylinder  flow. 
The  expressions  finally  developed  arej  _ 

*Tbo  definition  of  p  Is  given  In  equation  {DiiTT 


(Pr 

Pr)|9  "  (PrPr)/3=o 

and  a 

generalization  of 

(B31) 

1 

I  9  t 

clx 

’‘CHI.I 

'  "  Jl 

J  r- 

«  1 

j 

XI 

where 

Jj^  Is  given  by 

and 

r„ 

Pc 

(034) 


(U35) 


(B3G) 


where  the  subscripts  "s"  and  "c"  are  Introduced  to  distinguish  between 
streanwlse  and  crossflow  pressure  gradients;  It  should  be  noted  tlyit 
is  concerned  only  with  streanwlse  pressure  gradient  effects  and  E|^ 
with  crossflow  effects;  also  note  that  J|^  s  1.0  for  »  0  and 
-1.0  for  p  c  “ 

The  function  P  Is  given  for  either  streanwlse  or  crossflow  pressure 
gradients  bv  the  curve  of  figure  B4  as  a  single-valued  function  of  a 
paraneter  X  defined  by: 


(Z  T) 


E  -  - 
«  P 


m.B 


(ZT)^ 


(B37) 


for  streanwlse  pressure  gradients,  and  as 


.SL  _  (^  '^)in.o 


(B38) 


for  crossflow  pressure  gradients.  The  subscript  "n"  denotes  evaluation 
at  a  Mean  boundary  layer  enthalpy,  defined  by: 


in  ,B 


(‘S  *w) 


(B39) 
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(D40) 


aiui 


1 


m,o 


K  SL  +  ‘w)  +  -zoa  (ig  -  i^  sO 


f  Ps  ^  8 


Tne  aecond  ei^uallty  In  equatlon.4  (B37)  and  (B3B)  followa  from  the  condition 
of  constant  pressure  scrons  the  boundary  layer  (all  evaluations  are  made 
at  the  lo-ial  pressure).  Aifaln,  (B37)  throuRh  (B40)  are  the  generalisations 
of  equation  (B31)  and  (B32)  . 

While  figure  B4  and  equations  (H34-40)  were  developed  solely  on  the 
basis  of  providing  the  best  fit  to  the  available  data  within  the  framework 
of  the  form  of  equation  (n31)  the  obvious  similarity  of  (n39-40)to  the 
various  reference  enthalpies  appearing  in  the  literature  provides  some 
analytical  Justification  for  these  correlations. 

The  streaihwlse  pressure  gradient  parameter  0  ^  Is  herein  defined  asi 


d  Qn  11^) 
(I  (In  x) 


{B41) 


With  a  minor  modification  to  the  definition  of  Xeq^L  (to  be  discussed 
in  the  next  section),  It  Is  easily  shown  that  this  definition  of  ftg 
is  Identical  to  the  corresponding  parameter  of  reference  37.  In  the  present 
report  f,  la  evaluated  only  for  yawed  cylinder  flov,  in  which  case  p  o 
1.0,  also  consistent  with  the  results  of  reference  37.  Except  as  otherwise 
noted  In  the  text,  a  value  of  unity  was  used  for  all  leading  edge  theory 
calculations,  while  a  value  of  zero  (corresponding  to  E  =  1.0)  was  assumed 
for  lower  surface  theory  calculations. 

The  exponent  O  In  equations  (B35)  and  (1136)  Is  given  by 


O  “ 


(B42) 


which  is  also  based  on  fits  to  the  solutions  of  references  33,  34,  35,  and 
37.  The  accuracy  obtained  through  the  use  of  equations  (B33)  through 
(B43)  Is  illustrated  In  figure  B5,  wherein  solutions  from  references  33 
and  37  are  presented  In  terms  of  i  and  £  .  As  may  l>o  seen  by  comparing 
the  spread  of  the  Individual  numerical  solutions  with  the  Indicated  error 
band,  the  present  method  provides  excellent  agroement  with  all  solutions. 
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Application  to  Wonslailar  Flow 


lha  wquatloM  prwaantad  In  this  awotloa  have  all  bean  davalopad  fro« 
alHllaritjr  aolatloaa  to  tha  boundary  layar  aquations,  and  ao  ara  strictly 
applieabla  only  la  thosa  sltuationa  for  which  slnilarlty  applies.  Howarar, 
baaad  on  disoussioaa  givan  in  rafaraaeaa  37  and  38,  it  is  to  ba  aapactad 
that  tha  aaao  oorralations  could  provida  good  aatlaatas  for  nonslnllar 
flow  oondltlons  as  wall,  provided  their  straanwlsa  variations  ara  taken 
into  account.  This  has  bean  dona  by  incorporating  tha  factor 
appearing  in  aquation  (B38)  into  tha  integrand.  Tha  expression  for 
aqulvalant  dlatanoa  than  baconast* 


X 


oq,L 


(B43) 


It  is  easily  saaa  that  for  aiailar  flows,  wherein  J^,  is  constant,  (B43) 
radvioas  to  (B34) .  Ivslttatad  for  two-dinanslonal  flows  (f  B  1.0)  the 
use  of  CA43)  iB  (Bll)  results  in  a  definition  of  fig  Identical  to  that  of 
rafaranoa  37,  while  tha  use  of  (B34)  la  (B41)  provides  a  value  of  fi  g 
that  oorresponds  to  tha  "looal  slnilarlty"  approach  of  rafaranoa  33, 
wherein  tha  upatrasn  history  of  profile  affects  ara  naglactad  (that  is, 
tha  boundary  layar  profiles  ara  asaunad  to  adjust  Instantaneously  to  tha 
looal  pressure  gradient). 


Iquatlon  (B43)  assunas  tha  crossflow  paraaeter  E]^  to  ba  Independent 
of  straanwlsa  pressure  gradient  ef foots;  however,  the  presence  of  the  fig 
orossflow  tarns  in  the  resulting  definition  oi  fi  g  provida  a  coupling  between 
tha  transverse  and  straanwlsa  pressure  gradient  affects,  as  night' be 
axpaotad.  Tha  overall  affect  predicted  Is  in  qualitative  agreenent  with 
tha  results  of  rafaranoa  40;  imfortuaately,  the  difficulty  of  relating 
tha  correlating  paranatars  of  that  reference  to  the  present  systen  has 
so  far  prevented  quantitative  oonpcrlsons. 


Finally,  for  tha  generel  case  of  curved  straanllnas  (l.e.  ,  away  fron 
a  line  of  synnetry) ,  it  Is  assuaed  that  tha  foregoing  analysis  and  oorrala- 
tlona:  ara  valid  If  the  distance  paraaeter  x  is  taken  to  be  neasurad  along 
tha  atreanllne.  As  previously  noted,  tha  definition  of  fig  for  the  general 
ossa  Is  obscure,  and  all  oaloulations  herein  c'her  than  leading  edge  valves, 
have  been  aade  on  the  basis  of  fic-0 1),  which  oorresponds  to  the 
"saro  crossflow"  nethod  of  rafaranoa  41  as  far  as  three-dlaenslonal  effoots 
are  oonoarned.  All  present  oaloulations  do  retain  the  effect  of  strsanwlse 
j^rassure  gradients,  however. 


1iirbul*tit  Flow 


IlMr*  miw  oxaet  ealoulatioaa  of  turbuloot  boundorjr  loyor  flow,  oo 
tkat  o  dovolopMat  of  tho  tjpo  Juot  glwoa  io  not  poosi'blo.  Howovor,  oquo- 
tloo  (B2S)  otlll  Miwo*  to  Idontify  tho  inportont  poraaotora,  and  provldoa 
a  baaia  fo/  oonalatontly  eomiarli^  oxporiwontal  roanlta.  In  partlonlar, 
tho  tranafonatlon  Introduced  with  oqoationa  (BIO)  and  (Bll)  allowa  an 
onpirleal  inooaq^roaaiblo  akin  friotion  law  to  bo  uaad  la  place  of  tho 
(nenoxlatont)  oxaot  flow  aolotlona.  Aa  atated  in  tho  aoction  "Tranafematien 
of  tho  noaontuai  integral  oquaticm,"  tho  tranafonuktion  uaod  ia  baaed  on 
tho  work  of  Hager,  roforonoo  30.  A  diaouaoion  of  the  roaaonii^  behind 
tho  tranafomation  itaolf  ia  given  in  hla  paper.  The  preaent  netbad  departa 
fron  the  auggeationa  of  Hager,  however,  in  the  evaluation  of  the  varloua 
boundary  layer  paraaetcra.  ?1r.  Hanka  waa  guided  in  the  evaluationi  by  the 
valuea  of  the  correaponding  lanlnar  paraawtera,  an  iqtproaoh  whiob  uaa 
auggeated  by  the  very  auooeaaful  reaulta  of  the  firat  auoh  attenpt,  wherein 
tho  laminar  valuea  of  p.  were  -aaed  without  nodifloation  for  turbulent 
flow.  The  reaulting  prediotlona  were  in  excellent  agreement  with  leoently 
obtained  free  flight  data,  aoae  of  which  (notably,  that  from  the  X-IS 
program)  ?:<jre  not  ia  agreemaut  with  any  of  the  well  known  methoda. 

Hr.  Hanka  waa  alao  guided  by  the  requirements  of  a  design  project, 

and  ao  waa  ooaatralned  to  make  ocmaervative  approxlnatlona  where  approsl- 
xationa  were  required.  Thaa,  the  effeota  of  atreamwiae  and  tranaveiirae 
preaaure  g&'adienta  on  tho  turbulent  bowjidary  layer  profllea  wore  iaoluded  in 
the  oaloulatione ,  even  though  it  waa  known  that  the  effeota  wore  amll  and 
could  only  be  crudely  eatlmatsd.  The  available  evidence  indioated  that 
auoh  effeota  would  Incraaae  heat  tranafer,  ao  that  neglecting  them  would 
be  unoonaervative. 

InnoBpreaalble  flow  friotion  law.-  In  order  to  determine  and  a  a 
formula  for  akin  friotion  in  inoompreaaible  turbulent  flow  la  required. 

After  a  aurvey  of  propoaed  inoompreaaible  friction  formulae  a  minor  dcdlfloa" 
tion  of  the  Bohultn-Orunow  (ref.  42)  equation  waa  aeleotedi 


C 


( 


_ tarn _ 

[log30(®«  +  3000)]*’**®^ 


(B44) 


The  modifioatlon  that  waa  made  la  the  addition  of  the  oonatant  (3000)  to 
the  Beynolds  number  that  appears  in  the  denominator.  This  modifloatl<.>a  waa 
made  because  Mr.  Hanka  felt  that  the  high  valuea  of  predicted  by  \;he 
unmodified  equation  at  Reynolda  number  belcw  10^  were  not  reallatio  lii 
view  of  the  well  supported  predlstion  of  stability  theory  that  the  Income 
prepalble  laminar  boundary  layer  la  stable  at  Reynolda  leas  than  about 
60,000. 
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Thg  modified  and  unmodified  expressions  sre  shown  in  figure  B0, 
together  with  some  other  proposed  methods.  As  shown,  there  is  little 
difference  between  the  various  methods,  except  that  the  niaslus  equation 
falls  low  at  high  values  of  Reynolds  number.  Equation  (B44)was  originally 
selected  because  of  its  slight  conservatism,  although  any  other  expression 
could  have  been  used. 

The  form  of  equation  (044)  does  not  lead  itself  to  calculations  in 
the  framework  of  equation  (Q23)  due  to  the  variation  of  with  Reynolds 
number.  However,  comparisons  have  been  made  that  show  that  me4  Is  an 
adequate  approximation  for  evaluating'  geometric  effects.  For  example,  if  m 
is  evaluated  at  particular  values  of  Reynolds  number  using  equation  (B44) 
the  following  comparisons  are  obtained: 


Re 

e 

(1  (In  Cf) 

h 

cone 

^cylinder 

j  f. - i — 

(In  Re^) 

^flat  plate 

^cylinder*  m=4 

10“ 

3.4S 

1.20 

.90 

4.2x10* 

4.0 

1.17 

1.00 

10* 

4.34 

1.16 

1.00 

10* 

6.12 

1.11 

1.01 

Thus  the  effect  of  variations  In  »  la  seen  to  be  small.  Accordingly, 

■a4,0  has  been  selected  for  the  calculation  of  geometric  effects  (e.g., 
^cone'^^flat  plate^  herein.  However,  for  actual  calculations  of  C| 

equation  (B44)  was  used  as  there  can  be  considerable  error  In  a  •  constant 
approxlmatlonu  for  absolute  values  of  Cf.  An  example  of  such  a  friction 
law  Is  the  Blaslus  method,  for  which  a  >  4.  As  shown  In  figure  B6  the 
Blaslua  equation  falls  well  below  the  other  methods  at  high  Reynolds 
numbers. 

Denslty-Vlscoalty  producti  Reynolds  Analogy  Factor.-  As  already  noted 
the  reference  density-viscosity  product  for  turbulent  flow  Is  taken  to  be 
the  laminar  value.  This  basic  Identity  was  suggested  by  the  fact  that 
appears  only  in  connection  with  the  laminar  shear  terns  of  the  turbulent 
boundary  layer  equations.*  It  Is  also  assumed  the  effects  of  Prandtl 
number  and  Lewis  number  on  turbulent  flat  plate  heat  transfer  are  also 
Identical  to  the  laminar  values.  The  use  of  the  laminar  flow  Prandtl 
number  effect  la  common  practice.  The  Use  of  the  laminar  flow  l^wit  number 
effect  la  thought  to  be  a  conservative  upper  limit, 


eSee,  for  example,  equation  13  In  reference  30. 
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Jto^0r»nc«  ■tagnatlon  vlacoslty.-  Since  the  reference  etecnatlon  vle>  ! 

coelty  |io  doee  not  appear  in  the  laminar  equations,  no  information  regar- 
dlng  its  Evaluation  can  be  obtained  by  examining  the  laminar  solutions. 

The  reference  stagnation  viscosity  is  assumed  by  Mager  (ref.  30)  to  be 
the  viscosity  evaluated  at  stagnation  conditions.  For  real  gases  with 
the  viscosity  dependent  on  the  pressure  It  seems  more  realistic  to  consider 
the  local  flow  composition  rather  than  the  composition  corresponding  to 
stagnation  conditions.  Accordingly  |Iq  is  calculated  with  the  Sutherland 
law  and  |1|.  using  the  value  of  specific  heat  corresponding  to  Pj,  /II,. 

The  result  is : 


vhlch  expresses  the  well  known  tact  that  pressure  gradient  effects  on 
turbulent  heat  transfer  are  much  smaller  than  those  In  laminar  flow.  After 
an  examination  of  available  experimental  data  the  value 


Jf!!  *•  1 

-i - -  -  .48  <B48) 

Jl  -  1 


was  selected.  By  analt^y  to  equation  (B35)  it  is  assumed  that 

Jt  =  [l  +  .48  Fg  (B49) 


where  Fg  and  are  the  previously  described  laminar  values.  The  small 

exponent  ot  is  assumed  equal  to  the  laminar  value  (eq.  (ii4S)  )  although 
its  ultimate  effect  on  the  predicted  heating  rate  la  only  about  0.3%  for 

Pb  - 


similarly,  the  behavior  of  B  In  turbulent  flow  can  be  described  only 
qualitatively,  and  most  published  analyses  nealect  Its  effect.  However, 
its  effect  is  to  Increase  heating  rates,  and  was  therefore  Included  in 
the  present  method.  As  In  the  streamwlse  pressure  gradient  caa-^c  the 
values  actually  used  were  based  on  modifications  of  the  oorrespundlng  . 
laminar  correlations.  Unlike  the  streamwlse  parameter  however,  Rf, 

Is  strongly  Influenced  by  Mach  number,  as  evidenced  in  equations  (B38) 
and  (B40) ,  so  that  a  dual  modlfioatlon  la  indicated, 


Considering  first  the  case  for  aero  Mach  number  flow,  it  is  seen 
from  the  definitions  of  (B7)  that  the  upper  limit  on  B  is  6*/d  unless 
the  transverse  velocity  component  v  within  the  boundary  layer  exceeds 
the  esternal  value.  Laminar  solutions  (ref.  34)  rhow  that  these  velocity 
overshoots  do  not  occur  for  cold  wall  aero  Mach  number  flow,  hence  a 
oorrection  factor  of  the  following  form  is  suggested t 


Et,  0  -  I 

®L,  0  ■  ' 


(B50) 


The  precise  value  of  the  constant  C  in  this  expression  cannot  be  calcu¬ 
lated,  of  course.  Tor  the  previously  mentioned  upper  limit  case,  C  of 
course  is  equal  to  -1|  however,  in  the  interests  of  conservatism,  a  value 
of  C  <■  el  was  selected  to  represent  an  upper  limit.  Consistent  with 
equations  (B48)  and  (848),  them 
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I 


^T.O  ‘  * 

In  equations  (BSDEl^qIs  Just  Ei  evaluated  for  Mach  nunber  equal  to  zero. 
For  Mach  number  zero  equation  (B40)  reduces  to 


*m,o,0  2  ^ 


since 


<^S  ■  ‘e,  8L)  -  ® 


at  zero  Mach  number. 

The  effect  of  Mach  number  on  Eipwas  determined  from  observation  of 
emplrloal  trends  In  turbulent  yawed  cylinder  stagnation  line  heat  transfer 

data,  as 


(B5a> 


An  equivalent  form  that  Is  more  convenient  for  computer  applications  has 
been  used  for  all  calculations  In  this  report: 


.  T, 

*eq,  T,  0 

whsrsf  in  general 

’‘eq.  T.  0-j, 


eq,JU. 


eq,  L,  0 


m-p 


Cl  Jl 


’‘I 


eq,  L,  0  J 


eq,  L  J 


L  -  _ 

'>■  -I  [o,,.tL'“  "],j 


(B53> 
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It  10  •Miljr  dMKmstrsted  that  for  laflBlto  yawod  eyllmior  ataciMtloB 
ItM  flov,  (B59)  io  oxoetly  cqulvaloat  to  (B52).  For  other  typee  of  flow 
(BBS)  mad  (BBS)  ore  not  exootly  equlrolent.  However,  elaoe  (BBS)  lo 
booed  oa  yawed  oyllader  data,  (BBS)  and  (BBS)  are  equally  valid  aaauaptlono, 
and  (BBS)  haa  beea  found  to  be  wore  eonvenleat.  la  any  caoe,  the  final 
effect  ea  the  predicted  heatiag  rate  io  oaall. 


Conbined  Lanlnar  -  Turbulent  Method 

A  coapariaoB  of  the  equivalent  distance  expreosioaa  for  laainar  and 
turbulent  boundary  layers  shows  that  ia  general  the  two  values  are  not 
equal,  so  that  the  Reynolds  nuaber  based  on  the  equivalent  distance  will 
depend  on  the  boundary  layer  state.  This  inoonsisteacy  caa  be  avoided  by 
eaploying  the  following  deflnltica  of  a  reference  Reynolds  nuaber. 


R 


r 


a 

F  2 


*eq,L  ^ 


(B64) 


whers 


(B6B) 


When  terms  R  and  F  are  used  inequation  (B85)  there  results 
r  X 


(B66) 


vhsre  Cf,r  friction  coefficient  evaluated  for  the  reference 

Reynolds  nuaber,  Tbs  foraulas  used  in  the  present  report  are: 


^  .664 

C(,r,L  “  1/2 

R 


(B67) 


%r,T 


.370 


2.584 


(B68) 
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Vote  that  in  equation  (B56)  only  Cf  ^  depends  on  the  boundary  layer  state. 
As  a  result,  a  staple  relationship  Mtween  laolnar  and  turbulent  flow  can 
be  obtained  that  is  dependent  upon  the  reference  Beynolds  naaber  only. 

This  reiatinnshlp,  obtained  by  coabinlnK  equations  (B56),  (B57)  and 
(B58),  follows 


^  .186  _ 

“t  .332  [l<igj^  (8,  ♦  3000)]^‘°“ 


(B59) 


It  is  easily  shown  that  the  laoinar  fom  of  equation  (B56)  reduces 
identically  to  any  of  the  special  cases  prevlouiily  given.  For  exaople, 
considering  only  the  various  equivalent  distance  teias,  and  enploylng  the 
general  power  law  f om  for  C  corresponding  to  equation  (B25) ,  there 
results 


(B60) 


I'or  laainar  flow,  n  *  1,  and  all  terns  except  disappear  as 

desired.  For  turbulent  flow,  on  the  other  band,  equation  (B60)  together 
with  the  definitions  of  F^  given  In  equation  (B55),  becowes: 


‘eq,T,0 


>1 

In^l 


(B60a) 


The  tern  In  the  brackets  differs  from  the  previous  definition  of  x^g  ^ 
given  In  (B53)  only  by  the  factor  ' 

^eq.T.O 
^eq,L,0 

arising  fron  the  use  of  a  nonlnal  value  of  n.p  b  4  In  the  definition  of  F,^. 
For  all  cases  of  practical  Interest,  this  tern  will  have  a  negligible 
effect  on  heat  transfer  •*  on  the  order  of  one  percent  or  less. 

Use  of  a  skin  friction  law  of  the  fom  of  (B58)  In  the  heat  transfer 
equation  (B56)  has  the  effect  of  autonatlcally  Introducing  the  local  value 
of  n  In  (B60a),  so  that  the  proper  oonpresslbllity  effect  on  7^  Is  obtained. 


I7U 


Suimary  of  Method  end  Senple  Celculatlone 


Ihle  aection  euwuiriMa  the  eelouletion  procedure  for  the 
•ethed.  Table  B1  liata  epeelfle  ealuea  of  the  »arioue  parasetera  for 
aeveral  apecial  caaea.  Sinoe  the  calculation  dependa  on  the  functlona 
Meeff  •  exiat  only  aa  the  faired  curvea  of  fi#rurea  B2, 

B3,  and  B4  nuaerical  valuea  are  fivea  la  Tablea  B2,  B3  and  B4  froa  ehieh  the 
reader  can  eonetruct  the  neceasary  plota.  Alao  clven  la  Tablea  BS.  B3 
and  B4  are  alaple  curve  fit  expreaalona  that  are  ahoen  to  agree  cloaely 
with  the  valu-a?,  from  the  faired  curvea.  The  curve-fit  expreaalona  were 
not  diaoovered  until  after  the  bulk  of  the  analyala  waa  complete,  however, 
aad  were  not  uaed  for  any  of  the  comparlaona  preasnted  la  thla  report. 

It  la  aaaumed  la  the  following  dlacuaalon  that  the  following  uuaati- 
tlea  are  known: 


^6*  ^8*  ^e*  'vf 

"^S'*  "^0*  Pe*  Mw*  Ms*  “  ^  (^8* 

Pe*  Pw  PS*  =  f  Os* 

The  baalc  equatlona  are  Independent  of  unite,  ao  that  any  conalatent 
aet  of  unita  deaired  can  be  uaed. 

The  baalc  equatlona  to  be  aolved  for  each  caae  la  general  fora  are 


(Bel) 


where 


Cf,r,L 


.00 


(Bdi) 


(B63) 


(B64) 


(B65) 


(B60) 


(B67) 


(B68) 
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G,  =»  p  u  u 

L  "r  e 

(D09) 

Gt  =  Pr  Ur  Ua 

(B70) 

•'l-  [i  * 

(B71) 

Jt  -  [l  t  .49  r, 

(B72) 

-  ['  +  J’o 

(B73) 

(B74) 

Et.i)  -  [i  *  ■■’■>  ro'/^;]o^®’^ 

•;B76) 

Of  =  ,090  /  "-- 

V  Pw 

(B76) 
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If  only  luilnar  flow  is  r«qulr«d  (e.g<,  stagaatlOD  point) 


H 


(B7T) 


Th*  ganaral  procedure  for  «  given  cnee  is  na  follcernt 

a.  IxaBlne  equations  to  reduce  then  to  ninlnua  fora  oonsintent 
«lth  the  given  problem.  For  example,  the  flat  plate  values  of 

are  all  Independent  of  x  and  J  -  1,  so  that  x^q  >  x  in  all  cases. 

b.  Find  Pr  Mr  using  figures  B2  and  b3  or  Tables  BB  and  B3.  Recently 

the  folloiring  expresalons  have  bean  found  for  (Pq  Pr  ' 


(Pe  Me)eff 


Pb* 


,  1. 


86  -  ,  8S 


Pe*  ^8* 


(BTi) 


PrMr 


<Pe  Me)eff 


(PeMe)effl 
(Pw  J 


(B79) 


As  may  be  seen  from  the  values  tabulated  in  Tables  B2  and  B3, 

•q.  (B78)  and  (B70)  agree  with  tbs  curves  of  Figures  B2  and  B3  to  within 
about  3%,  corresponding  to  a  3  to  3  percent  error  la  the  predicted  heating 
rate,  and  so  are  considered  adequate  for  most  purposes. 


0.  find 


1. 

(Z  T)r 

from 

PrMr. 

Pe 

a. 

*r 

from 

<Z  T)*., 

Pe 

3. 

from 

(ZT)r 

4. 

ho 

from 

•s* 

(Z  T), 

and  equation  (B4B) 
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(Not*  that  la  raqulrcd  only  for  turbulent  flow  celculatione, ) 

Any  gas  properties  my  be  UMd  In  these  cslculetlons.  In  the 
present  report  the  gas  properties  of  references  2  and  26  were  used. 

d.  Deternlne  the  streanllne  dlveigence  pcremters  r  and  f.  For 
arbitrary  bodies  at  angle  of  attack  these  paraneters  are  often 
not  kncam,  although  their  product  A  my  '.ie  estlmted  frou 
the  pressure  distribution  or  oil  flow  patterns. 

The  values  of  r  and  f  for  several  special  cases  are  glwen  in 
Table  Bl. 

e.  Find  If  required.  Note  that  for  the  general  case,  exact 

evaluation  of  P^  reqUArea  an  iterative  solution,  since 


fis 


*6 


»eq.L[  “c) 

*  J  d(ln  x) 


(B80) 


and  the  tern  In  brackets  is  In  Itself  a  function  of  Pg  In  practice 
however,  a  finite  difference  Integration  along  the  streanllne  Is  perfomed 
for  the  bracket  paramter 


(B81) 


aaJ  the  local  value  of  Pg  at  x  »  nAx  can  be  evaluated  with  sufficient 
accuracy  (if  the  step  slxe  Is  smll)  by 


r*.] 


1  Q 

ngdanu^)!  r  x^  j 

-'x-nAx 

Lie  d(ln  x)4^  Ll'  X  J 

*  -*x»(n-l)Ax 


(082) 
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f.  Find  Eg,  Eq,  and  using  equations  (B37,  B38,  Fi39,  and 

B40).  A  discussed  In  connection  with  equation  (B50) ,  the  subscript  0 
denotes  evaluation  for  xero  Mach  number. 

g.  Find  r  from  £  using  figure  B4  or  Table  B4.  Recentl;/,  the 
following  expression  was  found: 


r  = 


(B83) 


The  error  in  heating  resulting  from  the  use  of  eq>  (B83)  rather  than  figure  B4 
la  less  than  1%  for  fi  <  10. 

h.  Evaluate  J,  E,  x^q  and  F^,  and  from  the  definitions  given 

earlier  in  this  section. 

I.  Find  lD,e,  lo.e/^e  Pe 

J.  The  heat  transfer  coefficient  H  can  now  be  calculated  from 
equations  (B61)  through  (B63).  Specific  values  of  the  various  boundary 
layer  parameters  are  given  In  Table  B1 . 
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TABLE  B2 


TABULATION  OF  VALUES  OF  p  ^ 


‘Oef'e'off 

P  P 

P  P 

w  w 

From  figure  (B2l 

From  equation  (B79) 

.15 

.226 

.226 

.20 

.290 

.296 

.25 

.354 

.362 

.30 

.414 

.426 

.35 

.485 

.486 

.40 

•  5*r-2 

.  544 

.5 

.641 

.650 

.6 

.739 

.744 

.7 

.818 

.826 

.8 

.892 

» 896 

.9 

.950 

.954 

1.0 

1.000 

1.1 

1.04 

equation  (B79) 

1.2 

1.08 

not  valid  for 

1.3 

1.12 

hot  wall  cases 

1.4 

1.15 

1.5 

1.18 

1.6 

1.21 

1.7 

1.25 

1.8 

i.26 

OCOO^C/I  I  loCOQOMJ^CJt 


TABLE  B3 


TABuLAriOir  OF  VALUES  OF  (n  u  ) 

eff 


Ps’^S’ 

P  P 
e  e 


.2 

.3 

.4 


<Ve>eff 

“A 


From  I’lgure  B3(a) 


From  equation  (B78) 


.335 

336 

.478 

.478 

.610 

.604 

.715 

.713 

.701 

.804 

.850 

.879 

.905 

.936 

.955 

.976 

1.000 

1.000 

TABLE  B4 

TABULATKW  OF  VALUES  OF  V 


2 

r 

From  figure  B4 

From  equation  (B83) 

.5 

.164 

.155 

.6 

.232 

.224 

.8 

.360 

,349 

1.0 

.475 

.460 

1.5 

.710 

.700 

2.0 

.925 

.903 

3.0 

1.27 

1.26 

4.0 

1.56 

1.56 

5.0 

1.83 

1.82 

6.0 

2.07 

2.07 

7.0 

2.30 

2.30 

8.0 

2.50 

2.51 

9.0 

2.70 

2.72 

10.0 

2.87 

2.90 

l83 


REFERENCE  43  36  37 

PRANDT.  NUMBER  CONSTANT  CONSTANT  VARIABLE 

SPECIFIC  HEAT  CONSTANT  CONSTANT  VARIABLE 

p  p  PRODUCT  CONSTANT  VARIABLE  VARIABLE 


Figure  Bl;  EFECT  OF  PRANDTLNUMBER  ON 
REYNOLDS  ANALOGY  FACTOR 


! 
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Figure  Bh  CORRELATION  OF  REFERENCE  DENSITY-VISCOSITY  PRODUCT  (M  -  0) 


(a)  EFFECTIVE  DENSITY-VISCOSITY  PRODUCT  AS  A  FUNCTION  OF  SOUNDARY 
LAYER  EDGE  AND  STAGNATION  ENTHALPY  CONDITIONS^ 


(b)  REFERENCE  DENSITY-VISCOSITY  PRODUCT  AS  A  FUNCTION  OF  EFFECTIVE 
EDGE  VALUES  AND  WALL  CONDITIONS.  M>0. 


Figure  B3:  DENS  ITY-V I SCOSIH’ PRODUCT  CORRELATIONS 

Jfl6 


2 

Figure  B4:  PRESSURE  GRADIENT  EFFECT  CORRELATION  PARAMETERS 
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Figure  B6:  COMPARISON  OF  INCOMPRESSIBLE 

TURBULENT  SKIN-FRICTION  FORMUUS 
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APPHODC  0* 


VOHSIKILAIt  BOUHDART  lAYXa  PRCXHtAM** 


Th«  lurpose  of  the  Konaiatlar  Boundary  Layer  Frograa  la  to  integrate  the 
laainar  boundary  layer  partial  differential  equations  'tHing  finite  difference 
■ethods,  but  without  the  use  of  sJjillarity  aasuaptiona. 

Nearly  all  published  exact  laadnar  boundary  layer  solutions  have  been  ob¬ 
tained  using  the  concept  of  sinllarlty.  These  solutions,  which  aust  be  obtained 
numerically,  require  that  the  viscous  flow  partial  differential  equations  be 
transfonsed  to  a  set  of  ordinary,  non-linear  differential  equations.  In  the 
transforaed  system  the  flow  properties  are  expressed  as  functions  of  a  single 
similarity  variable,  and  are  therefore  independent  of  chordwlse  location.  Un¬ 
fortunately,  the  necessary  transformation  reoulres  certain  flow  ei^itlons 
which  are  rarely  realised  on  realistic  configurations.  The  Boeing  Nonsimilar 
Program,  was  developed  to  avoid  such  limitations. 

The  Nonsimilar  Program  can  calculate  either  stagnation  or  non-stagnation 
houedary  layers  with  arbitral^  pressure  gradients,  with  or  without  mass  In¬ 
jection.  Three-dimensional  flow  effects  arm  calculated  using  the  tero  cross¬ 
flow  approximation,  {/ f[  <*  1,  which  implies  no  rotation  of  the  velocity  vsotora 
within  the  bounda:^  layer.  The  program  is  also  lijaltmd  to  attached  flow. 

The  program  described  herein  treats  air  in  ehemioal  equilibrium.  The 
program  can  be  applied  to  ideal  ipis  and  other  fluids  by  changing  tha  tabulated 
gas  transport  property  tables. 

The  program  is  capable  of  initiatiag  its  own  boundary  layer  solutions, 
given  only  external  flow  properties,  for  either  the  stagnation  point  or  sharp 
tip  cones  and  plates. 

Basic  Iquations.-  The  equations  solved  by  this  computer  progremi  are  the 
standard  boundary  layer  equations  of  state,  continuity,  x-mosmntun  and  energy. 
These  equations  are  presented  in  a  form  used  by  the  program  for  evaluation  at 
several  vertical  positions  in  the  boundary  layer  at  each  value  of  x  considered. 


*  This  Appendix  is  baaed  on  Appendix  C  of  Reference  l8. 


*#  Tl»l8  computer  program  was  developed  by  A.  L.  Nagel  and  R.  T.  Savage 
during  the  X-20  development  program. 
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EQJATIOH  OF  STATS 

P 

P  •  R  (*T) 

ooMTDrom 


(Cl) 


vbar*  8ub*oript  i-1  refers  to  y  •  y  -Ay 
x-MONERUN 


lx 


(C3) 
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Equation  (C2)  in  obtained  from  the  continuity  equation  by  Introducing 
equation  (C3)  and  (CU)  to  ellnlnute  )u/9x  and  *H/|X.  No  atomic  diffusion 
terms  are  required  in  equation (CU)becau.-.e  thl  ;  mode  of  energy  transport  has 
been  Included  in  the  Prandtl  number. 

Foivard  Integration.-  Since  v  Is  expressed  as  a  function  of  Input  data, 

V  can  be  determined  explicitly  at  each  point  in  the  boundary  layer  at  the 
Initial  or  start  position.  With  v  defined,  the  Inltlcd  u  and  H  profiles, 
}u/)x  and  dYi/ix  can  be  determined.  With  9u/9x  and  iVi/ix  determined,  the  pro¬ 
files  at  the  next  station  can  be  obtained  by  forward  Integration  using  the 
following  equations 
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At  erAh  point  In  the  boundary  ligrer,  the  elassleal  slnllaxaty  paraaeter 
la  ealenlated  using  the  following  relationship. 


Also  aalottlated  at  each  station  (x)  are  the  boundary  layer  displaosaent 
thlokness^S*,  heating  rate,  q^,and  shear  at  the  vall,r^. 


9he  heating  rate  and  shear  at  t^e  wall  are  oaleulated  using  the  energy 
and  wMsntus  integral  equations  rather  than  the  definitions  beoause  of  the 
^[reater  aoeuraoy  obtainable  with  the  Integrals.  Ae  use  of  the  definition  of 
q^  and  requires  yery  snail  Inoreasnts  to  obtain  aeeurate  yalues  of  the 

gradients  at  the  wall . 

For  problens  without  yortioity,  ■  0.  For  eases  with  yortloity, 
(|u/}y)^  is  Input  as  a  funetlon  of  x* 
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stagnation  Region  Caleolntloc 

Stagnation  point  profiles  are  obtained  by  nodlfylng  the  proeedure  Just 
discussed.  Input  profiles  are  corrected  by  Integrating  the  u  and  H  profiles, 
but  X  (not  equal  to  zero)  Is  not  Increased  during  the  Integration,  lihe  pro* 
files  are  assuned  correct  i^n 


<  .  1  Ug  for  all  1  (Clll 


irtiere  Is  the  Telocity  at  the  edge  of  the  Input  Telocity  profile  (l.e.  the 

Input  Telocity  farest  fron  the  vail.)  This  convergence  criterion  is  obtained 
vlth  the  velocity  slallarlty  stated  belov. 


lH| 

#x|e  ■  u  lx 


(C12) 


Experience  has  Indicated  that  an  enthalpy  profile  convergence  criterion  Is 
not  necessary. 

For  the  stagnation  region,  equatlcns  (C^)  and  (C6)  are  changed  to 


idiere  x^  *  Initial  x  location 
o 

AS- 


(C13) 


{Ci4) 


19>» 


a  fictitious  length 


Onoe  the  convergence  cilterle  have  bee:2  aatlsfied,  the  calculations  can 
be  stepped  forward  around  the  body  as  discussed  in  the  previous  section. 

Qas  Properties.-  As  stated  previously  the  progrea  treats  air  in  cheaical 
equllibriua  and  can  be  applied  to  air  as  an  ideal  gas  and  to  other  fluids  by 
changing  the  tabulated  gas  transport  property  tables. 

The  transport  properties  for  equlUbrlum  air  were  based  on  the  nine-species 
aodel  (Hg,  MO,  H,  0,  il+,  0+,  and  c  )  of  reference  2  and  conputed  using  the 

collision  integral  method  of  reference  W*.  The  transport  properties  are  built 
into  the  program  as  tabular  functions  of  enthalpy  and  pressure. 


T  = 


(C15) 


t 


H 


(C16) 


H 


(C17) 


(C18) 


195 


APPENDIX  D 


FORMULATION  OF  h8* 


The  boundary  layer  dlsplAcetncnt  thickness  heating  parameter,  defined  by 


(Dl) 


is  required  for  evalmting  the  effects  of  shallow  surface  waves  on  aerodynamic 
heating.  This  appendix  describes  the  methods  that  were  used  to  compute  the 
values  shown  on  Figure  6-45.  Only  laminar  flow  has  been  considered. 


Methods  for  estimating  the  laminar  displacement  thickness  were  developed  by 
Honks  and  Savage,  Reference  29»  by  correlating  the  boundary  layer  solutions 
first  presented  by  Cohen,  Reference  37.  The  closed  form  equation  for  6* 
that  is  presented  in  Reference  29  la  rewritten  here  and  forms  the  basis  for 
the  development  of  the  equations  used  to  define  the  product,  h6*.  For  zero 
pressure  gradient 


e 

where 


.* 

L  • 


-0,905  2  (TZ)^j  ♦ 

tm 


(TZ) 


ref 


0.73 


(n). 


(TZ) 

♦  - i 

(TZ), 


Results  shown  on  Figure  6-45  were  calculated  for  a  sharp  flat  plate.  For 
this  geometry 


r  •  1 


•  const 
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APPENDK  E 


SAMPLE  SOLUTIONS 

To  illustrate  the  use  of  the  charts  in  Sections  VI  and  VII,  example 
problems  were  prepared.  In  most  instances,  examples  have  been  chosen  that 
clearly  explain  the  procedure  required  to  reduce  the  curves  to  correct 
numerical  results.  The  applicable  text  section  number  precedes  each  example. 

A.  Section  VI  -  Reference  Condition 

Radius  Correction  -  Turbulent  Reference  Heat  Transfer  Coefficiert 

Example:  Determine  the  turbulent  stagnation  line  heat  transfer  coeffi¬ 
cient,  hgj  for  a  60®  swept  cylinder  having  a  leading  edge  radius  of 

R  =  0.1  ft.  The  flight  condition  is  VEL  =  10,000  fps  and  ALT  =  100,000  ft. 


Quantity 

Source 

Value 

Figure  6-4 

1.6  X  lO"^  Btu/ft^-sec-®R 

© 

Figure  6-7 

360,000 

© 

P 

^RT 

Figure  6-6 

1.64 

^SL,T 

© 

© 

2.62  X  lO"^^  Btu/ft^-sec-®I 

B.  Section  VI  -  Basic  Shapes 


1.  Sharp  Flat  Plate  and  Sharp  Cone 

The  computational  procedure  required  to  obtain  heating  rates  for  a  sharp 
flat  plate  or  a  sharp  unyawed  cone  is  similar.  The  procedure  is  illusti'ated 
using  a  sharp  flat  plate. 

a)  Laminar  Flow 

Fjcampie:  Detennlne  the  laminar  heat  transfer  rate  to  a  sharp  flat  plate 
at  an  angle  of  attack  of  20®,  at  a  distance  x  *  5  ft,  for  T„  =  1500  ®R  and 
Ta^/T^  =  0.87.  The  flight  condition  is  VEL  =  20,TCC  fps  and  ALT  =  200,000  ft. 


Quantity 

Source 

Value 

Figure  6-23 

0.275  (Altitude  effects 
neglected  for 
figui'e  6-23) 

® 

h 

0 

Figure  6-2 

3.56  X  10"3  Btu/ft^-sec-' 
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Vi. 

(D  X  © 

9.77  X  10"^  Btu/ft^-sec-®R 

© 

c 

X 

Figure  6-25 

0.4h5 

© 

\  =  5  f  t 

0 

© 

4.35  X  10'^  Btu/ft^-sec-“R 

© 

^T 

Figure  6-2^ 

:.o 

@ 

=  5  ft 

©  X  © 

4.35  X  10"^  Btu/ft^-sec-°R 

® 

T 

^AW 

Figure  6-1 

31,500  *R 

@ 

4 

=  @©-  1500) 

13.05  Btu/ft^-sec 

b) 

Turbulent  Flow 

To  obtain  turbulent  heating  rates  to  a  sharp  flat  plate  or  a  sharp  un¬ 
yawed  cone,  use  the  turbulent  reference  and  design  charts.  The  procedxire  is 
identical  to  that  shown  above. 

2. 

Hemisphere 

a) 

Laminar  Flow 

Ebcample.  Determine  the  laminar  heat  transfer  rate  to  a  1-foot  radius 
hemisphere  at  the  angular  location  B  =  60* .  The  fli^t  condition  is  VEL  * 
20,700  fps  and  ALT  =  200,000  ft.  Assume  *=  2000  *R  and  T^^/T^  =  0.95. 

Quantity 

Source 

Value 

Figure  6-9 

0.395  (Altitude  effects 
neglected  for 
figure  6-9) 

© 

h 

0 

Figure  6-2 

3.56  X  10"^  Btu/ft^-sec-*R 

© 

b  „  ^  „ 

60  * 

®  X  © 

1.4l  X  10"'^  Btu/ft^-sec-*R 

© 

^Aw/’^T 

Assume 

0.95 

© 

^AW 

Figure  6-1 

34,500  “R 

© 

• 

9 

-©  (©-  2000) 

Wall  temperaturf  corrections  can  be  made  using  figure  6-3* 
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b) 


Turbulent  Flow 


I 


Example:  Determine  the  turbulent  heat  transfer  rate  to  a  hemisphere 
using  the  conditions  given  in  the  example  for  lami.iar  flow. 


Quantity 

Source 

Value 

R 

’^MAX 

Figure  6-10 

1420 

© 

^MAX 

Figure  6-11 

36.25* 

'  ^MAX 

+23.75* 

Figure  6-12 

0.715 

© 

^^AX'^^SP 

Figure  6-13 

0.66 

© 

h 

o 

Figure  6-2 

3.56  X  10"^  Btu/ft^-sec-"R 

h 

© 

0 

X 

0 

1.53  X  10’^  Btu/ft^-3ec-*R 

^Aw/^T 

Assume 

0.95 

^AW 

Figure  6-1 

35,000  "R 

• 

9 

h(T^„  ’ 

=  ©  (  ©  -  2000) 

50.5  Btu/ft^-sec 

When  the  radius  of  the  hemisphere  is  different  from  R  =  1  ft,  tnc  value  of 
R  must  be  corrected  using  equation  (6.10)  before  using  figures  6-11  and 
^MAX 
6-13. 

3.  Swept  Infinite  Cylinder 
a)  Laminar  Flow 


Example:  Find  the  laminar  heat- transfer  rate  to  a  40®  swept  cylinder, 

1  foot  in  radius  at  an  angular  location,  60® .  The  flight  condition  is 

m  ==  20,700  fps  and  ALT  =  200,000  ft.  Assume  T^  =  1000  *R  and  =  0.94. 


Quantity 

Source 

Value 

© 

Figure  6-17 

0.222 

© 

h 

0 

Figure  6-2 

3.56  X  10  ^  Btu/ft^-sec-®R 

® 

®  X  © 

0.79  X  10‘3  Btu/ft^-sec-'R 
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0000000©  gH>  ^  ©©©  ©0© 


M 

“n 

Mjq  X  cos  Uo* 

17 

"l 

Figure  6-3 

1.05 

^  ^  =  60“ 

© 

X 

0 

0.83  X  10  ^  Bt.u/ft^-sec-“R 

T^  =  1000  “R 
^Aw/'^T 

Assume 

0.94 

T 

^AW 

Figure  6-1 

34,000  “R 

• 

^^^AW  ■ 

=  ©  (  ®-  1000) 

27*4  Btu/ft^-sec 

b)  Turbulent  Flow 


The  procedure  for  determining  turbulent  heating  rates  to  swept  cylinders 
similar  to  that  shown  above  for  laminar  flow. 

4.  Unswept  Infinite  Cylinder 

a)  Turbulent  Flow 


Example;  Determine  the  turbulent  heat  transfer  coefficient  to  a  1 
ot  radius  unswept  cylinder  at  an  angular  location,  0  =  60* .  The  flight 


ndltion  is  VEL  =  18,7000  fps  and  ALT  =  l80,000  ft. 


^S3vme  =  2000  ®R. 

w 


Quantity 


Source 


Value 


R 

^MAX 

^MAX 

Ad 

^lAx/^SP 

h 

o 

’‘^SP 

^  e=  6o“/^4ax 
^  60“ 


Figure  6-20 

5,500 

Figure  6-21 

33.6“ 

■  ^MAX 

+26.4“ 

Figure  6-13 

0.7 

Figure  6-2 

4.55  X  10"^  Btu/ft^-3ec"“R 

^SP  \(R  =  1  ft) 

4,-5  X  10  ^  Btu/ft^-sec-“R 

Figure  6-22 

0.66  Btu/ft^-sec-“R 

0x©x© 

2.10  X  10’^  Btu/ft^-sec-“R 

NOTE: 

When  the  radius  is  different  from  R  =  1  ft,  R  must  be  corrected  using 

MAX 
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equation  (6.10)  before  using  figures  6-13  and  6-21. 
4.  Sharp  Delta  Wing 
a)  Laminar  Flow 


Example:  Determine  the  laminar 

sharp  delta  wing  at  x  =  5  ft.  on  the 

centerline  on  the  rayline  e  /  e  = 

MAa 

of  20“  at  the  flight  condition  VEL  = 


heat  transfer  coefficient  to  a  75*  swept 
centerline ;  and  at  x  =  5  ft.  off  the 
0.5.  The  wing  is  at  an  angle  of  attack 

20/700  fps  and  ALT  =  240,000  ft. 


Quantity 

Source 

Value 

h-p,  ,/h 

Figure  6-30 

0.325 

© 

h 

0 

Figure  6-2 

1.62  X  10‘^  Btu/ft^-sec-“R 

© 

^SD,L 

@  X  © 

••li  p 

5.26  X  10  Btu/ft  -sec-“R 

Figure  6-25 

0.445 

© 

^x  =  5  ft 

© 

© 

0 

2.34  X  10  Btu/ft  -sec-“R 

© 

70“  ^x  =  5  ft 

Figure  6-3I 

1.065 

© 

^x  =  5  ft 

©  X  © 

2.50  X  10"^  Btu/ft^-sec-“R 

A'  75“ 

© 

II 

Figure  6-32 

1.01 

^x  =  5  ft 

©  X  ® 

-4  /  2 

2.52  X  10  Btu/ft  -sec-* R 

The  wall  temperature  correction  is  obtained  from  figure  6-24. 


To  determine  turbulent  heat  transfer  coefficients  to  a  sharp  delta  wing 
requires  the  use  of  turbulent  reference  and  design  charts.  The  procedure  is 
identical  to  that  required  to  obtain  laminar  heat  transfer  coefficients. 


C.  Section  VI  -  Combinations  of  Basic  Shapes 


1.  Unyawed  Blunt  Cone 
a)  Laminar  Flow 


Example:  Calculate  the  laminar  heat  transfer  rate  to  a  blunted  cone 
having  a  hemispherical  nose  of  R  0.5  ft  and  a  serai -vertex  angle  of  30“  at 
S/R  =  4.  The  flight  condition  is  VEL  =  20,700  fps,  and  ALT  =  200,000  ft. 
Assume  T,^  =  2000  “R  and  T^^^/T^  =0.95. 
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©  ©00  ©©©©©©  g  0©©©©©© 


Quantity 

Source 

Value 

^^,l/^sp^s/r=  4 

Figure  6-27 

0.198 

h 

0 

Figure  6-2 

Ip 

3.56  X  10"^  Btu/ft^-sec-“R 

^^c,l^s/r  =  4 

-  "ot 

5.04  X  10‘^  Btu/ft^-sec-*R 

^S/R 

®  X  @ 

.995  X  10“^  Btu/ft^-sec-*R 

^Aw/^T 

Assume 

0.95 

'’’^AW 

Figure  6-1 

34,500  “R 

• 

<1 

^^^AW  ' 

*©(©-  2000) 

32.3  Btu/ft^-sec 

b)  Turbvilent  Flow 

Example:  Determine 

the  turbulent  heat  transfer  rate  to  an  unyewed  blunt 

using  the  conditions 

shown  in  the  example  above. 

Quantity 

Source 

Value 

^  R  =  1  ft 

Figure  6-10 

1420 

*MAX 

=  0.5  ft 

n® 

710 

)  s/r=4 

^  ^MAX 

Figure  6-28 

8.5 

R 

r 

©  X  @ 

0O35 

Figure  6-29 

1.23 

\ 

From  Step  ©  in 
laminar  heating 
example 

0.995  X  10"^  Btu/i't"-sec-" 

s 

©  X  © 

1.193  X  10”^  Btu/ft^-sec-* 

Assune 

0.95 

^AW 

Figure  6-1 

34,500  "R 

• 

<1 

"(’’ab  -  V 
-  ©  (  @  -  2000) 

38.8  Btu/ft^-sec 
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2 .  Blunt  Delta  Wing 
a)  Laminar  Flow 

Example:  Determine  the  laminar  heat  transfer  coefficient  to  a  70*  swept 
blunt  delta  wing  at  the  centerline  location  s/r  =  IJ.  Also  determine  the 
leading  edge  stagnation  line  value,  and  its  amount  of  shift  due  to  lower 
surface  effects,  as  well  as  the  heating  on  the  lower  surface  at  N/R  =  4.  The 


flight  condition  is  VEL  »  20,700  fps  and  ALT 

=  200,000  ft 

at  an  angle  of 

attack  of  30“ .  The  nose  and 

leading  edge  radii  are  1  ft. 

Quantity 

SouiMe 

Value 

^^d,l/^sp^s/r=17 

Figure  6-34 

0.101 

© 

h 

o 

Figure  6-2 

3.56  X 

10‘^  Btu/ft^-sec-”R 

® 

^SP 

^SP,(i^l  ft)  "" 

h^  3*56  X 

10"^  Btu/ft^-sec-“R 

^BD,L 

©  X  © 

3.59  X 

lO"^  Btu/ft^-sec-“R 

© 

^SL 

Figure  6-38 

65“ 

© 

■^EFF 

Figure  6-39 

54.5“ 

^^sl/^o^Aeff 

Figure  6-16 

0.362 

® 

^■SL,A  eff 

®x@ 

1.29  X 

lO""^  Btu/ft^-sec-“R 

@ 

^^d,l'^^sp^n/r=4 

Figure  6-36 

0.126 

© 

®x® 

4.49  X 

10  Btu/ft  -sec-“R 

b)  Turbulent  Flow 

Example:  Solve  for  the 

turbulent  heat 

transfer  coefficients  on  a  bl\uit 

delta  wing  using  the  conditions  given  in  the  example  for  laminar  flow  except 
let  the  flight  condition  change  to  VfiL  =  18,700  fps  and  ALT  =  l80,f'''0  ft. 


Quantity 

Source 

Value 

^^d,l/^sp^s/r  =  17 

Figure  6-34 

0.101 

© 

h 

0 

Figure  6-2 

4.55  X 

10"^  Btu/ft^ -sec -“R 

® 

^SP 

*^SP  (R=1  ft)  ""  ^0 

4.55  X 

10‘^  Btu/ft^ -sec -“R 

© 

^D,L 

©X® 

4.59  X 

10"^  Btu/ft^ -sec -*R 
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© 

R 

Figure  6  -10 

2450 

MAX 

© 

Figu.’.*e  6-40 

37 

© 

©X© 

91,000 

^  s/r  =  17 

Figure  6-29 

2.6 

©X© 

1.19  X  10"^  Btu/ft^-sec-*R 

^SL 

Figure  6-38 

65® 

© 

^^EFF 

Figure  6-39 

54.5* 

© 

Figure  6-18 

1.22 

@ 

Figure  6-4 

2.82  X  10”^  Btu/ft^-sec-®R 

® 

eff 

@x@ 

3.44  X  10“^  Btu/ft^-sec-®R 

NOTE: 


When  the  nose  radius  is  different  frcan  R  =  1  ft,  the  value  of 
corrected  by  equation  (6.10)  before  proceeding  with  step  (^. 


must  be 


At  the  present,  no  method  for  evaluating  turbulent  heating  between  the 
centerline  and  the  leading  edge  is  available.  Estimates  must  be  made  by 
fairing  a  curve  through  the  end  points. 


3.  Flap  Surfaces 
a)  Laminar  Flow 


The  increase  in  laminar  heat  transfer  across  a  flap  can  be  directly 
obtained  using  figure  6-4l. 

b)  Turbulent  Flow 

Example:  Determine  the  increase  in  turbulent  heating  across  a  10® 
positively  deflected  flap  vrhere  the  Initial  surface  is  at  an  angle  of 
attack  of  20®  for  the  fligtit  condition  VEL  =  18,700  fps  and  ALT  =  l80,0C0  ft. 


Quantity 

Source 

Value 

Figure  6-4l 

3.3 

Figure  6-42 

0.925 
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@  (V^^T  ®*®  3.05 

D.  Section  VI  -  Surface  Condition  Effects 
Surface  Roughness 

Example!  Determine  the  Increase  In  laminar  surface  heating  on  a  flat 
plate  at  an  angle  of  attack  OL=  20®,  for  the  fli^t  condition  VEL  =  20,700 
fps  and  ALT  =  2^0,000  ft.  Assume  the  wave  hei^t  parameter  W/r  =  25  and 


R  = 

0.1  ft. 

Quantity 

Source 

Value 

Q 

h  8* 

0 

Figure  6-^5 

1.26  X  10  ^  Btu/ft-3ec-®R 

© 

h/h 

Figure  6-23 

0.275 

© 

h 

0 

Figure  6-2 

1.62  X  lO"^  Btu/ft^-sec-®R 

8* 

®  /((i)x®) 

0.0283  ft 

R/8* 

0 

• 

® 

3.53 

© 

Figure  6-^3 

>1.7 

E. 

Section  VI  -  Application  To  Complex  Configurations 

Control  Surface  Gaps 

a)  Laminar  Flow 

Example:  Determine  the  maximum  heat-tremsfer  coefficient  on  the  lower 
surface  fin  and  control-surface  gap  for  the  fli§^t  condition  VEL  *  20,700  fps 
and  ALT  =  2^0,000  ft.  The  vehlc) ?  is  at  an  emgle  of  attack,  Of=  20® .  The 
surfaces  creating  the  gap  have  a  circular  cross-section  of  radius  =  1  ft  and 
are  far  enough  apart  to  prevent  choking. 

Quantity 

Source 

Value 

© 

**0,^ 

Figure  6-^7 

0.0332 

© 

h 

0 

Figure  6-2 

1.62  X  10’^  Btu/ft^-sec-®R 

© 

^G,L 

®x© 

0.538  X  10”^  Btu/ft^-sec-“R 

b)  Turbulent  Flow 


The  solution  for  turbulent  b.at-transfer  coefficients  in  a  control 
surface  gap  requires  the  iice  of  tiu:bu3.eat  refsrence  and  design  charts.  The 
procedure  is  identical  to  that  for  laminar  flow. 


207 


F.  Section  VII  -  Use  of  Ground-Test  Data  in  Design 

Two  methods  are  used  to  extrapolate  wind-tunnel  data  to  fli^t  conditions. 
The  first  method  involves  the  use  of  canpressibility  parameter  charts  and 
extrapolation  factor  charts.  Flat  plates,  sharp  cones.  Infinite  swept 
cylinders,  sharp  delta  wings,  and  control  surface  gaps  fall  under  this 
method  of  extrapolation.  The  procedure  in  cwnputing  the  results  is  the  same 
in  all  cases  and  thereftjre  only  one  example  is  provided  that  is  applicable  to 
all  of  the  geometries  mentioned  above.  This  procedure  is  illustrated  in  the 
examples  under  the  subheading,  "General  Method."  The  second  method  is 
illustrated  in  the  examples  under  the  subheading  "Method  for  Deflected  Flaps." 
The  extrapolation  from  wind  tunnel  to  fli^t  for  laminar  flow  on  a  hemisphere 
is  a  special  case  (A  =  l)  and  can  be  obtained  directly  from  figure  7“25- 

1 .  General  Method 
a)  Laminar  Flow 

Example:  Laminar  data  is  available  from  a  wind  tunnel  on  a  fiat  plate 
at  an  angle  of  attack,  a=  20“,  at  x  =  o.4  ft.  When  normalized  to  the 
stagnation  point  heat- transfer  coefficient  for  a  1  ft  radius  sphere  the  ratio 
is  0.38.  The  total  enthalpy  of  the  turn  .el  is  2  x  10'  ft  /sec'^  and  the  test 
section  Mach  number  is  32.  Determine  the  heating  rate  ratio  for  the  laminar 
design  flight  condition,  ’TEL  =  20,700  fps  and  ALT  =  240,000  ft. 


Quantity 

Source 

Value 

Q 

^^,L^^oW  - 

’  X  c  1  ft 

0.38  X  •- 

0.24 

© 

S-LT 

i  ~V^/2 

FLT 

2  X  10® 

@ 

(  2  X  10^) /(  2  X  10®) 

0.1 

® 

1 

Figure  7“19 

0.854 

© 

^^,L^^o^FIT  A  fP,L 
^FP,L 

Figure 

1.079 

© 

^^,L'^^0^FLT 

®  /(0  X  ©) 

0.261 

© 

^^x^^o^FLT 

/TN  1  X  =  1  f  t 
^  yx  =  .4  ft 

0.413 

b )  Turbulent  Flow 


To  detemine  turbulent  heating  rate  ratios  use  the  turbulent  reference, 
design,  and  extrapolation  charts.  The  procedure  is  identical  to  that  for 
laminar  flow. 


208 


?.  Method  for  Deflected  Flaps 
a)  Laminar  Flow 

Example:  Laminar  heating  data  from  a  Mo,=  10  wind  tunnel  show  a  100$ 
increase  in  heating  for  a  10*  positively  deflected  flap  whose  initial 
surface  is  Reflected  20®,  i.e.,  angle  of  attack,  (X=  20*.  Compute  the 
expected  change  in  heating  for  the  same  geometry  at  the  laminar  reference 
condition,  VEL  =  20,700  fps  and  ALT  =  2h0,000  ft.  If  data  show  a  20f>  decrease 
in  heating  in  the  wind  tunnel  for  a  10®  negatively  deflected  flap,  compute 
the  expected  change  for  the  flight  condition  stated  above. 


Quantity 

Source 

For 

6=  +10*.  a=  20* 

® 

^^'^^IL^FLT 

Figure  7-3O 

1.29 

© 

Given 

2.0 

© 

^V^n^FLT 

0 

X 

® 

2.58 

For 

5=  -10*.  a=  20* 

® 

^^'^^IL^FLT 

Figure  7-3O 

0.82 

© 

Given 

0.8 

© 

^V^\tJfLT 

© 

X 

© 

0.656 

b)  Turbulent  Flow 


For  turbulent  flow  the  procedure  is  the  same  as  shown  above  except  an 
additional  step  is  required  to  convert  to  expected  turbulent  values.  This  is 
dene  by  using  figure  6~k2  in  conjunction  with  the  calcr.iated  laminar  value. 
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